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Presentazione del lavoro di 
questa tesi 



La spettroscopia e lo studio della materia e delle sue proprieta realizzato 
analizzando la luce, il suono, o le particelle che vengono emesse, assorbite 
o riflesse dalla materia che si vuole descrivere. La spettroscopia puo anche 
essere definita come lo studio dell'interazione tra luce e materia [2TJ. Storica- 
mente la spettroscopia era riferita a un tipo di scienza in cui la luce visibile 
veniva usata per studi teorici sulla struttura della materia e per analisi quali- 
tative e quantitative. Recentemente, comunque, la definizione ha assunto un 
carattere piu generale con lo sviluppo di tecniche sperimentali che utilizzano 
non solo la luce visibile, ma anche molte altre forme di radiazione elettro- 
magnetica e non: microonde, onde radio, raggi X, elettroni, fononi (o onde 
sonore) e altre. Tali tecniche vengono utilizzate nel campo della fisica e del- 
la chimica per l'identificazione di sostanze attraverso lo studio dello spettro 
emesso o assorbito da esse; in astronomia la spettroscopia viene fortemente 
utilizzata e molti telescopi sono dotati di spettrografi utilizzati per misurare 
la composizione chimica e le propieta fisiche di oggetti astronomici. 

Le eccitazioni elettroniche in particolare sono alia base degli spettri di 
eccitazione che vengono piu comunemente studiati. Negli ultimi vent'anni ci 
sono stati numerosi miglioramenti delle tecniche sperimentali utilizzate per 
lo studio delle eccitazioni elettroniche in materia condensata; d'altra parte 
molti degli aspetti fisici legati a tali esperimenti richiedono, per essere de- 
scritti, opportuni modelli teorici. Appare dunque interessante sviluppare tali 
modelli e verificarne la validita utilizzando i dati sperimentali piu recenti. Lo 
scopo e quello di disporre di strumenti in grado di prevedere i risultati degli 
esperimenti prima di eseguirli. Infatti l'approccio teorico offre il vantaggio 
di essere meno dispendioso, sia in termini economici che in termini di risorse 
umane e di tempo. 

In ogni esperimento di spettroscopia un campione viene perturbato e si 
misura la risposta del sistema a tale perturbazione, quindi in generale non 
e sufficiente calcolare le proprieta di stato fondamentale del sistema per in- 
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Pirett Inverse Absorption 

photoenussion phutuemission 




Figura 1: Rappresentazione schematica delle eccitazioni coinvolte in diffe- 
rent tipi di spettroscopia: fotoemissione diretta, fotoemissione inversa e 
assorbimento. 



terpretare e prevedere i risultati degli esperimenti. La foto-emissione diretta 
ed inversa e l'assorbimento possono essere considerati come il prototipo dei 
fenomeni che si vorrebbero descrivere in questo contesto (fig. [1]). II lavoro 
della presente tesi e focalizzato sullo studio dello spettro di assorbimento di 
sistemi isolati che presentino effetti legati alio spin. 

La tesi si apre, Capitolo 1, con un'introduzione alle equazioni che sono 
alia base dello studio di un sistema fisico tramite il formalismo della mecca- 
nica quantistica e vengono illustrate quali siano le maggiori difhcolta che si 
incontrano nel tentativo di risolverle. Sono quindi esposte nei due capitoli 
successivi le due teorie che vengono utilizzate per superare tali difficolta nella 
descrizione dello stato fondamentale del nostro sistema: la Density Functional 
Theory (DFT), Capitolo 2, e la Many Boduy perturbation Theory (MBPT), 
Capitolo 3. Nel Capitolo 4 vengono quindi illustrate le approssimazione in- 
trodotte nell'ambito della DFT, che costituisce la teoria che utilizzeremo per 
i calcoli eseguiti nell'ambito di questa tesi. 

La descrizione dello stato fondamentale del sistema e il punto di partenza 
per poterne descrivere lo spettro di eccitazione ma, come gia detto, non e 
sufheiente a comprendere la fisica che sta alia base degli esperimenti di spet- 
troscopia. Nel Capitolo 5 della tesi vengono quindi esposte le estensioni delle 
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(a) densita dello stato fondamentale (b) densita del primo stato eccitato 



Figura 2: Molecola test del BeH. Raffronto tra la densita dello stato fonda- 
mentale e quella del primo stato eccitato. Le immagini sono state costruite 
a partire dall'interpretazione delle funzioni d'onda di KS. Solo l'immagine 
dello stato fondamentale potrebbe, in base a giustificazione formale, essere 
interpretata come la densita reale del sistema. Possiamo comunque osservare 
come la densita elettronica attorno all'atomo di sinistra, atomo di idrogeno, 
sia diminuita poiche la nube elettronica e ora distribuita su di uno stato ec- 
citato. Quelle rappresentate in figura sono due isosuperfici a valore costante 
della densita; i colori sul piano indicano invece il valore della densita sul piano 
stesso, il rosso indica le zone di massima densita, il blu le zone a densita piu 
bassa. 



due teorie presentate necessarie per lo studio delle energie di eccitazione di 
un sistema ed in particolare l'equazione di Bethe-Salpeter, estensione della 
MBPT, e la Time Dependet Density Functional Theory (TDDFT), estensio- 
ne della DFT. Come per la parte di stato fondamentale le approssimazioni 
necessarie per l'utilizzo della TDDFT vengono esposte nella parte finale del 
capitolo. 

Nella seconda parte della tesi ci concentriamo invece in modo piu det- 
tagliato su quello che e il formalismo di Casida per la TDDFT, esposto nel 
Capitolo 6, formalismo atto alio studio di sistemi isolati e che e in particolare 
legato alia descrizione degli effetti di spin, presenti sia nello studio di sistemi 
a spin perfettamente equilibrato, sia nello studio di sistemi con stato fon- 
damentale spin polarizzato, polarizzazione collineare. Nel capitolo vengono 
in particolare messe in luce le differenze tra i due tipi di sistemi nell'ambito 
dell'approssimazione eseguita. 

II lavoro della presente tesi e stato in particolare mirato alio studio di 



viii 



Presentazione del lavoro di questa tesi 




(a) densita a spin up (b) densita a spin down 

Figura 3: Molecola test del BeH. Raffronto tra la componente a spin up e 
quella a spin down della densita dello stato fondamentale. 



sistemi con stato fondamentale spin-polarizzato. Le modifiche necessarie per 
lo studio di sistemi spin-polarizzati sono state implementate all'interno del 
software ABINIT [2], che in precedenza permetteva di studiare al livello 
TDDFT solo sistemi non polarizzati. E stata poi scelta la molecola di BeH 
(fig. 121 13]) per testare l'implementazione fatta confrontando i risultati da noi 
ottenuti con altri calcoli teorici e con dati sperimentali, Capitolo 7. 



Parte I 
Theoretical background 
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Convenzioni utilizzate nel capitolo 1 

n: numero di elettroni 

N: numero di ioni 

m: massa elettrone 

M: massa ione 

e: carica dell'elettrone 

r: coordinata elettrone 

R: coordinata ione 

p,q: indice elettroni 

P,Q: indice ioni 

io: abbreviazione per ione 

el: abbreviazione per elettrone 

T: energia cinetica 

V: potenziale esterno 

H: Hamiltoniana 

i: v 73 ! 

h: costante di Plank ridotta 
funzione d'onda molti corpi 
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Capitolo 1 

La fisica della materia 



1.1 Lo studio dei sistemi a molti corpi 

Studiare le proprieta di una molecola o di un materiale qualsiasi scegliendo 
come punto di partenza la meccanica quantistica significa cercare di risolvere 
l'equazione che regola il comportamento di tutte le particelle che costituiscono 
tale sistema. Si tratta piu in specifico dell'equazione di Schroedinger per il 
moto accoppiato di atomi ed elettroni, che, definite le seguenti grandezze 

n t2y 2 N *2y2 

el ^ 2m lon 2s 2M 

p=i p=i 

e 2 * e 2 



Vel—el / i i ^ion- 

*s \r v — r *s Rp — R 

n N 2 
Vion—el = ~ / / 7^ r 

*— ' *s Rp — rJ 

p =i p=i i ^ pi 

= tf(Ri,...,Rjv,ri,...,r n ,i) , 



risulta: 



(^T ei + T ion + V;;_ e i + Vi on _ ion + Vion-elJ * = ^"^T ■ (I- 1 ) 

Cosi formulata tale equazione rappresenta un problema di difficolta non 
risolvibile, perlomeno con i mezzi matematici che abbiamo oggi a disposizione. 
Questo per diversi motivi: innanzi tutto perche e un'equazione differenziale 
alle derivate parziali e perche deve descrivere il comportamenti di un numero 
di particelle interagenti che puo essere dell'ordine del numero di Avogadro 



6 



La fisica della materia 



(10 23 particelle!). Di questo tipo di equazione esiste una soluzione esatta, al- 
meno in 3 dimensioni, solo per lo studio di una singola particella nel vuoto o 
immersa in potenziali che presentino particolari simmetrie o infine per lo stu- 
dio di due particelle interagenti. In altre parole per quanto riguarda la fisica 
della materia siamo in grado di trovare una soluzione esatta a questa equazio- 
ne solo per quanto concerne lo studio di un elettrone libero o per l'elettrone 
deH'atomo di idrogeno (l'unico elemento che presenta un solo elettrone). 

L'avvento dei computer ha offerto una possibilita alternativa al calcolo 
analitico per la soluzione del problema. D'altra parte, come mostrato da 
Walter Kohn nelle sue Nobel Lecture [6] anche con il supporto di software 
numerici la possibilita di studiare la funzione d'onda molti corpi resta limitata 
a sistemi che presentano un numero di particelle ridotto. Cercando di studiare 
sistemi piu complessi ci si trova infatti di fronte a quello che Kohn definisce 
un "Exponential Wall" (Muro Esponenziale) e che di fatto limita le attuali 
possibilita, nel caso di soluzione diretta dell'equazione, a sistemi contenenti 
circa 10 elettroni. 

Per poter utilizzare la meccanica quantistica per i nostri scopi dovremo 
dunque ricorrere a nuovi strumenti matematici che ci permettano di tra- 
sformare l'equazione scritta e quindi introdurre diverse approssimazioni cer- 
cando in qualche modo di mantenere il controllo, durante i calcoli, di tutte 
le approssimazioni eseguite. La teoria che illustreremo all'interno di questo 
lavoro costituira proprio quell'insieme di approssimazioni e teoremi che ci 
permetteranno di giungere, partendo dall'equazione ( II. ip . ad un problema 
risolvibile. 

Come vedremo, un ruolo chiave nella teoria verra rivestito da equazioni 
di Schroedinger a particella singola che costituiscono, come gia detto, l'unico 
caso in cui siamo in grado di ottenere una soluzione. Un ruolo chiave sara 
svolto inoltre dalle grandezze integrate, in particolare la densita del sistema 
e la funzione di Green. Le equazioni che troveremo saranno di tipo auto- 
consistente. Cio significa che per studiarle un ruolo dominante sara svolto dai 
calcoli numerici, poiche solo grazie all'avvento dei moderni computer siamo in 
grado di risolvere tali equazioni che, altrimenti, rimarrebbero non risolvibili 
cosi come l'equazione iniziale. 

La teoria non verra presentata secondo l'evoluzione storica che ha avuto, 
piuttosto secondo quello che e lo schema logico che ho utilizzato per studiarli 
e comprenderli. 

Le due teorie piu importanti che andremo ad illustrare nei prossimi capi- 
toli saranno la teoria della funzione di Green, che indicheremo come MBPT 
(Many-Body Perturbation Theory), e la teoria del funzionale densita, che 
indicheremo come DFT (Density Functional Theory). 



1.2 L'approssimazione di Born-Oppenheimer 
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I teoremi di cui ci serviremo per lo sviluppo della teoria saranno: il teo- 
rema di Hohenberg-Kohn e il teorema di Kohn e Sham per la DFT, e il 
loro equivalente per l'estensione al caso dipendente dal tempo (TD-DFT); il 
teorema di Gell-Mann e Low e il teorema di Wick per la MBPT. 

Le prime due approssimazioni che introdurremo saranno l'approssima- 
zione di Born-Oppenheimer, che ci servira per disaccoppiare il moto degli 
elettroni da quello degli ioni e l'approssimazione di potenziale statico per il 
campo di Coulomb (approssimazione che abbiamo in realta gia introdotto 
scrivendo l'equazione di Schroedinger), poiche considerando i ritardi anche le 
equazioni della fisica classica diventano troppo complesse, immaginate quindi 
quelle della meccanica quantistica. Le altre approssimazioni principali di cui 
parleremo saranno la LDA (Local Density Approximation) per la DFT e la 
A-LDA (Adiabatic-LDA) per il caso dipendente dal tempo; l'approssimazione 
GW (chiariremo inseguito il significato) per la formulazione Many-Body. 

Le equazioni piu importanti che verranno nominate per le due teorie sa- 
ranno infine le equazioni di Kohn e Sham per la DFT ed il loro equivalente per 
il caso dipendente dal tempo (TDDFT); le equazioni di Hedin e l'equazione 
di Bethe-Salpeter per la MBPT 

1.2 L'approssimazione di Born-Oppenheimer 

L'approssimazione di Born-Oppenheimer consiste nel considerare disaccop- 
piati il moto degli elettroni e il moto dei nuclei che costituiscono il nostro 
sistema partendo dalla considerazione che essi presentano masse molto diffe- 
rent i e quindi evolvono su di una scala di tempo molto differente. Dal punto 
di vista matematico cio significa che possiamo idealmente separare la nostra 
equazione in due parti, la prima per il moto degli elettroni, nella quale le 
posizioni degli atomi vengono considerate dei parametri e la seconda che pos- 
siamo risolvere per trovare la posizione degli atomi una volta che abbiamo 
ottenuto la funzione d'onda degli elettroni. Per il momento ci concentreremo 
solamente sull'equazione del moto degli elettroni: 



La ragione che ci induce a concentrarci, per il momento, sulla risoluzione 
dell'equazione per gli elettroni nasce dal dal fatto che gran parte delle pro- 
priety del sistema dipenderanno proprio dalla soluzione di questa equazione, 
mentre per quanto riguarda gli ioni ci bastera semplicemente conoscere la loro 
posizione di equilibrio (e non la funzione d'onda) principalmente per deter- 
minare i parametri dell'equazione elettronica. Se volessimo calcolare anche 
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il moto degli ioni questo potrebbe ad esempio essere fatto in approssimazio- 
ne armonica per studiare le proprieta vibrazionali del sistema. All'interno di 
questa tesi comunque ci occuperemo solo dello studio delle eccitazioni elettro- 
niche e dunque gli ioni, una volta trovata la posizione di equilibrio, saranno 
considerati fermi. 

1.3 II ruolo dell'interazione, le grandezze in- 
tegrate e le funzioni d'onda a particella 
singola 

L'equazione resta comunque un problema formidabile. La maggiore difficolta 
che deve essere affrontato e dovuta al fatto che il nostro sistema presenta un 
numero di particelle molto elevato e che queste sono tra loro interagenti. La 
conseguenza matematica e dunque quella di avere un'equazione che presenta 
un numero molto elevato di variabili che devono essere studiate contempora- 
neamente. Entrambe le teorie che abbiamo citato in precedenza, la DFT e la 
teoria MBPT, sono state elaborate appunto per superare tale problema. 

II problema del numero elevato di variabili viene superato da entrambe 
le teorie grazie all'utilizzo delle "grandezze integrate" , cioe di funzioni in cui 
tutte le variabili, tranne poche, sono state integrate. Si tratta della densita 
per la DFT e della funzione di Green per la MBPT. L'idea e quella di non 
voler necessariamente calcolare la funzione d'onda dell'mtero sistema ma di 
calcolare quantita piu maneggevoli, che siano pero in grado di descrivere tutte 
le proprieta del sistema. 

Anche in questo modo pero tali teorie restano inutilizzabili poiche la pre- 
senza dell'interazione rendera impossibile anche il calcolo di queste due gran- 
dezze. Nell'affrontare questo problema le due teorie presentano un approccio 
completamente differente. 

La DFT riesce in modo elegante a spostare il ruolo dell'interazione al- 
l'interno di un potenziale esterno che agisca su un sistema di particelle non 
interagenti. A partire da questo, che sottolineiamo in principio non ha senso 
fisico, vengono poi ricavate le grandezze fisiche ed in particolare la densita e 
l'energia dello stato fondamentale. II pregio maggiore di questa teoria resta la 
semplicita delle equazioni che, a partire da essa, possono essere costruite. In 
particolare grazie a questo tipo di approccio si riesce a superare il problema 
del "Muro Esponenziale" di cui abbiamo parlato in precedenza [6]. II limite 
maggiore invece e costituito dal fatto che la forma di tale potenziale non e 
noto e che dunque la teoria necessiti di un input esterno per poter essere 
utilizzata. 



1.3 II ruolo dell'interazione, le grandezze integrate e le funzioni 
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La MBPT affronta invece il problema cercando di giungere alia soluzio- 
ne attraverso un approccio di tipo perturbativo a partire dalla soluzione del 
sistema non interagente. II punto di partenza di tale approccio e quello di 
considerare l'interazione come la grandezza perturbante del sistema. La teoria 
viene poi riscritta sotto forma di un set di equazioni che descrivono le grandez- 
ze fisiche caratteristiche del sistema. In questa formulazione si partire dalla 
soluzione del sistema non interagente e si calcolno via via le grandezze fisiche 
del sistema in modo iterativo. Cosi facendo ad ogni iterazione si aggiungono 
le correzioni alia soluzione del sistema non interagente dovute all'mterazione. 
Questo teoria offre il vantaggio, dal punto di vista dell'autore, di essere piu 
facilmente interpretabile grazie anche ai diagrammi di Feynman. Lo svan- 
taggio e in primo luogo quello di richiedere strumenti matematici molto piu 
raffinati e di presentare calcoli molto piu complessi della DFT. 

Concludiamo questo capitolo anticipando che cio che viene fatto in pra- 
tica e utilizzare equazioni che sono il risultato dell'incrocio delle due teorie, 
spesso andando al di la di cio che e stato formalmente giustificato. Cio di cui 
parleremo in questa tesi infatti si avvicina spesso a quelli che sono le fron- 
tiere della ricerca e a volte i modelli utilizzato trovano giustificazione solo a 
posteriori nel fatto che sono in grado prevedere in modo corretto i risultati 
sperimentali. 
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Capitolo 2 

La teoria del funzionale densita 



2.1 II teorema di Hohenberg-Kohn 
2.1.1 Dimostrazione del teorema 

In questo capitolo presenteremo quelle che sono le basi della teoria del fun- 
zionale densita. II punto di partenza di tale teoria e il teorema di Hohenberg- 
Kohn, il quale afferma che, dato un sistema fisico di particelle in cui l'in- 
terazione e fissata, esiste una relazione biunivoca tra il potenziale esterno 
applicato a tale sistema, la funzione d'onda dello stato fondamentale e la 
densita di particelle dello stato fondamentale 

V ext Po • 

Dimostriamo il teorema scegliendo come modello un sistema di particelle 
descritto dalla seguente Hamiltoniana 

H = f + V ext + w ; 

dove T rappresenta Foperatore energia cinetica, w e l'interazione e V ext il 
potenziale esterno. Immaginiamo ora che l'interazione sia fissata, mentre 
che il potenziale esterno possa essere libero, con la sola condizione pero che 
lo stato fondamentale del sistema risulti non degenere. Delle relazioni tra 
potenziale esterno, funzione d'onda dello stato fondamentale e densita dello 
stato fondamentale due sono, in principio, banali: la funzione d'onda puo 
essere infatti ricavata a partire dal potenziale esterno risolvendo l'equazione 
di Schroedingei0 e la densita dalla funzione d'onda 

Po = (#o|p|*o). (2.1) 

^ffcrmiamo qui che la relazione c banale poiche nclla dimostrazione del teorema siamo 
semplicemente interessati alia sua esistenza in principio. Resta invece non banale utilizzarc 
in pratica tale relazione 
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Restano dunque da dimostrare le due relazioni inverse. 

Partiamo dalla prima delle due implicazioni. Dimostreremo per assurdo 
che fissare per il nostro sistema uno stato fondamentale significa fissare il 
potenziale esterno a meno di una costante additiva, che pero non riveste 
alcuna importanza dal punto di vista fisico poiche l'energia di un sistema 
puo sempre essere stabilita a meno di una costante. Ipotizziamo dunque che 
esistano due differenti potenziali esterni che diano come soluzione lo stesso 
stato fondamentale 

(f + w + V ext )\V ) = E \y ) 
(f + w + V' ext )\%)=E' \%) , (2.2) 

con |#o> = TO- 
Se ora sottraiamo la prima equazione alia seconda otteniamo 

(V ext -V' ext )\^ ) = (E -E' )\^ ) . 

Applicando ad entrambe le equazioni l'operatore (^o| e assumendo che la 
funzione d'onda sia normalizzata ad uno otteniamo 

V ext (x) - V' ext (x) = (E - E' Q ) , (2.3) 

il che dimostra appunto, assumendo che le funzioni d'onda siano eventual- 
mente uguali a zero solo in regioni dello spazio di misura nulla, che i due 
potenziali possono differire al piu di una semplice costante additiva. 

Dimostriamo ora la seconda implicazione, vale a dire che fissata la den- 
sita risulta stabilita anche la funzione d'onda. A questo scopo dimostriamo 
che due funzioni d'onda differenti non possono generare la stessa densita, in 
altre parole che fissato j^o) 7^ l^o) risultera po 7^ Pq- P er 1 & dimostrazione 
utilizziamo il fatto che |^o) e lo stato fondamentale del sistema e dunque lo 
stato a minima energia: 

e = ro < ro#ro , 

dove H = T + w + V ext e l'Hamiltoniana associata al potenziale esterno che 
ha |\l/o) come ground state. Riscriviamo ora l'equazione come 

(%\H' - V' ext + V ext \%) = E' + J d 3 r p'(r) [V ext (r) - V' ext (r)] , (2.4) 

dove H' = T + w + V' ext e abbiamo utilizzato il fatto che l'operatore poten- 
ziale e un operatore moltiplicativo in spazio (r). Otteniamo cosi la seguente 
relazione 

E <E' + J d 3 r p'(r) [V ext (r) - V' ext (r)] . (2.5) 



2.1 II teorema di Hohenberg-Kohn 
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Ripetendo lo stesso ragionamento partendo dalla funzione d'onda \^' Q ) 
otteniamo l'equazione 

E' <E + J dh p(r) [V' ext (r) - V ext (r)] . 

Se infine ipotizziamo che p(r) = p'(r) e sommiamo le due equazioni otteniamo 

Eq + E' < E + E' , 

il che e ovviamente impossibile. 

Cio che abbiamo sin qui dimostrato e in realta che le seguenti mappe 

F : V ext -> (2.6) 

sono iniettive. Affinche esista invece una corrispondenza biunivoca tra funzio- 
ni d'onda, densita e potenziali abbiamo bisogno che esse siano anche suriet- 
tive. La questione potrebbe apparire superflua poiche se vogliamo utilizzare 
tale teoria per descrivere sistemi fisici reali (e se crediamo nella validita dell'e- 
quazione di Schroedinger) e ovvio che utilizzeremo funzioni d'onda che sono 
associate ad un qualche potenziale esterno. Alio stesso modo restera dunque 
ovvio che utilizzeremo, per la densita, funzioni che abbiano un significato fi- 
sico e che dunque debbano essere ottenibili a partire da una qualche funzione 
d'onda. In realta, poiche proprio la densita costituira, come vedremo, il pun- 
to di partenza di tale teoria gionevole chiedersi quale debba essere lo 
spazio delle funzioni da cui sia necessario partire nello svolgere i calcoli. Resta 
dunque non banale la suriettivita della seconda delle due mappe rispetto alio 
spazio delle funzioni. 

Lo studio di tale problema necessita la formulazione di ulteriori teoremi ed 
una discussione di questo fatto esula dagli scopi di questo lavoro. Altro punto 
di cui non ci occuperemo all'interno di questa tesi e quello dell'estensione 
di tale teorema a sistemi che presentano uno stato fondamentale degenere. 
Per chi fosse interessato alia questione e possibile trovare una discussione 
approfondita del problema nella referenza [2]. 

2.1.2 Lo scopo del teorema di Hohenberg Kohn e il 
ruolo della densita 

Come abbiamo gia accennato in precedenza la densita del sistema vuole essere 
la grandezza da cui partire per studiare un sistema a molti corpi nell'ambito 
della teoria del funzionale densita. Si tratta infatti di una grandezza che, a 
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differenza della funzione d'onda, e funzione di una sola variabile e che dunque 
puo essere piu facilmente controllabile. L'idea per poter utilizzar il teorema 
appena dimostrato e quella di formulare un principio di minimo partendo 
proprio da questa grandezza. Utilizzando infatti la mappa inversa della fun- 
zione G possiamo costruire, in linea teorica, una funzione d'onda a partire 
da una data densita. Variando la densita nello spazio delle funzioni ad una 
variabile ci muoveremo nello spazio delle funzioni d'onda individuando tut- 
te quelle che sono lo stato fondamentale di un qualche sistema descritto da 
un'Hamiltoniana del tipo 

H\p] = f + w + V ext \p] , 

dove abbiamo indicato che il potenziale esterno sara il funzionale della densita 
stabilito dall'inversione della funzione F ■ G. Se ora costruiamo la seguente 
grandezza, 

E = {9\p]\H\V\p]) , 

dove pero questa volta l'Hamiltoniana e quella che descrive il nostro sistema 
fisico e in cui dunque il potenziale esterno e fissato, otteniamo un funzionale 
il cui minimo corrispondera con il prendere la funzione d'onda dello stato 
fondamentale e dunque la funzione d'onda e la densita che, all'mterno dello 
schema di Hohenberg-Kohn sono legate all'Hamiltoniana di partenza e che 
hanno significato fisico. 

Abbiamo dunque raggiunto uno dei nostri scopi, cioe quello di costruire 
una teoria che possa partire semplicemente dalla densita del sistema piut- 
tosto che dalla sua funzione d'onda e quindi dalla soluzione dell'equazione 
di Schroedinger a molti corpi. Per trovare infatti lo stato fondamentale del 
nostro sistema non dovremo far altro che minimizzare, rispetto alia densita 
il funzionale appena costruito. II problema di questo approccio consiste pero 
nel fatto che il funzionale passa di fatto ancora tramite la costruzione della 
funzione d'onda molti corpi e che il teorema di Hohenberg-Kohn non fornisce 
alcun metodo pratico per costruire tale funzione. II problema della costruzio- 
ne della funzione d'onda resta in ultima analisi legato a quello del fatto che 
le particelle sono interagenti. II modello che riesce a superare questo ostacolo 
verra spiegato nel prossimo paragrafo ed e stato formulato un anno dopo il 
teorema appena dimostrato. Come abbiamo detto nell'introduzione, sebbe- 
ne la grandezza fondamentale per questa teoria resti la densita del sistema, 
all'interno di questo modello assumeranno un ruolo fondamentale le funzio- 
ni d'onda di particella singola che sono le uniche funzioni d'onda che siamo 
in grado di calcolare a partire dall'equazione di Schroedinger. Sottolineiamo 
pero che tali funzioni d'onda non avranno di fatto nessun significato fisico ma 
saranno semplicemente un artifizio utilizzato per raggiungere i nostri scopi. 



2.2 Lo schema di Kohn e Sham 
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2.2 Lo schema di Kohn e Sham 

2.2.1 Perche scegliere di lavorare con un sistema non 
inter agent e 

Lo schema di Kohn e Sham nasce dall'idea di utilizzare come supporto per i 
nostri scopi un sistema di particelle non interagenti, in modo che, nei calcoli 
che dovremo fare, ci troveremo a lavorare con equazioni di Schroedinger a 
particella singola che saremo quindi in grado di risolvere. La scelta di lavora- 
re con le funzioni d'onda di un sistema non interagente, anziche direttamente 
con la densita del sistema, sembrerebbe una complicazione della teoria sin 
qui presentata poiche in questo modo ci troveremo a lavorare con grandez- 
ze a molte variabili e dunque dovremo nuovamente affrontare il problea del 
numero elevato di particelle. D'altra parte questa scelta offre il vantaggio di 
essere un punto di partenza migliore per introdurre delle approssimazioni che 
consentano di utilizzare in pratica la teoria presentata. In particolare il suc- 
cesso di questa scelta e in gran parte dovuto al successo dell'approssimazione 
LDA di cui parleremo nel capitolo 4 di questa tesi. 

Anche in questo tipo di approccio la densita resta comunque la grandezza 
fondamentale dello schema e il punto di partenza della teoria. 



2.2.2 Le equazioni di Kohn e Sham 

L'idea e dunque quella di utilizzare il teorema di Hohenberg-Kohn per un 
sistema ausiliario in cui l'interazione e nulla. In questo modo otterremo a 
partire da una densita di prova una funzione d'onda di particelle indipendenti 
^ s [p] soluzione di un'Hamiltoniana con potenziale esterno V^ ea;t [p] entrambi 
costruiti a partire dalle mappe F s e G s per un sistema non interagente. 

(f + v s [p])^s[p]=E^ s [p] . (2.7) 

L'apice "s" indica appunto il fatto che facciamo riferimento ad un sistema in 
cui w = 0. 

Vogliamo ora provare a riscrivere il funzionale della densita da cui rica- 
viamo Fenergia del sistema nel suo minimo partendo dalla funzione d'onda 
^sfp]. L'idea per tale funzionale e quella di fare il valor medio dell'Hamil- 
toniana non interagente rispetto alia nuova funzione d'onda e di inserire poi 
dei termini correttivi per tener conto dell'interazione e della differenza tra la 
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funzione d'onda reale e quella non interagente. 

E.[p]=(*M\T + v ext \^ s [p))+ [ [ l r) _ p( jl ) rf3r dV + ^M 



(2i 



T s [p} + U ext [p] + E H [p] + E xc [p] 



Cerchiamo di analizzare il funzionale scelto. II primo termine e costituito 
dall'energia cinetica delle particelle non interagenti (T s [p]) poiche la speranza 
e quella che non differisca di molto dall'energia cinetica del sistema reale. II 
secondo termine e costituito dall'energia dovuta al potenziale esterno (U ext ) 
e poiche questa dipende solo dalla densita sara identica a quella del sistema 
interagente. II terzo termine e costituito dall'energia di Hartree (Eh[p\) cioe 
l'energia dovuto a quello che e i campo medio degli elettroni. E stato inserito 
poiche e Tunica parte dell'interazione che sappiamo scrivere in modo esatto 
sotto forma di funzionale della densita e di cui quindi sappiamo calcolare 
l'energia; anch'esso dipende unicamente dalla densita e corrisponde all'energia 
di campo medio del sistema interagente. L'ultimo termine infine viene definito 
come la differenza tra il funzionale costruito per il sistema interagente e i 
termini sin qui considerati del funzionale appena costruito. Esso in particolare 
tiene conto della differenza tra l'energia cinetica dei due funzionali e di quella 
dovuta tra l'energia di Hartree e l'energia dovuta all'interazione. 

E xc [p] = (T [p] - T s [p] ) + (E w [p] - E H [p] ) , 

dove con E w abbiamo indicato il valor medio dell'interazione w rispetto alia 
funzione d'onda del sistema interagente. Scritta in particolare l'espressione 
dell'interazione come 



w 




d 3 rd 3 r' 1 — ft(r)ft(r')4>(r')4>(r) , (2.9) 
possiamo ora defmire la "funzione di correlazione a coppie" 

g{l y M= (^l-^l , (2 , ) 

p(r)p(r') p(r) 
e scrivere in modo esplicito il valor medio {w (r, r')) 

il che ci da un'idea della differenza di questo rispetto all'energia di Hartree. 

Poiche il funzionale cosi costruito e identico a quello del sistema intera- 
gente se ora cerchiamo il minimo rispetto alia densita otterremo, come in 
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precedenza, la densita e l'energia di stato fondamentale del nostro sistema. 
Cio che invece non avra significato fisico sara la funzione d'onda. 

Possiamo ora sfruttare il fatto che stiamo studiando un sistema di parti- 
celle non interagenti e scrivere la funzione d'onda come il prodotto di funzioni 
d'onda di particella singola. In questo modo in particolare la densita del si- 
stema non sara altro che la somma dei moduli quadri delle funzioni d'onda. 
Pensando ora al nostro funzionale come a un funzionale composto scegliamo 
di minimizzarlo rispetto alle funzioni d'onda anziche rispetto alia densita. 
Poiche pero le funzioni d'onda rispetto a cui vogliamo minimizzare il nostro 
funzionale non possono variare liberamente nello spazio delle funzioni d'onda 
(che e il prodotto tensoriale degli spazi di Hilbert di particella singola) ma 
devono essere ortonormali dovremo cercare un minimo vincolato. Utilizzan- 
do la teoria dei moltiplicatori di Lagrange l'equazione che dovremo risolvere 
sara: 

(2.12) 

II fatto che le funzioni d'onda sono ortogonali ci permette ora di pensare all'e- 
quazione appena scritta come ad un set di n equazioni differenti. Scriviamo 
dunque in modo esplicito tali equazioni, poiche esse costituiscono in effetti 
il set di equazioni a particella singola che siamo in grado di risolvere e che 
cercavamo. A questo scopo utilizziamo la seguente relazione tra la derivata 
funzionale rispetto alle funzioni d'onda e quella rispetto alia densita 

5 dp 5 5 

ipi — , (2.13) 



Sip* dip* 5p 5p 

in modo da poter fmalmente scrivere l'equazione a particella singola nel 
seguente modo: 

w Ts + ^% ^ Eh + uext + Exc ^ = XlA 

(f + v H [p] + V ext + v xc [p\) A = ■ (2.14) 

Siamo dunque giunti alia fine delle nostre fatiche, almeno per quanto con- 
cerne la formulazione della teoria che dovremo poi utilizzare. Cerchiamo di 
analizzare i termini di questa equazione e capire come possiamo pensare di 
utilizzarla ai nostri scopi. Come volevamo e un'equazione a particella singola, 
abbiamo un termine per l'energia cinetica, un termine per il campo medio, 
il potenziale esterno e un potenziale che tiene conto del fatto che il risultato 
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della nostra equazione deve essere un set di funzioni d'onda che generi una 
densita identica a quella di un set di particelle reali poste all'interno dello 
stesso potenziale esterno presente nella nostra equazione. Otteniamo dun- 
que che, se sommiamo in un unico termine il potenziale esterno, il potenziale 
di Hartree e il termine di scambio e correlazione otteniamo esattamente il 
potenziale v s introdotto nell'equazione (12. 7p . 

2.3 Alcune considerazioni 
2.3.1 Lo spin 

Finora abbiamo parlato di DFT prendendo in considerazione la densita e 
abbiamo mostrato come Fenergia dello stato fondamentale del sistema possa 
essere espressa come funzionale di questa grandezza. Tutte le dimostrazioni 
sin qui eseguite restano di fatto valide se il potenziale esterno puo essere 
espresso nella forma scelta per la dimostrazione del teorema, ovvero nel caso 
in cui non siano presenti campi magnetici. Se volessimo includere gli effetti di 
un campo magnetico, perlomeno l'accopiamento del campo con il momento 
magnetico dovuto alio spin 



+ / d 3 r p'(r) [V ext (r) - V' ext (r)] - m[p'](r) ■ [B(r) - B'(r)] , (2.16) 



il che porterebbe pero ad avere un potenziale che agisce in modo non locale 
sulla funzione d'onda nello schema di KS [2]. 

Quello che dunque viene solitamente fatto b sviluppare la teoria scrivendo 
un funzionale della densita e della magnetizzazione (o della matrice densita 
di spin, il che e equivalente) in modo da poter lavorare con dei campi locali. 
Non riportiamo qui la dimostrazione di questa estensione dei teoremi visti; 
nel corso del nostro lavoro faremo comunque uso di questo funzionale 




d 3 r [V ext {r)p{r) - B(r) • m(r)] 



(2.15) 



dovremmo riscirivere il passaggio ( 12. 4p nel seguente modo 



{V' \H' - V' ext + V ext \%) = E' 



E = E[p a p] 



(2.17) 



poiche tale teoria viene utilizzata anche per lo studio di sistemi che presentano 
uno stato fondamentale spin-polarizzato. 



2.3 Alcune considerazioni 



19 



Poiche in realta all'interno di questo lavoro ci limiteremo a studiare sistemi 
collineari e perturbazioni con campi magnetici diretti secondo l'asse z sara 
sufficiente studiare il seguente funzionale 



in cui Fenergia e cioe funzionale delle sole componenti diagonali della matrice 
densita di spin. Vedremo piu in dettaglio quale sia il significato di tale scelta 
nel capitolo 4. 

Se infine volessimo inserire anche l'accoppiamento del campo magnetico 
con la corrente elettronica dovremmo sviluppare ulteriormente la teoria. Cio 
esula pero dagli scopi di questo lavoro. 

2.3.2 La teoria costruita 

Abbiamo costruito una teoria che possa essere utilizzata per studiare un si- 
stema di particelle interagenti. Come volevamo siamo riusciti a semplificare 
l'equazione di Schroedinger in modo da dover studiare solo equazioni a parti- 
cella singola, partendo da una grandezza integrata (la densita). Le quantita 
che siamo in grado di calcolare sono l'energia dello stato fondamentale del 
sistema e la sua densita. Non possiamo invece dire niente, almeno secondo 
giustificazione formale, della funzione d'onda. 

La cosa piu notevole del nostro risultato credo consista nel fatto che nei 
passaggi sin qui svolti non abbiamo introdotto nessuna approssimazione. Le 
uniche approssimazioni fin qui utilizzate sono quella di Born-Oppenheimer (di 
cui abbiamo gia parlato) e il fatto che stiamo utilizzando un'interazione di ti- 
po statico. Ovviamente sarebbe tutto troppo bello se potessimo accontentarci 
di queste due approssimazioni. La nostra teoria presenta infatti un elemento 
oscuro che e appunto il potenziale di scambio e correlazione (o se preferi- 
te l'energia di scambio e correlazione). Per una discussione sulla questione 
rimandiamo nuovamente al capitolo 4 del presente lavoro. 

Una considerazione a parte merita la questione dello studio dello spettro 
di eccitazione del sistema. Fino ad ora infatti ci siamo limitati a parlare 
di energia di stato fondamentale e abbiamo sviluppato un metodo pratico, 
lo schema di Kohn e Sham, per poter calcolare tale grandezza. Ritornan- 
do al teorema di Hohenberg-Kohn pero possiamo osservare come, poiche il 
potenziale esterno puo essere scritto come funzionale della densita a partire 
da questo, tramite la soluzione dell'equazione di Schroedinger potremmo, in 
principio, riuscire a calcolare tutti gli stati eccitati del nostro sistema. In 
pratica avremmo bisogno di sviluppare un metodo, equivalente alio schema 
di Kohn e Sham, per il computo delle energie di eccitazione del sistema. Nel 



E = E[p t , Pi ] 
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capitolo sullo studio dello spettro di eccitazione vedremo come tale metodo 
pratico e in effetti costituito dalla Time-Dependent DFT (TDDFT) per quel 
che riguarda lo studio delle energie di eccitazione per stati che presentano 
lo stesso numero di particelle dello stato fondamentale. Non ci occuperemo 
invece dello studio delle energie di eccitazione di stati con un numero di par- 
ticelle differenti dal quello del stato fondamentale. Per chi fosse interessato 
al problema rimandiamo come sempre alia referenza [2] e all'articolo [7]. 



Capitolo 3 

La teoria della funzione di 
Green 

3.1 Introduzione 

In questo capitolo cercheremo di introdurvi a quella die e la "Many Body 
Perturbation Theory" o se preferite "Teoria della funzione di Green". Per la 
comprensione del capitolo sara richiesta, oltre che la conoscenza della mecca- 
nica quantistica secondo la formulazione di Schroedinger anche il formalismo 
matematico della seconda quant izzazione e dunque la definizione di operato- 
re di campo. Cercheremo di delineare tutti quei passaggi che portano alia 
costruzione delle formule matematiche che verranno poi utilizzate all'interno 
di questo lavoro e, per quanto possibile di riportare anche una dimostrazione 
dei teoremi che via via citeremo. Poiche d'altro canto questa teoria e di fatto 
molto piu elaborata della DFT non ci sara possibile presentare in modo sem- 
plice e lineare tutti i suoi risultati. Saremo quindi spesso costretti a rimandare 
il lettore ad altri libri specialistici o a pubblicazioni sull'argomento. In parti- 
colare all'interno di questa tesi faremo spesso riferimento al testo di Fetter e 
Walecka [I], "Quantum Theory of Many-Particle Systems" che costituisce il 
libro su cui ho studiato tale teoria. 

3.2 La funzione di Green 
3.2.1 Definizione 

La MBPT cerca di descrivere il nostro sistema di particelle interagenti a par- 
tire da una grandezza integrata, la funzione di Green, cosi come la DFT 
partiva dalla densita del sistema. Se pero e certo chiaro a tutti cosa sia la 
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densita e sicuramente piu difficile spiegare cosa sia la funzione di Green. A 
dire il vero esistono perlomeno due differenti definizioni di questa grandezza 
e, non avendo ben chiaro quale sia il suo significato fisico, risulta ancora piu 
complicato capire perche le differenti definizioni possano portare alia stessa 
grandezza. La scelta che faremo all'interno di questo lavoro sara quella di 
presentare innanzi tutto la funzione di Green secondo la definizione che e alia 
base della MBPT. A partire da questa cercheremo di capire quale possa essere 
il suo significato fisico e quale sia il motivo per cui sia stata scelta tale gran- 
dezza per descrivere il nostro sistema. II riferimento all'altro modo in cui e 
possibile definire tale funzione verra fatto unicamente in un secondo momento 
per aiutarci nella comprensione del significato fisico di tale grandezza. 
La funzione di Green e dunque definita nel seguente modo: 



z'G(l,2) 



(*o|T 




l*o) 




(*o|*o) 




(*o|T 


4(1)^(2) 


l*o) 





(3.1; 



(^ol^o) 



Nella seconda linea abbiamo introdotto la convenzione 1 = (xi,ii,0i) di cui 
faremo uso a volte all'interno di questa tesi. Questa grandezza e in particolare 
la funzione di Green T-ordinata a particella singola. E dunque definita come 
il valor medio, rispetto alio stato fondamentale j^o) di un sistema descritto 
da un'Hamiltoniana indipendente dal tempo. Possiamo immaginare sia la 
stessa Hamiltoniana che abbiamo utilizzato nel presentare la DFT e dunque 
l'Hamiltoniana del nostro problema many-body introdotto nel primo capitolo 
di questa tesi 

H = f + w + V ext , (3.2) 

dove ancora T rappresenta l'operatore energia cinetica, w l'interazione e V ext 
il potenziale esterno del sistema, nel nostro caso quello degli ioni. II valor 
medio e fatto sugli operatori di campo in descrizione di Heisenberg (vedi pa- 
ragrafo successivo). I due operatori sono sotto l'operatore di T-ordinamento, 
o ordinamento temporale, che e definito nel seguente modcQ 



T 



^(l)ft(2) =0(t 1 -t 3 )^(l)^(2)-0(t 2 -t 1 )^(2)^(l) . (3.3) 



In altre parole i due operatori vengono T-ordinati facendo agire prima 
quello il cui argomento temporale e precedente. II segno meno tra i 2 operatori 



Piu in generale il T-ordinamento di N operatori di campo c definito come 
T^}{t^)<ij}{tx)^}{t2)\ = (—I) p tfj(ti)ilj(t2)ip(t3) con P ugualc al numcro di permutazioni 
eseguite e con ti > t 2 > t 3 
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e dovuto al fatto che stiamo studiando un sistema di fermioni e die quindi gli 
operatori di campo sono definiti secondo regole di anti-commutazione^|. 

3.2.2 Rappresentazione di Schroedinger, Heisenberg e 
interazione 

Nella formulazione piu semplice della meccanica quantistica, che chiamiamo 
rappresentazione di Schroedinger, l'evoluzione temporale del sistema e de- 
scritta dall'evoluzione delle funzioni d'onda o stati del sistema. Come noto 
infatti sono le funzioni d'onda a descrivere tutte le proprieta del sistema. 
Quando ci si trova pero a dover studiare sistemi con un numero molto elevato 
di particelle, come abbiamo gia detto, le funzioni d'onda diventano oggetti 
difficili da controllare e si deve dunque ricorrere ad altro. II formalismo ma- 
tematico che viene utilizzato per questo scopo e appunto quello della seconda 
quantizzazione. In esso parte centrale viene svolta dagli operatori e in par- 
ticolare dagli operatori di campo, l'idea e dunque quella di spostare anche 
l'evoluzione temporale del sistema su questi oggetti. Vediamo come. 

Nella formulazione di Schroedinger le funzioni d'onda sono rappresentate 
dalle soluzioni dell'equazione 

in^\* s {t)) = H\* s (t)) , (3.4) 

dove il pedice S indica che e la rappresentazione di Schroedinger. L'operatore 
evoluzione temporale in questa descrizione risulta: 



(3.5) 



Se ora vogliamo una rappresentazione in cui l'evoluzione temporale venga 
tolta dagli stati (che chiamiamo rappresentazione di Heisenberg) possiamo 
definire 

\* H ) = e il ?\* H (t)) , (3.6) 
da cui si vede facilmente che \^h) e indipendente dal tempo: 

ih^\* H ) = ■ (3.7) 

Se ora vogliamo imporre che le due formulazioni siano equivalenti dobbiamo 
richiedere che i valori medi degli operatori di Heisenberg rispetto a tali stati 



La teoria della funzione di Green e in realta piu gcncrale e puo essere utilizzata anche 
per lo studio di bosoni. All'intcrno di questa tesi ci limitercmo pero a studiare la funzione 
di Green fcrmionica 
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siano uguali ai valori medi degli operatori di Schroedinger rispetto agli stati 
di Schroedinger. 

(* H \dH\*H) = {*s\d\*s) , (3-8) 
da cui utilizzando l'inversa della formula (13. 6p si vede facilmente che 

H = e l ^Oe~ 1 ^ , (3.9) 



dove gli operatori senza pedice sono quelli in descrizione di Schroedinger. 

La maggior difficolta del problema many-body e costituita dal fatto che 
le particelle siano tra loro interagenti. Un punto di partenza migliore per 
superare tale difficolta sarebbe quindi costituito da una descrizione che sepa- 
ra l'evoluzione temporale dovuta all'interazione da quella dovuta alia parte 
restante dell'Hamiltoniana. Vediamo come. 

Immaginiamo la nostra Hamiltoniana suddivisa in due parti, entrambe 
indipendenti dal tempo, 

# = #o + #i , (3.10) 

e definiamo lo stato come 

\9i(t))=e iS F\* s (t)) , (3.11) 

il pedice / sta per "interazione" . Vediamo ora qual'e la dipendenza dal tempo 
egli stati cosi costruiti 



iTi^jit)) = - H e l ^\q s (t)) + e l ^[H + Hi}\V s (t)) 

(3.12) 

. Hr\ t . H n t 

= e^H ie -^\y s (t)) , 



dove abbiamo usato l'inversa della (13. lip e il fatto che e % ~^~ commuta con 
Hq. Lo stato e dunque ancora dipendente dal tempo ma si vede che la sua 
dipendenza sara dovuta alia seconda parte dell'Hamiltoniana, moltiplicata 
per i due esponenziali. Richiedendo come in precedenza che anche per questa 
formulazione risulti 

(tf/ld/i*/) = (v s \d\y s ) , (3.i3) 

otteniamo 

. Hnt ~ . Hat 

7 = e l -tOe- l -t . (3.14) 

Dunque l'evoluzione dello stato e dovuta semplicemente all'operatore H\ in 
descrizione di interazione mentre l'operatore evolve secondo il solo operatore 
Hq. Se ora scegliamo Hq = T + V ext e Hi — w abbiamo ottenuto il nostro 
scopo: l'evoluzione libera e stata spostata sugli operatori mentre l'evoluzione 
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dovuta all'interazione resta sugli stati e in particolare puo essere racchiusa 
nell'operatore U(t, t ) che diventa il termine rilevante della teoria. 

Vediamo dunque come esprimere proprio tale operatore studiando Fevo- 
luzione di uno stato 



!*/(*)> 



e^\^ s (t)) = 

.Hgt , H(t-t ) 

e s e « e " 



.Hgt 

e l r e 



fi |*s(*o)> 



(3.15) 



*j(*o)> 



E possibile infine mostrare come [T] l'operatore di evoluzione temporale cosi 
scritto puo essere ricondotto ad un'espressione del tutto simile a quella del- 
l'operatore che descrive Fevoluzione degli stati di Schroedinger, con la sola 
differenza che in questo caso l'Hamiltoniana che descrive Fevoluzione e dipen- 
dente dal tempo, poiche l'operatore H\ e in descrizione di interazione e che 
quindi "tutto e piu complicato" . II risultato che viene dimostrato in Ref. [I] e 
il seguente: 



U(t,t c 



n=0 
f 



dt n f[H 1 (t 1 )]...H 1 {t n )] 



to 



(3.16) 



L 'ultima riga e un'espressione formale. In particolare tale espressione descrive 
Fevoluzione di un qualsiasi stato in presenza di un'Hamiltoniana dipendente 
dal tempo, come in effetti e lo stato in descrizione di interazione essendo 
soluzione della f)3.12p . 



3.2.3 Perche la funzione di Green 

La scelta di utilizzare la funzione di Green come grandezza integrata da cui 
partire a studiare il nostro sistema puo essere motivata da diverse ragioni; uno 
dei motivi principali resta pero il fatto che a partire da essa e possibile ricavare 
il valore di aspettazione rispetto al ground state di un qualsiasi operatore a 
particella singola e che b inoltre possibile calcolare Fenergia totale del sistema. 

La prima delle due affermazioni e facilmente verificabile. Consideriamo 
un generico operatore A a particella singola (per comodita di notazione sce- 
gliamo un operatore indipendente dal tempo, per operatori dipendenti dal 
tempo la dimostrazione sarebbe identica) e scriviamone la sua rappresenta- 
zione in spazio (x, a) prima in rappresentazione di Schroedinger, quindi in 
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rappresentazione di Heisenberg: 



A=J2 y"d 3 xid 3 x 2 A a , |8 (x 1 ,x 2 )^ 1 (x 1 )Vi CT2 (x 2 ) 



AH(t) = ^2 / d 3 xid 3 x 2 A CTl)CT2 (xi,x 2 )V)t o _ i (xi,t)^ iCT2 (x 2 ,t) 



(3.17) 



Se ora facciamo il valor medio rispetto alio stato fondamentale in descrizione 
di Heisenberg della seconda riga e inseriamo una variabile t 2 — > t + otteniamo 
(e scegliendo t\ = t) 



(V \A H \W ) 



d 3 xid 3 x 2 lim t4 CT1i(72 (xi, x 2 )G' (7liCT2 (-x.it, x 2 t 2 ) 



(3.18) 

II simbolo t + indica un tempo che e un infinitesimo successivo all'istante t 
ed e stato utilizzato per garantire che la funzione di Green sia ben definita. 
Se infine vogliamo studiare un operatore diagonale in spazio x dovremo ba- 
nalmente inserire un limite x 2 — > x in modo analogo a quanto abbiamo fatto 
per l'istante T. II computo dell'energia totale del sistema risulta invece piu 
complicata poiche per ottenere tale grandezza dobbiamo calcolare anche il 
valor medio di un operatore a due particelle quale e l'interazione che, in linea 
di principio, richiederebbe l'utilizzo di una funzione di Green a due particelle. 
L'ostacolo puo essere pero aggirato utilizzando l'equazione di Schroedinger; 
riportiamo il risultato che e possibile trovare sulla referenza pQ: 



E 



d 3 x lim lim 

t 2 ->i+ x 2 ->x 



d_ 
di 



h 



o (x 



G a>a (x.t, x 2 t 2 ] 



(3.19) 



dove ho = T + V ext Da questa Fenergia totale del sistema nel suo stato 
fondamentale puo essere facilmente ottenuta come la somma dei valori medi 
degli operatori che compongono la nostra Hamiltoniana. 



3.3 II teorema di Gell-Mann e Low 

3.3.1 Introduzione: l'accensione adiabatica 

Abbiamo gia anticipato che delle tre rappresentazioni di cui abbiamo parlato 
in precedenza quella che utilizzeremo sara la rappresentazione d'interazione. 
L'idea e quella di utilizzare infatti tale rappresentazione per legare in qual- 
che modo la funzione d'onda soluzione dell' Hamiltoniana completa a quella 
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soluzione dell'Hamiltoniana non interagente facendo evolvere quest 'ultima in 
descrizione di interazione. 

Partiamo dalla seguente Hamiltoniana: 

H(t) = H + e~ £ ^Hi , (3.20) 

dove nel nostro caso Hi e l'interazione e H q = T + V ext e l'Hamiltoniana del 
sistema non interagente, e e una costante positiva arbitraria. Notiamo che 
per t ^> e' 1 risulta H(t) ~ H e che per t = invece H(0) = H + Hi. 
Utilizziamo ora la descrizione di interazione^) e scriviamo Fevoluzione di uno 
stato dal tempo to & 1 tempo t: 

!*/(*)> =&(Mo)|*j(*o)> , (3-21) 

dove per ora lo stato I 1 !// (to)) non e speciflcato. Utilizziamo quindi nuovamen- 
te la deflnizione (13.111) e scegliamo come stato di partenza per la descrizione 
di Schroedinger lo stato fondamentale del sistema non interagente. Infine 
scegliamo come istante iniziale un tempo to ^ e 1 m modo che a tale istante 
l'evoluzione temporale dello stato fondamentale del sistema sia (quasi) quella 
libera e otteniamo 

|*/(*o)> = e^ |^ s (to)) ~ |$ ) , (3.22) 

cioe l'autostato a energia minima del sistema non interagente (potremmo 
partire da un autostato qualsiasi, o anche da una qualsiasi combinazione 
lineare di tali autostati, ma poiche siamo interessati ala ricerca dello stato 
fondamentale del sistema...). Se ora facciamo il limite per to —> — oo le 
relazioni scritte diventano esatte ed otteniamo: 

|tf/(t)> =f) £ (t,-oo)|$o> , (3.23) 

siamo cioe in grado di scrivere uno stato in descrizione di interazione in un 
istante t arbitrario. Cio che ci chiediamo ora e cosa accade se facciamo il 
limite per e — > 0, poiche siamo interessati a far sparire la dipendenza da e 
(non ha senso che i nostri risultati dipendano da una costante arbitraria) e 
poiche in questo caso l'Hamiltoniana si riduce appunto a quella che vogliamo 
studiare9. 

La risposta a tale domanda e tutt'altro che banale e costituisce il risultato 
del teorema di Gell-Mann e Low. 

3 E possibile dimostrare che la p,16[) e valida anche nel caso in cui l'Hamiltoniana 
dipenda dal tempo 

4 Parte essenziale della teoria e che i due processi di limite vengano eseguiti ncU'ordinc 
che abbiamo indicato, e cioe prima to — > +oo e poi e — > 0. L'inversione dei due limiti 
porterebbe ad una formulazionc del tutto inutile. 
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3.3.2 Enunciato 

L'enunciato del teorema di Gell-Mann e Low e dunque: 

Se la seguente quantita esiste a tutti gli ordini in teoria perturbativa 

lim ft(°.-°°)l*-> , (3.24) 
««{*o|f/,(0,-oo)|$„) 

e un autostato dell'Hamiltoniana H — H + Hi 

Per la dimostrazione di tale teorema rimandiamo alia bibliografia ppj. Se- 
gnaliamo pero come in effetti il teorema non parla di stato fondamentale 
dell'Hamiltoniana ma di semplice autostato. Gia nel paragrafo precedente 
abbiamo accennato al fatto che in effetti tutto cio che e stato fatto flnora vale 
per un qualsiasi autostato del sistema. II problema diventa dunque come indi- 
viduare lo stato fondamentale. In precedenza abbiamo suggerito di scegliere 
come stato di partenza lo stato fondamentale del sistema non interagente (<3>o 
appunto). II fatto che lo stato fondamentale non interagente porti alio stato 
fondamentale interagente non e in realta garantito in alcun modo dal teore- 
ma; questo punto rappresenta anzi una delle maggiori difficolta della MBPT 
ed e dovuta alia scelta di considerare l'interazione come una perturbazione 
del sistema. 



3.4 Interpretazione fisica della funzione di Green 
3.4.1 Analizziamo l'espressione 

Cerchiamo ora di dare un'interpretazione fisica alia funzione di Green ana- 
lizzandone l'espressione in descrizione di interazione. Partiamo dunque dal 
considerare uno stato \^j(tf)) a cui aggiungiamo una particella nel punto 
(xT) in modo da ottenere lo stato ipl^'t^l^iit')). Facciamo quindi evol- 
vere lo stato dal tempo t' al tempo t utilizzando l'operatore di evoluzione 
temporale in questa descrizione (13. 16[) U(t, t')i(}j(x!t')\^fi(t')). 

Ora per un istante qualsiasi t > t' consideriamo qual'e la sovrapposizione 
dello stato cosi ottenuto rispetto alio stato che si ottiene creando una parti- 
cella nel punto (xt), ip\(?d)\^> i(t)) . II risultato sara appunto la funzione di 
Green per t > t' . Per t < t' e possibile fare un ragionamento analogo a quello 
sin qui svolto, considerando la creazione di una buca all'istante t fatta poi 
evolvere fino all'istante t'Jfl 

5 Questa situazione puo anche essere interpretata come una particella creata al tempo 
t' che evolve indietro nel tempo. L'interpretazione c dovuta a Feynman in un suo famoso 
articolo 
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La funzione di Green descrive quindi il propagarsi di una particella (o di 
una buca) creata all'interno del gas elettronico e che quindi interagisce con 
esso. Piu in particolare essa ci da la probability di trovare una particella (o 
una buca), creata nel punto (x'i'), nel punto (x£). 

3.4.2 L'equazione del moto e confronto con la defini- 
zione classica di funzione di Green 

Per scrivere l'equazione del moto della funzione di Green non dobbiamo far 
altro che applicare l'operatore di derivazione rispetto al tempo alia definizione 
(13. ip . svolgere quindi in modo esplicito la derivazione dei vari termini (la 
funzione 6{t) e l'operatore di campo), utilizzare l'equazione del moto per gli 
operatori in descrizione di Heisenberg 

^P-^-iniHJUt)] (3.25) 

e infine calcolare il commutatore dell'Hamiltoniana con l'operatore di campo. 
Non riportiamo qui la dimostrazione che, come sempre, puo essere trovata 
nella Ref. pQ e scriviamo direttamente il risultato: 



ot 



G(l,2) = M(l,2) -ih w(l,3)G(l,3;2,3 + ) . (3.26) 



La notazione qui introdotta e l'estensione di quella gia utilizzata in preceden- 
za: 3 + = (x 3 ,cr 3 ,t^). Se eliminiamo l'ultimo termine di sinistra dell'equazione 
ritroviamo quella che e la definizione di funzione di Green utilizzata nell'e- 
lettromagnetismo classico per la soluzione delle equazioni delle onde e che 
puo essere facilmente interpretata come il campo generato da una sorgen- 
te di tipo deltiforme (nello spazio e nel tempo). II termine che compare in 
piu non rappresenta nient'altro che una "correzione" a tale definizione. Tale 
correzione e dovuta al fatto che il campo fermionico (a differenza di quel- 
lo elettromagnetico classico) "interagisce con se stesso". Alia luce di questa 
analogia possiamo capire l'interpretazione che in precedenza abbiamo fatto 
della funzione di Green. II creare una particella nel punto (x't') corrisponde 
al porre una sorgente puntiforme nello spazio e nel tempo che poi si propaga 
all'interno del gas elettronico. 

Ora l'analogia si arresta qui poiche nei due casi le funzioni di Green 
vengono utilizzate in modo differente. 

Per quanto riguarda l'elettromagnetismo classico viene utilizzato il princi- 
pio di sovrapposizione per ottenere il campo in presenza di una certa densita 
di carica attraverso l'integrale di convoluzione della funzione di Green con 
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tale densita; oltre che alia densita sorgente la soluzione dell'equazione viene 
poi definita imponendo le condizioni al contorno nello spazio e nel tempo. 

Nel caso fermionico non ha piu senso parlare di sorgente del campo (o per- 
lomeno non e in questo modo che si procede). Cio a cui siamo interessati e la 
funzione di Green stessa che descrive il modo in cui il campo si propaga men- 
tre interagisce con se stesso. Ancora una volta la difficolta principale nasce 
dalla presenza dell'interazione e quindi del termine in piu che in precedenza 
non compariva. Vediamo infatti come la funzione di Green a una particella, a 
causa della presenza dell'interazione sia legata a quella a due particelle; alio 
steso modo per una funzione di Green a n particelle si ricava un'equazione del 
moto che conterra una funzione di Green a n + 1 particelle. Per superare il 
problema dovremo dunque cercare di esprimere in modo diverso l'interazione 
che come sempre e la causa di tutte le nostre difficolta. Vedremo nei paragrafi 
successivi come cio sara fatto. 



3.5 La rappresentazione di Lehmann 

Prima di procedere con la presentazione della teoria many-body ci voglia- 
mo soffermare sull'esporre quella che e la rappresentazione di Lehmann della 
funzione di Green. In questa sezione le funzioni d'onda vanno intese in rap- 
presentazione di Heisenberg. L'idea e semplicemente quella di partire dalla 
funzione di Green e di inserire tra la coppia di operatori di campo l'operatore 
identita scritto come: 

I = J2\^j)^j\ , (3-27) 
j 

dove sono gli autostati dell'Hamiltoniana completa[§. 

iG(l,2)=J](*o|^(l)|*j)<*j|^(2)|*o)e(t 1 -t 2 )- 

<^o|fc(2)|* J ><* J |^(l)|W )e(t 2 - ti) . (3.28) 

Notiamo che nella sommatoria non compare il termine per J = poiche il 
valor medio sullo stato fondamentale di un operatore di campo e nullo dato 
che l'Hamiltoniana scelta commuta con l'operatore numero di particelle. Ora, 
sfruttiamo il fatto che le funzioni d'onda sono autostati dell'Hamiltoniana 



6 Abbiamo prcso la funzione d'onda dello stato fondamentale normalizzata 
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completa otteniamo: 



iG(l,2) = ^<* |^(x 1 )|*j)(W J |^(x 2 )|*o>e(t 1 -t 2 )e 



^- 1 (Sj-S )(il-t2)_ 
J 



Scrivendo infine la funzione G come 



1 /" +oc p-M*-*') 
Q(t -t') = lim / (3.29) 

2717 r?->0 UJ + 11] 



e applicando la trasformata di Fourier otteniamo: 

^o|^t(x 1 )| ^ J )(^ J |^(x 2 )|^ ) 
-uj + ir) 



G niTa {Xi,X2,U) = > — + 

t-y 1 U -Uj + 17] 



(3.30) 

UJ + UJj — If] 



dove abbiamo definito Ticjj = (Ej — Eq)- Questo modo di scrivere la funzione 
di Green nello spazio delle frequenze ci permette di individuarne i poli come 
la differenza tra l'energia dello stato fondamentale e quella di un qualsiasi 
stato che presenti una particella in piu o in meno di tale stato. II fatto 
che consideriamo solo gli stati che abbiano una particella in piu o in meno 
dello stato eccitato dipende dal fatto che stiamo considerando gli elementi di 
matrice degli operatori di campo che creano o distruggono una particella. 

La funzione di Green potrebbe costituire quindi, in principio, un metodo 
per ottenere le eccitazioni del sistema, per stati con N ± 1 particelle, a partire 
dalla DFT e in particolare dal teorema di Hohenberg-Kohn, dato che la fun- 
zione di Green e calcolata come valor medio rispetto alio stato fondamentale 
del sistema. D'altra parte, ancora una volta, il teorema di Hohenberg-Kohn 
non fornisce un metodo pratico per la costruzione della funzione d'onda e 
dunque per il calcolo dello spettro di eccitazione del sistema. Lo schema di 
Kohn e Sham non e infatti qui applicabile dato che l'espressione necessita del- 
la funzione d'onda molti corpi del sistema. In effetti non e a tutt'oggi ancora 
stato sviluppato un metodo pratico, per la DFT, che permetta di accedere alle 
energie di eccitazione che con N ± 1 particelle. Per ulteriori approfondimenti 
sull'argomento si veda la referenza [2]. 
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3.6 II teorema di Wick 
3.6.1 Introduzione al teorema 

Ora che abbiamo trovato un modo di esprimere lo stato fondamentale del 
sistema interagente a partire dallo stato fondamentale del sistema in cui w = 
vogliamo utilizzare questo risultato per calcolare la funzione di Green. Dalla 
definizione di funzione di Green otteniamo: 



(*0 


T 


'$ H m H {2) 






<*0 





D ($ \U e (+oo,0)T ^(1)^(2) t> e (0,-oo)|$ > 

lim : : " 

e->o D 



($o|E4(+oo,0)l7 6 (0,-oo)|$o) 



(3.31) 



Definiamo quindi l'operatore S E = U £ (+oc, — oo) e riscriviamo gli operatori 
di campo in rappresentazione di interazione e otteniamo: 



lim 



U e (+oo, t)4(l)U e (t u t 2 )i>i{2)U e {t 2 , -oo) |$ ) 



(%\S E \$o) 



(3.32) 



dove gli operatori unitari U e (+oo, 0) e U £ (0, — oo) sono passati all'interno del 
T-ordinamento poiche indipendenti dal tempo. Scrivendo infine esplicitamen- 
te tutti gli operatori di evoluzione temporale e facendo alcuni "giochetti con 
gli indici" otteniamo quello che puo essere definito "The most useful result of 
quantum field theory" [TJ: 



+2^ /-A" l f +OQ 
iG a p(xt x ,yt y ) ( ~f- ) — / dti...dt n 



($o|T Htih) . . .H^Qi)} a {x,t x )i> It p(y,t y ) |$ ) 



($o\S\$ ) 



iG(l,2) 



^l(l)V/(2) |$o) 



($o|5|$o) 



(3.33) 



dove nella prima riga, avendo scritto in modo esplicito l'operatore di scatte- 
ring S abbiamo usato una notazione differente per i due estremi non integrati 
della funzione di Green (x, a,t x ) e (y,f3,t y ) per evidenziarli. Notiamo in 
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particolare die nell'ultima espressione e stato svolto il limite per e — > 0, la 
cui convergenza e garantita dal fatto che gli operatori di evoluzione tempo- 
rale compaiono sia al numeratore che al numeratore, in analogia a quanto 
accade nella dimostrazione del teorema di Gell-Mann e Low, cancellando la 
divergenza della fase neU'espressione. Risulta quindi S = lim^o S e . 

Scopo del teorema di Wick e esattamente quello di valutare il T-ordinamento 
di un numero arbitrario di operatori di campo e dunque di valutare in modo 
esplicito l'ultima formula scritta. 

3.6.2 Enunciato 

Prima di scrivere l'enunciato del teorema dobbiamo introdurre due nuovi 
operatori, il primo e lW-ordinamenti (o ordinamento normale) di operatori 
di campo, mentre il secondo e la contrazione di coppie di operatori di campo. 

L 'ordinamento normale del prodotto di un numero arbitrario di operato- 
ri di campo non e altro che un riordinamento di tali operatori in cui tutti 
gli operatori di distruzione, ovvero quegli operatori che applicati alio stato 
fondamentale di un sistema danno come risultato zero, vengono spostati a 
destra. In altre parole 

= {-1) P $41M 2 , (3.34) 

dove con ip e ip^ abbiamo indicato appunto se un operatore e di distruzione 
o di creazione. P indica il numero di permutazioni eseguite per ottenere la 
formula finale in questo esempio. 

La contrazione di una coppia di operatori di campo e invece definita come 
la differenza tra il T-ordinamento e l'ordinamento normale di tale coppia di 
operatori: 

i — i 

= T[%(4] - N$$\ . (3.35) 

Poniamo l'accento sul fatto che sia il T-ordinamento che l'N-ordinamento 
sono definiti a partire dagli operatori di campo e che possono essere facilmen- 
te estesi a qualsiasi operatore poiche un qualsiasi operatore puo essere scritto 
utilizzando gli operatori di campo (anche se nel caso del T-ordinamento si in- 
contrano dei problemi dovuti al fatto che esso non e ben definito per operatori 
a tempi uguali). L'operatore di contrazione e invece definito solo per coppie 
di operatori di campo, la definizione non h valida per una coppia qualsiasi di 
operatori. 

Definiti questi operatori il teorema di Wick puo essere enunciato nella se- 
guente maniera: 



N 
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"II T-ordinamento di un numero qualsiasi di operatori di campo e uguale 
alia somma dell'ordinamento normale di tali operatori, piu la somma su tutte 
le possibili coppie deirordinamento normale di tali operatori in cui una coppia 
e pero "contratta", piu la somma su tutte le possibili doppie coppie di ope- 
ratori dell'ordinamento normale di tali operatori in cui due coppie sono pero 
"contratte", piu..." continuando cosi fino al termine in cui tutti gli operatori, 
tranne al piu uno, sono contratti. 

Tale teorema afferma in altre parole che 



i 1 i — i i — i 

T[^i4i>..] =jv[^^..]+^iv[^W..]+ ^[## •••]+••■ • ( 3 - 36 ) 

rn n n 



L'idea del teorema e quella di spostare via via verso destra tutti gli ope- 
ratori di annichilazione che poi, applicati alio stato fondamentale del sistema 
non interagente danno appunto risultato nullo. In questo modo restano dei 
termini dovuti al fatto che la commutazione di due operatori e in generale 
differente da zero. Poiche anche i termini in cui compaiono le somme sulle 
contrazioni sono iV-ordinati gli unici termini non nulli, una volta applicato il 
valor medio sullo stato fondamentale, saranno quelli in cui tutti gli operatori 
sono contratti a coppie. E dunque evidente che un espressione in cui compare 
un numero dispari di operatori di campo sara identicamente nulla. 

In definitiva cio che resta da calcolare e il prodotto di coppie di operatori 
di campo "contratte". Per calcolare il valore di una contrazione partiamo 
dalla considerazione che la contrazione di una coppia qualsiasi di operatori e 
un numero e non un operatore. Questa affermazione puo essere dimostrata 
osservando che la differenza tra il T-ordinamento e l'ordinamento normale di 
una coppia di operatori di campo o e nulla o e pari al commutatore di tali 
operatori che sono entrambi numeri [T|. Se ora applichiamo Foperazione di 
valor medio sullo stato fondamentale a entrambi i termini della defmizione di 
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operatore di contrazione otteniamo facilmente che: 

i 1 

^(l)^(2) = z</(l,2) 
i 1 

^(l)^t(2) = -i<7(l,2) 
i 1 

^(1)^(2) = 
i 1 

^(1)^(2) = , 

dove g(l, 2) la funzione di Green del sistema non interagente, ovvero l'ordine 
zero dell'espressione ricavata nella sezione precedente. Essa e infatti il valor 
medio rispetto alio stato fondamentale del sistema non interagente di una 
coppia di operatori di campo. E implicito che tutto il formalismo e stato 
sviluppato per operatori in descrizione di interazione, poiche il nostro scopo 
era quello di calcolare il valor medio di una serie di operatori in descrizione 
di interazione. 

3.7 Analizziamo la funzione di Green con il 
teorema di Wick 

3.7.1 La prima equazione di Hedin 

Abbiamo finalmente a disposizione tutti gli strumenti matematici necessari 
per studiare la funzione di Green. Partiamo innanzi tutto dall'espressione 
( I3.33P scrivendo in modo esplicito il termine: 

Hi(ti)= / d3z i rf3 z 2 w(z 1 ,z 2 )5(t 1 -t 2 )^(z 1 tt)^(z 2 ttMz 2 t 2 Mzit 1 ) 

01, (T2 

(3.38) 

A questo punto le cose diventano "abbastanza semplici" poiche non dobbiamo 
far altro che contrarre in tutti i modi possibili gli operatori di campo che com- 
paiono ad ogni ordine di tale espressione, ricordando che la contrazione di una 
coppia da come risultato la funzione di Green libera. Nell'eseguire le contra- 
zioni partiamo sempre dall'operatore di campo ip^ a (xt x ) e terminiamo sempre 
con il termine ^g(yt y ) poiche sono evidentemente gli estremi della funzione di 
Green completa e sono le uniche due grandezze non integrate nell'espressione. 
II risultato sara il prodotto di convoluzione di una serie di funzioni di G(l, 2) 
e di w(l, 2) che possono essere rappresentati con i diagrammi di Feynman. 



(3.37) 
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Tralasciamo i dettagli di questo procedimento e ci limitiamo a segnalare 
che svolgendo lo sviluppo sia per il numeratore che per il denominatore e pos- 
sible dimostrare che tutti i termini che compaiono al denominatore cancellano 
in modo esatto alcuni termini che compaiono al numeratore ed in particola- 
re vengono cancellati tutti quei termini che contengono diagrammi di vuoto, 
dove i diagrammi di vuoto sono definiti da una serie di contrazioni che parte 
e termina su due punti dello spazio N = (z N , a N , t^) e N' = (zn>, cat', ijv) 
entrambi integrati. 

Resta dunque intuitivo che a tutti gli ordini in teoria perturbativa po- 
tremo scrivere il risultato come una prima funzione di Green che parte dal 
punto (x, a,t x ) ed arriva nel punto 1 = (zi,<7i,£i), piu una serie di termini 
piu un'ultima funzione di Green che parte da un punto iV = (z N ,a^,t N ) e 
che termina nel punto (y,/3,t y ). Raccogliendo su tutti gli ordini queste due 
funzioni di Green e definendo "cio che rimane" come la self energia del sistema 
otteniamo la seguente equazione 

G(l, 2) = (/(l, 2) + g(l, 3)S H (3, 4)#(4, 2) , (3.39) 

dove stiamo integrando sulle variabili ripetute, il pedice alia self energia ci 
ricorda che il potenziale di Hartree e in essa incluso. L'ultimo passo che 
ci manca per scrivere l'equazione di Dyson e per la funzione di Green, che 
costituisce la prima equazione di Hedin, e quello di studiare il termine appena 
defmito. 

Studiando i diagrammi di Feynman e possibile mostrare che la Self-energia 
cosi defmita puo essere scritta come una serie di infmiti termini identici che 
si ripetono e quindi 

(1, 2) = E^(l, 2) + E^(l, %(3, 4)E^(4, 2)+ 

E^(l, %(3, 4)E^(4, 5)^(5, 6)S^(6, 2)... (3.40) 
= E^(1,2)+E^(1,%(3,4)Eh(4,2) . 

Per la definizione di E^(l,2) rimandiamo alia referenza PQ. Inserendo l'e- 
quazione appena scritta all'interno della (I3.39P e usando l'equazione stessa 
otteniamo finalmente: 

G(l, 2) = g(l, 2) + g{l, 3)E^(3, 4)G(4, 2) . (3.41) 

3.8 Equazione del moto con la self-energia 

Ora che abbiamo costruito l'equazione di Dyson per la funzione di Green 
possiamo cercare di guardare da un nuovo punto di vista quale sia il significato 
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fisico delle grandezze sin qui definite. A questo scopo ci serviremo, oltre che 
di tale equazione, della rappresentazione di Lehmann della funzione di Green, 
per confrontare la funzione di Green libera e quella interagente. 

Iniziamo con l'osservare che la funzione di Green libera, cioe soluzione 
dell'equazione 







S(1,2) = M(1,2) , 
puo essere facilmente scritta in rappresentazione di Lehmann come 



# CTlCT2 (xi,x 2 ,a;) = 

i 



ipi*! (xi)^? CT2 (x 2 ) V4 2 (x 2 )^ CTl (xi) 



uj — Ui + irj 



co + Ui + irj 



(3.42) 



(3.43) 



dove hui = €i — e e e, son gli autovalori di singola particella dell'Hamilto- 
niana h . Partendo dall'espressione (13.301) abbiamo utilizzato come base le 
funzioni d'onda dell'Hamiltoniana non interagente e abbiamo quindi sfruttato 
il fatto che esse sono il prodotto di funzioni d'onda a particella singola. Ve- 
diamo quindi che la funzione di Green non interagente e facilmente calcolabile 
(poiche e possibile calcolare le funzioni d'onda di particelle non interagenti) ed 
e piu immediata la sua interpretazione fisica: i suoi poli sono le auto-energie 
di singola particella del sistema e a numeratore troviamo il prodotto delle 
funzioni d'onda e quindi la matrice densita del sistema. 

Cio che vogliamo fare ora e cercare di ottenere una formulazione simile 
per la funzione di Green completa. A questo scopo applichiamo l'operatore 

d - __ 

ih- /^o( x i) all'equazione di Dyson della funzione di Green completa f)3.4ip 

ot 

ottenendo: 



ih~T- - ho(x.i) 
ot 



G(l,2)= M(l,2) + ^(1,3)G(3,2) 



(3.44) 



dunque 







(^--/ i0 (x 3 ))5(l,3)-^(l,3) 



G(3,2) = M(l,2) 



(3.45) 



dove abbiamo utilizzato l'equazione (I3.42p . 

E possibile dimostrare [3 Ej che se inseriamo l'espressione di Lehmann 
(I3.30p nella (I3.45P e riscriviamo il numeratore definendo le ampiezze di Leh- 
mann come: 

' {Eo\^{x)\Ej) if coj>fi 
(Ej|^(x)|E ) Mu)j<n 



(3.46) 
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otteniamo un'equazione del tutto simile a quella per il caso non interagente, 



(hcuj - h (xa))6(l, 3) - hT? H (l, 3; u>j) f J(T3 (x 3 ) = 



(3.47) 



con la differenza che ora anziche il prodotto delle funzioni d'onda di particella 
singola abbiamo il prodotto delle funzioni /j CT (x) e che questa volta compare 
la self energia del sistema. huj = (Ej — E ) risulta essere l'autostato di queste 
funzioni rispetto all'equazione di Schroedinger appena scritta e coincide con 
i poli della funzione di Green. Nelle ultime due equazioni e stata inoltre 
introdotta la notazione 1 = (xi,<Ti) in cui, rispetto alia notazione abituale, 
non e piu inclusa la variabile temporale. 



3.9 Le equazioni di Hedin 

Lo studio della funzione di Green tramite lo sviluppo della self-energia a 
ordini successivi sin qui presentato porta alia formulazione di equazioni che 
non possono essere risolte se non nel caso del gas uniforme di elettroni poiche 
si ottengono delle espressioni molto complesse dal punto di vista matematico. 
Anche in questo caso pero la soluzione puo essere trovata solo per i primi 
ordini perturbativi e di nuovo e tutt'altro che banale. Un ulteriore problema e 
costituito dal fatto che abbiamo scelto di scrivere come correzioni perturbative 
delle quantita che non siamo sicuri siano trascurabili o piccole rispetto ad altri. 
Cio si traduce in quella che puo essere considerata una delle maggiori difficolta 
di tutte le teorie di campo, quale e la MBPT, e che e costituita dalla presenza 
di termini di ordine elevato che, anziche che essere trascurabili, risultano 
essere divergenti. Un'esempio e costituito dal secondo ordine dell'energia del 
gas uniforme appunto e piu in particolare dal termine che viene solitamente 
indicato [T] come E™ 9 . 

Un approccio alternativo a quello sin qui introdotto, sempre basato sulle 
grandezze della funzione di Green e della self-energia e quello delle cinque 
equazioni di Hedin. L'idea e quella di trovare delle quantita che insieme 
alia funzione di Green e alia self-energia formino un set di equazioni chiuso. 
Tale set di equazioni ancora una volta risulta essere troppo complesso per 
essere risolto in modo esatto ma ci sono tentativi di risoluzione, attraverso 
l'uso dei calcolatori, in modo iterative L'idea per superare il problema dei 
termini divergenti che compaiono nella MBPT e quella di non scegliere come 
"grandezze fondamentali" quelle libere, #(1,2) e w(l,2), ma le grandezze 
"rivestite" G(l,2) e W(l,2). 

Prima di procedere con la formulazione del set di equazioni di Hedin segna- 
liamo che il punto di partenza sara Fespressione gia scritta per la funzione di 
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Green grazie all'ausilio del teorema di Gell-Mann e Low. In questo caso pero 
noteremo come, a differenza si quanto sin qui fatto, l'operatore S(+oo, — oo) 
sia costruito a partire da un potenziale esterno anziche dall'interazione w. 
Questo poiche utilizzeremo il teorema di Gell-Mann e Low pensando l'Hamil- 
toniana separata in due parti con Hq costituito dall'Hamiltoniana completa 
(interazione inclusa) e Hi — <p, che rappresenta un potenziale esterno appli- 
cato al sistema. Ovviamente la funzione di Green cosi definita sara differente 
da quella sinora definita, ritroveremo la funzione di Green sin qui studiata 
semplicemente ponendo il potenziale esterno pari a zero. 



3.9.1 Introduzione alle equazioni 

Partiamo come annunciato con lo scrivere la funzione di G+reen secondo la 
formula (I3.33|) 



iG v {±,2) 

dove l'operatore S e defmito come 
T[S] = T 



(* |T ^t(l)^( 2 ) |* ) 



(^ol^l^o) 



(3.48) 



exp ( -- V f (l + )V(1)^(1) 



(3.49) 



e lo stato rispetto a cui facciamo il valor medio e questa volta lo stato 
fondamentale dell'Hamiltoniana completa. 

Eseguiamo ora alcuni "giochetti" che ci permettono di scrivere una rela- 
zione tra la funzione di Green a particella singola e la funzione di Green a 
due particelle; scriviamo innanzitutto la variazione della funzione di Green 
cosi definita rispetto al potenziale esterno (p 

^G v (l,2) = L ; J iG v (l,2 y*?l II I] 01 , (3.50) 



(*o|S|*o> 



<^o|5|^ ) 



dove abbiamo utilizzato 
T[SS] -- 



(3.51) 



s ^t(i+Mi)Mi; 

quindi utilizziamo la definizione di funzione di Green a due particelle ottenedo 
5GV(1,2) 



M3) 



-h^GJl, 3; 2, 3+) + h^GJl, 2)G V (3, 3" 



(3.52) 
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Possiamo finalmente valutare la nostra espressione nel caso in cui il poten- 
ziale esterno sia uguale a zero in modo da avere delle grandezze confrontabili 
con i risultati sin qui ottenuti per via perturbativa0. In particolare sostituendo 
l'espressione ricavata per la funzione di Green a due particelle nell'equazione 
del moto per la funzione di Green (I3.26P otteniamo: 



itifo ~ ^o(xi) 



G(l,2) + i u>(l,3)G(3,3 + )G(l,2) 



ihw(l,3) S -^± = h6(l,2) . (3.53) 



Utilizziamo ora la seguente identita matematica per riscrivere l'equazione 
ottenuta: 

e confrontandone il risultato con l'equazione ( I3.45P otteniamo 

^(1,2) = -t5(l,2) w(l,3)G(3,3 + )+ 

-»^1,3)G(1,4) ^" 1(4 ; 2) (3.55) 

= 5(1, 2)v H (l) + z w(l, 3)G(1, 4)r(4, 2; 3) . 

II primo termine e evidentemente il potenziale di Hartree, mentre il secondo 
termine costituisce un'espressione in forma funzionale (e non piu perturbativa 
quindi) dei termini della self-energia oltre il potenziale di Hartree. Nell'ultima 
linea abbiamo gia definito la funzione di vertice come 

Potremmo ora procedere con il costruire un'equazione di Dyson per la fun- 
zione vertice e ottenere cosi un set di 3 equazioni con tre funzioni incognite: 
il propagatore, la self-energia e il vertice (l'interazione libera e nota). Cio che 
invece si preferisce fare e partire dall'interazione "rivestita" anziche da quella 
libera e costruire il set di cinque equazioni che dimostreremo nel prossimo 
paragrafo poiche, come gia detto, l'utilizzo dell'mterazione libera porta dei 
problemi di divergenza nella teoria. 



7 D'ora in poi tutte le funzioni di Green saranno intese valutate a potenziale esterno 
nullo c quando indichcrcmo una variazione rispetto al potenziale sottintenderemo di fare 
la variazione e poi valutare l'espressione a campo nullo 
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3.9.2 Dimostrazione delle equazioni 

In questo paragrafo presenteremo la teoria cosi come l'abbiamo trovata in 
letteratura. Poiche d'altra parte crediamo che questo modo di procedere 
(sebbene elegante dal punto di vista matematico) impedisca in qualche mo- 
do di vedere la fisica che dietro di essa e nascosta, nel paragrafo successivo 
cercheremo di svolgere quelle considerazioni a nostro avviso necessarie per 
capire il significato delle grandezze che compaiono. Come altrove nell'esposi- 
zione della MBPT saremo dunque costretti a presentare la teoria in modo non 
lineare, crediamo d'altronde che sia questo il modo migliore di guardare ad 
essa e che entrambi i livelli di analisi siano necessari per poterla comprendere. 

Iniziamo dunque con il defmire un "potenziale esterno efficace" pari alia 
somma della perturbazione introdotta piu il potenziale di Hartree 

V(l) = <p(l)+v H (l) ■ (3.57) 

Definiamo quindi la funzione dielettrica come (vedremo piu avanti il perche 
di questa scelta) 

1 ' ] 6<p(2) (3.58) 
= 5(1,2)+ 3)11(3, 2) , 

dove per la seconda linea abbiamo utilizzato la defmizione classica di funzione 
di risposta, ricordando che la densita e pari alia parte diagonale della funzione 
di Green 

, . SG(1,1+) 

n(1 - 2 '=- < -^r ' (3 - 59) 

Otteniamo quindi l'interazione efficace come 

W(l,2) = e- 1 (1,3M3,2) 

= w(l,2)+ w(l, 3)11(3, 4)w(4, 2) . [ ' } 

Utilizziamo ora le regole di derivazione a catena nell'equazione scritta per la 
funzione risposta ()3.59|) . o polarizzazione del sistema, e scriviamo: 

nn 2)= _j *C(M + )*V(3) 
1 ' ) SV(3) S(p(2) (3.61) 

= 6(1,2) + n(l,3)w(3,4)n(4,2) , 



dove abbiamo definito la polarizzazione ridotta come la derivata della densita 
rispetto al potenziale efficace; in questo modo abbiamo "messo in evidenza" il 
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ruolo del campo di Hartree. Inserendo l'ultima equazione scritta nella f)3.60p 
otteniamo 

W(l, 2) = w(l, 2) + w(l, 3)n(3, A)W(A, 2) (3.62) 

Utilizzando ora nuovamente la relazione f)3.54p insieme alle due definizioni 
di funzione risposta e funzione risposta ridotta otteniamo 

HU(1, 2) = - * G(l, 3)G(4, l + )r(2; 3, 4) 
ftn(l, 2) = -iG(l, 3)G(4, l + )f (2; 3, 4) . 

Le due grandezze differiscono solo per il fatto che una e considerata rispetto 
al potenziale libero mentre l'altra rispetto al potenziale efficace. Alio stesso 
modo differiscono la funzione T e la funzione T che abbiamo qui introdotto. 

Se ora riprendiamo l'equazione per la self-energia scritta nel paragrafo 
precedente e utilizziamo la regola di derivazione a catena a la definizione di 
interazione efficace appena introdotta otteniamcH 

= zG(l,4)W(3,l+)r(4,2;3) . 

Siamo cosi quasi giunti alia fine dei nostri sforzi. Ci manca ora solo di 
ricavare l'equazione per il vertice. Questa puo essere ricavata partendo dal- 
la definizione di funzione di vertice, utilizzando l'inverso dell'equazione del 
propagatore 

^(1, 2) = G~\\, 2) + 5(1, 2)(V H (1) + ^(1)) + £*(1, 2) , (3.64) 

dove abbiamo esplicitato nuovamente il potenziale esterno per maggior chia- 
rezza, nell'equazione per il vertice si deve poi immaginare, come sempre, di 
mettere il potenziale pari a zero. Utilizziamo ora le regole di derivazione a 
catena insieme con la relazione (]3.54p . In questo modo si ottiene (scriviamo 
qui il risultato per la funzione ridotta, ma si potrebbe fare lo stesso con la 
funzione T completa) 



f(l; 2, 3) = 5(1, 2)6(1, 3) 



E*(2,3) 



5V(1) 

5(1,2)5(1,3)+ (3.65) 
If(l;4,5) G (6,4) G (7,5)§!ff ' 



8 Indichiamo qui con £* solo il secondo termine della self energia, ovvero la self-energia 
meno il potenziale di Hartree Y,* H — Vh, per questa ragione non comparira piu, d'ora in 
poi, il pedice H accanto all'operatore della self-energia. 



3.9 Le equazioni di Hedin 



43 



Scriviamo dunque le equazioni che abbiamo sin qui ricavato, introducendo 
per la scrittura finale la funzione di Green di Hartree in modo da utilizzare 
solo il secondo termine della come gia fatto in precedenza. 

g] I 1 (l,2)=g- 1 (l,2) + V H (l,2) . (3.66) 
Le equazioni di Hedin risultano quindi cosi scritte: 

G(l, 2) = g H (l, 2) + g H (l, 3)£*(3, 4)G(4, 2) (3.67) 

W{\, 2) = w(l t 2) + w(l, 3)6(3, 4)W(4, 2) (3.68) 

nfi(l, 2) = — if (1; 3, 4)G(2, 3)G(4, 2) (3.69) 

ft£*(l, 2) = zT(4; 1, 3)G(3, 2)W(2, 4) (3.70) 
f(l;2,3) = 5(1,2)5(1,3) + 

^- 1 f(l;4,5)G(6,4)G(5,7)^|^ . (3.71) 

3.9.3 La fisica dietro le equazioni di Hedin 

Prima di procedere oltre con l'esposizione della teoria sin qui esposta vorrei 
soffermarmi ad osservare le grandezze che abbiamo definito e cercare di capire 
quale sia il loro significato fisico. A questo scopo faremo uso, in modo decisivo 
delle regole di cambio base in seconda quantizzazione per un operatore e del 
concetto stesso di operatore che abbiamo precisato in appendice. 

La prima domanda che mi sono posto dopo aver studiato questa teoria 
e perche la polarizzazione (e di conseguenza la funzione dielettrica) e l'inte- 
razione "rivestita" cosi definite dovrebbero coincidere con la loro defmizione 
nell'elettromagnetismo classico. A questo scopo osserviamo piu con calma le 
grandezze che compaiono. Innanzi tutto abbiamo un sistema di particelle che 
interagiscono e l'operatore che rappresenta l'interazione e il seguente 



w= / rf 3 xrf 3 y ^ t (x)^ t (yV(x,yMy)^(x) , (3.72) 




che e evidentemente un operatore a due particelle, dove iy(x, y) = . II 
campo medio e il potenziale di Hartree ed e rappresentato dall'operatore ad 
una particella 

v H = J J d 3 xd 3 y ^t( x ) ^(y) ^( x ) . (3.73) 

Ipotizziamo infine, cosi come facciamo nella costruzione delle equazioni di 
Hedin, che sia presente un potenziale esterno applicato a tale sistema: 

<p= [ rf 3 x <^(x)^ t (x)^(x) . (3.74) 
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Immaginiamo ora che il sistema si trovi nel suo stato fondamentale e di ap- 
plicare una variazione del campo esterno 5ip(x) il campo di cui risentira una 
particella qualsiasi del sistema sara dovuto alia somma del campo applica- 
to piu la variazione dell'interazione con tutte le altre particelle dovuta alia 
variazione della funzione d'onda delle stesse a causa del campo applicato. 

6V eff = 60 + 5w . (3.75) 

Dove l'operatore Sw deve rappresentare in qualche modo il campo dovuto alia 
variazione delle funzioni d'onda delle particelle^]. 

Cioe che nelle equazioni di Hedin e approssimare tale operatore con la 
variazione del campo di Hartree. Possiamo guardare all'approssimazione fatta 
come ad un modo di trascurare gli effetti quantistici di scambio e correlazione 
e tenere solo il campo medio. Se ora scriviamo esplicitamente la variazione 
del potenziale di Hartree 

8v H = J J d 3 xd 3 y ^t( x ) j( x ) 9 (3.76) 

vediamo facilmente come la funzione di risposta, che lega la variazione di 
densita al potenziale esterno all'ordine lineare rispetto alia densita, possa 
essere anche utilizzata per definire l'interazione rivestita e quindi la funzione 
dielettrica. Scriviamo in modo esplicito quanto affermato: 

8V e ff ~ 50+ w(l,2)x(2, 3)50(3) , (3.77) 

da cui si capisce il perche della funzione dielettrica definita come la variazione 
del potenziale efficace rispetto al potenziale esterno. 

Abbiamo parlato di approssimazione per l'interazione e abbiamo chiarito 
in che senso l'interazione e stata considerata in modo approssimato. Questo 
non ci deve pero in alcun modo far pensare che la soluzione delle equazioni 
di Hedin sia in qualche modo approssimata! Tali equazioni sono infatti state 
costruite per permetterci di calcolare la funzione di Green e la sua self-energia 
e sono queste le due grandezze a cui dobbiamo guardare. L'interazione che 
compare (e con essa la polarizzazione) sono solo un mezzo per ottenere questo 



Formalmente non ha senso parlare di variazione di un operatore poiche esso e definito in 
modo unico e non dipende dalla funzione d'onda. Quello che viene fatto in pratica e cercare 
di costruire un operatore che possa tener conto della variazione delle funzioni d'onda del 
sistema e che risulti essere un operatore a particella singola. Nella presente teoria si riesce 
a includere solo una parte di questo effetto, il campo medio; se ad esempio guardiamo 
alia DFT in questa teoria tale effetto e viene incluso in modo satto per il sistema non 
interagcnte come variazione del potenziale di scambio e correlazione. II prczzo da pagare e 
pero che in questo caso la funzione d'onda molti-corpi non sara quella csatta. 



3.9 Le equazioni di Hedin 
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scopo. L'idea di utilizzare un'interazione COS1 fclttcl nasce dalla necessita di 
una grandezza fisica che possa essere un punto di partenza migliore dell'in- 
terazione libera, ma che non deve necessariamente essere esatta. Gli effetti 
quantistici che non compaiono in tale grandezza saranno quindi inclusi altrove 
ed in particolare questo ci porta a guardare alia funzione di vertice. 

L'interpretazione fisica di questa grandezza appare meno immediata delle 
altre (perlomeno all'autore di questa tesi) ma puo essere in qualche modo 
delineata perlomeno in modo intuitivo. L'idea e che stiamo studiando un 
sistema in cui compaiono due grandezza, il campo fermionico degli elettroni 
e l'interazione. Per entrambe abbiamo definito "un propagatore" (possiamo 
immaginare l'interazione libera come il propagatore del fotone, sebbene in 
realta il campo venga trattato qui in modo classico) e la sua self-energia (per 
l'interazione essa sara costituita dalla polarizzazione). Appare quindi naturale 
la presenza di un termine che "leghi" tali "campi" e che, come detto, dovra 
in qualche modo supplire al fatto che stiamo studiando il problema in modo 
semi-classico (o perlomeno con l'approssimazione appena indicata per quanto 
concerne l'interazione). Per l'interpretazione della funzione di vertice come 
inclusione degli effetti quantistici trascurati nella definizione dell'interazione 
rivestita rimandiamo alia bibliografia [U [5] . 
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Approssimazioni 



4.1 Introduzione 

Fino ad ora abbiamo presentato la formulazione esatta di due teorie per lo 
stato fondamentale di un sistema di elettroni interagenti. Se pero vogliamo 
applicarle per lo studio di sistemi fisici dobbiamo necessariamente introdurre 
alcune approssimazioni. La necessita di eseguire tali approssimazioni nasce 
in modo differente per le due teorie: nel caso della DFT la necessita nasce 
dal fatto che non e noto il potenziale di scambio e correlazione e che dun- 
que dovremo in qualche modo trovare una buona approssimazione per tale 
potenziale. Per la MBPT il problema nasce invece dal fatto che i calcoli che 
compaiono sono troppo complessi per poter essere risolti in modo esatto. 

Le due approssimazioni che solitamente vengono utilizzate sono l'appros- 
simazione locale o "Local Density Approximation" (LDA) per la DFT e l'ap- 
prossimazione "GW" per la MBPT. La prima consiste nello scrivere il fun- 
zionale dell'energia del sistema a partire dal funzionale del gas omogeneo che 
puo essere calcolato a partire dalla MBPT e da calcoli Montecarlo, la secon- 
da nel risolvere le equazioni di Hedin ( I3.67P approssimando la funzione T con 
l'identita e dunque la Self-Energia con l'espressione £ = GW, da cui il nome 
dell'approssimazione. 

Nei prossimi paragrafi ci concentreremo soprattutto sulle approssimazio- 
ni necessarie all'utilizzo della DFT poiche i calcoli che faremo all'interno di 
questa tesi saranno svolti nell'ambito di tale teoria. La teoria MBPT restera 
comunque sempre presente sia come uno dei possibili punti di partenza per 
la costruzione di approssimazioni per i funzionali del potenziale e dell'energia 
della DFT sia, come vedremo nell'ultimo capitolo di questa tesi, come pos- 
sibile punto di partenza per correzioni alle approssimazioni svolte nei calcoli 
da noi eseguiti. 
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4.2 L'approssimazione locale 
4.2.1 Approssimazione LDA 

L'approssimazione locale costituisce una notevole semplificazione della teoria 
presentata e si basa sull'assunzione che la densita di energia in un punto dato 
del nostro sistema dipenda solamente dalla densita del gas elettronico in quel 
punto, da cui il nome di approssimazione locale. Facendo questo presupposto 
e possibile approssimare l'espressione dell'energia di scambio e correlazione 
come funzionale della densita nel seguente modo 

EM = J dh p(r)e xc [p(r)} * j d 3 v p{v)e^{p{v)) , (4.1) 

dove e^° m (p) e una semplice funzione di una variabile che puo essere ricavata 
dall'energia del sistema con densita uniforme p 

e h x ° c m (p) = - £ yf I - , (4-2) 

dove V e il volume del sistema. II potenziale di scambio e correlazione in 
questo modo risulta 

^■ c H(r) = ^^-4r(p(r)) + p(r)^f^ . (4.3) 
5p{r) dp{r) 

Per esprimere il funzionale dobbiamo quindi scrivere qual'e l'energia del 
sistema omogeneo in funzione della densita. A questo scopo scegliamo in- 
nanzi tutto di separare l'energia di scambio e correlazione in due parti, un 
termine per l'energia di scambio ed un termine per l'energia di correlazione 
appunto. II primo termine di approssimazione locale puo essere ottenuto a 
partire dall'espressione dell'energia di scambio dell'interazione che sappiamo 
puo essere scritto in modo esatto in funzione della matrice densita di Kohn e 
Sham del sistema [2J. 

E? a \p] = 4 / d 3 r I 'dh'^^f . (4.4) 



dove 7 S (r, r') = (r)^^')]^ ) e appunto la matrice densita del sistem£0- 

Sviluppando infatti all'ordine lineare la matrice densita rispetto alia densita 



1 Vcdiamo come la matrice densita diffcrisca dalla funzione di Green poiche quest 'ultima 
c definita cpon gli operatori in rapprcsentazione di Heisenberg e quindi c dipendente dal 
tempo. 



4.2 L'approssimazione locale 
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del sistema si ottiene [2] 

e« S = -dopT con ^o = ^0) V3 • (4.5) 

Definendo la quantita = j|- (cioe il raggio della sfera che contiene una 
quantita unitaria di materia), possiamo riscrivere l'espressione come [6]: 

0.458 e 2 . . 

e x = • 4.6 

r s 

Per la parte di correlazione una delle prime parametrizzazioni proposte in 
questa approssimazione e invece la seguente [5|: 

44 e 2 

(4.7) 



r s + 7.8 a 

dove a e il raggio di Bohr. Abbiamo dunque 

0.458 e 2 0.44 e 2 



r s + 7.8 a 



Molte altre parametrizzazioni piu sofisticate sono comunemente note pQ e 
sono ottenute da fits di calcoli ab-initio virtualmente esatti. 

4.2.2 Approssimazione LSDA 

Fino ad ora abbiamo parlato di parametrizzazioni per il funzionale E xc [p\; al- 
l'interno di questa tesi faremo pero principalmente uso del funzionale E xc [p^, pj 
per lo studio di sistemi spin polarizzati. II funzionale E xc [p], e dunque anche 
la sua approssimazione LDA, tratta in modo simmetrico la componente e 
la componente pdownarrow della densita e dunque se provassimo a costruire da 
questa i potenziali di scambio e correlazione vl c [p] e v^ c [p] otterremmo 

t r I £%Ec[Pt + gj SE xc \p] 5p 

vUp] = ^^=—Jir v - [p] (49) 

v i ln] _ SE xc [ Pt + Pi ] _ 5E xc [p] 5p _ 

Appare dunque evidente che per sistemi con magnetizzazione di spin sia 
necessario introdurre uno schema differente per approssimare il funzionale. 
L'idea e quella di partire dall'energia del sistema omogeneo completamente 
polarizzato, dove per completamente polarizzato intendiamo che la grandezza 

C = (4.10) 
P 
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sia pari ad 1, il die significa die tutti gli elettroni hanno spin orientato nella 
stessa direzione. Calcolata l'energia del sistema omogeneo imponendo ( = 1, 
si trovano, in approssimazione LDA, le seguenti parametrizzazioni 



e xc (p, = e xc (p, C = 0) + (e xc (p, C = 1) - e xc {p, ( = 0))g(() , (4.12) 



4.3 Perche l'approssimazione LDA funziona 
cosi bene 

4.3.1 Introduzione 

L'approssimazione eseguita (LDA) e basata sull'idea che il funzionale del- 
l'energia di scambio e correlazione non differisca di molto, nella forma, da 
quello del sistema omogeneo. Ci si aspetterebbe dunque die un approccio di 
questo tipo possa esser ragionevole per sistemi in cui la densita varia lenta- 
mente. In realta applicazioni di questa approssimazione sono state utilizzate, 
con ottimi risultati, anche in sistemi che sono molto lontani dall'avere un 
comportamento della densita lentamente variabile, come possono essere ad 
esempio le molecole che studieremo all'interno di questa tesi. 

II successo ottenuto puo comunque essere spiegato dal fatto che l'approssi- 
mazione locale soddisfa una regola di somma che riguarda "la densita di buca 
di scambio e correlazione" e che la stessa densita di buca in approssimazione 
LDA, pur essendo di molto differente da quella reale, presenta una media sfe- 
rica molto simile a quest'ultima. Nei prossimi paragrafi spiegheremo dunque 
piu in dettaglio questi punti, partendo dal teorema di connessione adiabatica 
e utilizzando nella dimostrazione il teorema di Hellmann Feynman [T5] . 

Concludiamo questa introduzione segnalando pero che pur essendo un ot- 
timo punto di partenza anche per lo studio di sistemi isolati l'approssimazione 




(4.11) 



Local 



dove g(() e una funzione di interpolazione che vale 




(4.13) 



2 o per lo studio di sistemi in presenza di campo magnctico 



4.3 Perche l'approssimazione LDA funziona cosi bene 
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LDA ha l'inconveniente di non presentare il corretto andamento asintotico del 
potenziale di scambio e correlazione; questo inconveniente sara una delle fonti 
di problemi nel calcolo dello spettro di eccitazione di alcune molecole. 



4.3.2 II teorema di connessione adiabatica 

Partiamo dal considerare un'Hamiltoniana H\ cosi costruita 



Hx = f + Y,j dh ^(r)^(r)VUr) 



/Wr'w^i'lgWftMiW^M , (4.14) 



A 



e per il valore A = 1 scegliamo di prendere v\ = V ext , il potenziale ester no 
applicato al nostro sistema fisico. 

Ora sfruttiamo il teorema di Hohenberg-Kohn per il sistema descritto dalla 
seguente Hamiltoniana e fissiamo idealmente il potenziale v\ per tutti i valori 
di A differenti da 1 in modo che risulti 



r) = E^Al^W^(r)|vI/ A ) = p 1 (r) = p o (r) , (4.15) 



In questo modo il potenziale vo e evidentemente vo = vks = V ext + vh + v xc - 
Utilizziamo ora il teorema di Hellmann-Feynman, che corrisponde alia 
prima riga della (I4.16p . e scriviamo la variazione dell'energia dello stato fon- 
damentale del sistema, intesa come il valor medio rispetto alia funzione d'onda 
dell'Hamiltoniana H\, come 



= \J J d 3 rdh> w(v,r'){^\p(v)p(v') - 6{r - r')p(rM) 

d f 

(4.16) 

L'energia dello stato fondamentale del sistema fisico, che e quella dell'Ha- 
miltoniana a A = 1, puo ora essere scritta come 

E{1) = E(0) + J d\^^- = T s [p] + j dh V ext {v)p{v) 

+ l -J J d 3 rdh' w(r,r')p(r)p(r') + E xc [p] , (4.17) 
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dove abbiamo definite- 



E *c=\j dxj j d 3 rcZV w(r,r') 

m\ ip(r) - p(r)) (p(r') - p(r')) |*°) - p(r)5(v - r')] . (4.18) 

Definendo ora la funzione di correlazione a coppie g\ come estensione della 
(12.101) nel caso in cui la funzione d'onda dipenda da A otteniamo 

E *c = \j j dhdh' w(r, t')p(t) P (t') J d\[g x (r, r', [p]) - 1] . (4.19) 

4.3.3 La buca di scambio e correlazione 

Per capire il successo dell'approssimazione LDA ci resta ora da capire cosa 
significhi eseguire tale approssimazione per la funzione g\ 
A questo scopo definiamo un'ulteriore grandezza 



MO 



[ l d\[g x (r,v', [p])-l] 
Jo 



(4.20) 



che puo essere interpretata come la densita di una buca di scambio e correla- 
zione [2]. Inserendo questa espressione nel funzionale per l'energia di scambio 
e correlazione e dunque possibile interpretare tale energia come un'interazione 
particella-buca di scambio e correlazione. 

Ora e possibile mostrare che, espandendo tale densita in termini di multi- 
polo rispetto alia distanza r — r' solo il termine di monopolo, cioe il termine 
con andamento di tipo sferico, da contributo all'energia di scambio e corre- 
lazione^) e che e sempre il termine di monopolo il solo a dare contributo alia 
regola di somma 

' 'd 3 r' &c (r,r') = -l . (4.21) 

II successo dell'approssimazione LDA nasce proprio dal fatto che anche 
per quest'ultima g xc ha un andamento di tipo sferico e soddisfa alia regola 
di somma valida per il caso esatto. In altre parole pur essendo q xc molto 
differente da g^ A per sistemi in cui la densita p varia rapidamente, risultano 
molto simili (14. ip le due grandezze g^ A e g®° c , dove con 00 abbiamo preso il 
termine di monopolo della buca di scambio e correlazione. 



3 Questo e vero solo nel caso di un'interazione che dipende soltanto da |r — r'|. 



4.4 Oltre l'approssimazione locale 
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(b) Termine di monopolo dclla buca di scambio e correlazione. Media sfcrica. 



Figura 4.1: Confronto tra l'andamento della densita esatta di scambio e 
correlazione e la sua media sferica. 



4.4 Oltre l'approssimazione locale 
4.4.1 L'equazione di Sham-Schluter 

L'equazione di Sham-Schluter costituisce un possibile punto di partenza per 
la costruzione di un potenziale di scambio e correlazione nella teoria DFT 
a partire dalla MBPT. Essa stabilisce infatti una relazione esatta tra il po- 
tenziale di scambio e correlazione e la self-energia del sistema. II punto di 
partenza per giungere a tale equazione e dato dall'osservazione che la parte 
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diagonale della funzione di Green coincide con la densita del sistema nello 
stato fondamentale che, come sappiamo, puo essere calcolata nel formalismo 
del funzionale densita. 

-iG(l,l + )=p(l) . (4.22) 

Partendo da questa osservazione l'idea e quella di definire la funzione 
di Green soluzione dell'equazione del moto di KS come GksO-i 2). Se ora 
cerchiamo di stabilire un legame tra la funzione di Green esatta e quella 
appena definita osservando le due equazioni del moto a cui queste soddisfano 
otteniamo 



G(l, 2) = G KS (1, 2) + G KS (1, 3) (£(3, 4) - 5(3, A)v xc (A) j G(A, 2) . (4.23) 

Ricordiamo che della funzione di Green GksO-^) s °1° ^ a parte diagonale ha 
significato fisico. Utilizzando ora proprio il fatto che il termine diagonale di 
entrambe coincide con la densita del sistema otteniamo 

G KS (1, 3) (V(3, 4) - 5(3, 4K C (4)^ G(4, 1) = . (4.24) 

Quella scritta e l'equazione di Sham-Schluter. Come sempre quando ab- 
biamo a che fare con grandezze a molti corpi essa non puo essere risolta in 
modo esatto, e necessario quindi trovarne delle soluzioni approssimate. Una 
possibilita e quella di trovare una soluzione linearizzata di tale equazione (il 
che significa approssimare nella (I4.24|) G(4, 1) ~ G^s(4, 1)); soluzione che 
puo essere utilizzata ad esempio per ottenere un potenziale di scambio e cor- 
relazione con il corretto andamento asintotico in modo da correggere l'errore 
prodotto dall'approssimazione LDA [S]. 



Capitolo 5 



Spettro di eccitazione. 
L'accoppiamento particella 
buca. 

5.1 Introduzione 

Fino ad ora ci siamo occupati di delineare le due teorie per lo studio dello sta- 
te fondamentale del nostro sistema ed abbiamo mostrato come da entrambe 
possiamo ottenere l'energia totale. Abbiamo inoltre gia parlato delle ener- 
gie di eccitazione del sistema per stati che presentano una particella in piu 
o in meno dello stato fondamentale del sistemi. Non abbiamo invece detto 
nulla sullo spettro di eccitazione che riguarda stati che presentino lo stesso 
numero di particelle dello stato fondamentale poiche lo studio di tali energie 
comporta la necessita di tener conto di un importante effetto fisico, l'accop- 
piamento particella buca, che le teorie si qui formulate non sono in grado 
di descrivere in modo esauriente. Avremo dunque la necessita di sviluppar- 
le ulteriormente entrambe; andremo dunque verso quelle che sono le attuali 
frontiere delle ricerca, facendo in particolare riferimento, per quanto concerne 
lo stato dell'arte alle referenze [7J |3]. 

Per quel che riguarda la teoria many body lo spettro di eccitazione puo in 
realta gia essere ottenuto a partire dalla rappresentazione di Lehmann della 
funzione risposta (o polarizzazione) . Studiando pero tale grandezza vedremo 
che compariranno delle funzioni "a quattro punti" , che sono legate quindi alia 
funzione di Green a due particelle. Per questo motivo, se vogliamo scrivere 
un'equazione del tipo "equazione di Dyson" per questa grandezza avremo 
bisogno di definire una funzione di risposta a quattro punti. II risultato di 
questo tipo di approccio sara l'equazione di Bethe-Salpeter. 
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Per quel che riguarda la teoria del funzionale densita, sebbene in principio 
a partire dal teorema di Hohenberg-Kohn potremmo accedere alio spettro 
di eccitazione del sistema, di fatto abbiamo bisogno di costruire un metodo 
pratico per accedere a tale spettro. Come gia anticipato tale metodo sara 
costituito dalla TDDFT o meglio dalla Time-Dependent Density-Functional 
Response Theory (TD-DFRTq). Dovremo dunque dimostrare nuovi teoremi 
che possano estendere la teoria al caso dipendente dal tempo. 

Prima di procedere pero alio sviluppo delle due teorie vogliamo intro- 
durre la teoria della risposta lineare e mostrare come, a partire da questa, 
sara possibile ricavare lo spettro di eccitazione del sistema attraverso la sua 
rappresentazione di Lehmann. 

5.2 La teoria della risposta lineare 

Consideriamo un sistema descritto da un'Hamiltoniana H indipendente dal 
tempo. Immaginiamo quindi di perturbare il sistema accendendo, all'istante 
t = t Q > 0, una perturbazione esterna dipendente dal tempo 5V ext (t) che 
possa essere considerata piccola. Scopo della teoria della risposta lineare e 
studiare come varia una grandezza fisica, aH'ordine lineare, dopo Faccensio- 
ne della perturbazione esterna. L'evoluzione di uno stato in descrizione di 
Schroedinger sara data dalla solita funzione esponenziale per la parte di Ha- 
miltoniana indipendente dal tempo e dall'operatore descritto nella (I3.16P per 
la perturbazione che e dipendente dal tempo. Risultera dunque: 

\*s(t)) = e- iAt / n U(t,t )Mt )) 

= e - iAt/h \y(t )) + e- lflt/h [ di! 5V^ xt (t')\^(t )) + 

dove abbiamo sviluppato all'ordine lineare rispetto alia perturbazione esterna 
l'operatore di evoluzione temporale U(t, 0). t e l'istante a cui scegliamo di 
accendere la perturbazione esterna, mentre lo stato \I/ (to) e autostato del- 
l'Hamiltoniana H ed b indipendente dal tempo: \&(to) = ^(0). Scriviamo 
ora il valor medio di una grandezza qualsiasi che venga rappresentata da un 
operatore 0(t) rispetto alio stato perturbato 

<d(*)> = mo)\o 6 (t)Mo))+ [ t dt , (m\[dH(t),5vf(t')\m) + --- , 

J t 

(5.2) 

^^Per ragioni pratiche d'ora in avanti ci limiteremo a parlarc di TDDFT 



(5.1) 



5.2 La teoria della risposta lineare 
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dove abbiamo spostato l'operatore di evoluzione temporale sull'operatore 0(t) 
in modo che ora l'operatore sia scritto in descrizione di Heisenberg rispetto 
all'Hamiltoniana totale H. In questo modo oltre alia possibile dipendenza dal 
tempo del'operatore, l'operatore dipendera dal tempo perche scritto in tale 
rappresentazione (0&(t)) 

Di particolare interesse e la risposta dell'operatore densita ad un poten- 
ziale esterno che possa essere scritto nella forma: 

5V ext (t) = I d 3 * SV^x, t)p CT (x) . (5.3) 

In questo modo la formula della risposta lineare puo essere scritta nel seguente 
modo 

/+oo p 
dt' / d^Q{t-t')8V^t') 
-oo J 

<*(0)|^(x,*),^(x , ,O]|*(0)> , (5.4) 

dove l'integrazione sul tempo e stata estesa a tutto l'asse reale sfruttando il 
fatto che la perturbazione e nulla per t' < to e inserendo la funzione Q(t — t') 
per i tempi t' > t. La funzione di risposta puo essere defmita come (torniamo 
qui ad utilizzare la notazione 1 = "x.\,ti, aj): 

X R (l ! 2) = (M/(0)|[ / 3(l),p(2)]|vl/(0))e(t 1 -t 2 ) . (5.5) 

In questo modo risulta infatti 

5p(l)= X R (l, 2)V ext (2) , (5.6) 

dove stiamo integrando sulle variabili ripetute. Troviamo cosi la defmizione 
classica di funzione risposta come la variazione della densita rispetto ad una 
perturbazione esterna applicata. 

Quella che abbiamo qui definito h in realta la funzione di risposta "ritarda- 
ta" . Fino ad ora avevamo invece lavorato con grandezze "T-ordinate" , poiche 
queste sono le uniche che possiamo sviluppare utilizzando in teorema di Wick 
e quindi i diagrammi di Feynman. La polarizzazione o funzione di risposta 
T-ordinata infatti, se calcolata esplicitamente utilizzando la defmizione (13.611) 
risulta pari a 

X (l,2) = <*(0)|T[ ( yp(l)^(2)]|*(0)> , (5.7) 

dove 5p = p — p° e l'operatore fluttuazione rispetto alia densita media p° = 
(tf(0)|fltf(0)>. 

I due tipi di grandezze possono in realta essere ricollegate (insieme con 
le grandezza "avanzate") utilizzando la rappresentazione di Lehmann delle 
stesse. Vediamo come nel prossimo paragrafo. 
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5.2.1 La rappresentazione di Lehmann della funzione 
risposta 

Riprendiamo la definizione di rappresentazione di Lehmann (I3.30P della fun- 
zione di Green e svolgendo passaggi del tutto analoghi a quelli gia svolti per 
tale grandezza otteniamo la rappresentazione di Lehmann delle due funzioni 
risposta. La sola differenza sara questa volta data dal fatto che il termine 
1 = non compare poiche viene cancellato in mo do esatto sia nella funzione 
T-ordinata che in quella ritardata. Scriviamo dunque i due risultati 

r, , , _v-(*o|n(x)|* / )(* / |n(x')|* ) (* |n(x')|*7)(*/|n(x)|* > 
X l x ) x ) — / 4 



lo — Uj + %r\ uj + oji + ir\ 



(5.1 



, (*o|n(x)|*j)<*j|n(x')|*o> (^ |w(x')|^/)(^/|w(x)|^ > 



U - Ui + 17] UJ + Ul - 17] 



(5.9) 

Vediamo come la funzione ritardata e la seconda funzione T-ordinata differi- 
scano solo per il segno della parte immaginaria a denominatore. 

Lo spettro di eccitazione oui scritto a denominatore della (15. 8p e della 
(15. 9p e costituito dagli autovalori di tutti quegli stati che hanno lo stesso 
numero di particelle del ground-state (poiche l'operatore densita commuta con 
l'operatore numero di particelle) relativo all'energia dello stato fondamentale: 
ujj = Ej — Eq Per poter ottenere dunque tali valori non dobbiamo far altro che 
calcolare i poli di tale grandezza. Vediamo infine come la funzione ritardata e 
la funzione T-ordinata differiscano solo per il segno della parte immaginaria 
a denominatore. 



5.3 L'approccio many-body 

5.3.1 Perche la funzione di Green a particella singola 
non e sufficiente 

L'idea per poter studiare la funzione risposta e quella di scrivere un'equazione 
di tipo Dyson per tale grandezza; nelFambito dell'approccio many-body delle 
due grandezze studieremo la funzione risposta T-ordinata. A questo scopo 
riscriviamo l'equazione di Hedin (I3.63|) per la polarizzazione inserendo al posto 
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del vertice T la sua equazione di Hedin (13.651) nella versione completa 

r(l; 2, 3) = 5(1, 2)5(1, 3) + ±r(l; 4, 5)G(6, 4)G(7, 5) 

[5(2,3K(3) + S^(2,3)] . (5.10) 



5G(6, 7) 
In questo modo otteniamo 

n(l, 2) = G(l, 2)G(2, 1) + G(l, 3)G(1, 4) 

^ 6) [5(2, 3K(3) + S*(3, 4)] G(5, 7)G(6, 8)T(7, 8; 2) . (5.11) 

Se ora cerchiamo di ricavare un'equazione di tipo Dyson per la polarizzazione, 
definendo la grandezza IIo(l, 2) = G(l, 2)G(2, 1) come ordine zero del nostro 
sviluppo, ci accorgiamo immediatamente che non e possibile ottenere il no- 
stro obiettivo poiche le due funzioni di Green che compaiono a sinistra della 
derivata funzionale non sono chiuse, cosi come non lo sono le due funzioni di 
Green che compaiono a sinistra della stessa derivata funzionale. Se cosi fosse 
potremmo infatti riscrivere l'equazione appena ricavata utilizzando la (I3.63P 
e la definizione appena introdotta di Ho : 

n = n + ivai , (5.12) 

con la grandezza K (il kernel dell'equazione) definita da 

K(l, 2; 3, 4) = [5(1, 2)v H (2) + E*(l, 2)] . (5.13) 

II problema e che compare una grandezza sconosciuta, un kernel a quattro 
punti; fatto che in qualche modo e legato alia forma dell'operatore di in- 
terazione, che e un operatore a due particelle e che porta alia comparsa di 
grandezze a quattro punti come avevamo gia visto nell'equazione del moto 
per la funzione di Green (I3.26P in cui compare la funzione di Green a due 
particelle. Nello sviluppare la teoria per il ground-state siamo riusciti a su- 
perare questa difficolta introducendo la self-energia del sistema e quindi le 
equazioni di Hedin. Sebbene anche in questa formulazione compaia lo stesso 
termine a quattro punti nell'equazione per il vertice, termine che costituisce 
la fonte di tutte le difficolta che si incontrano del tentativo di risolvere tale 
set di equazioni in modo esatto, per lo studio del ground state e possibile 
fare approssimazioni (l'approssimazione GW) che ci permettano di lavorare 
semplicemente con grandezze a due punti. 
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Cercando di calcolare lo spettro di eccitazione del sistema non e invece 
piu conveniente continuare a formulare la teoria in termini di grandezze a due 
punti. La ragione fisica di questa necessita puo essere individuata nel fatto che 
quando eccitiamo un sistema dal suo stato fondamentale oltre a promuovere 
un elettrone in uno stato eccitato, creiamo una buca (hole) nel sistema e che 
l'interazione tra i due oggetti diventa decisiva nello studio del sistema. Non e 
piu dunque possibile relegare il ruolo dell'mterazione coulombiana schermata 
particella-buca in secondo ordine, all'interno della self energia. II problema 
in altre parole dipende dal fatto che, poiche non siamo in grado di risolvere la 
teoria in modo esatto dobbiamo cercare di "mettere in evidenza" le grandezze 
che hanno un ruolo dominante dal punto di vista fisico. 

5.3.2 L'equazione di Bethe Salpeter 

L'idea per costruire l'equazione di Bethe-Salpeter e quella di partire dalla 
funzione di Green a due particelle. A questo scopo introduciamo una nuova 
grandezza che definiamo come la differenza tra la funzione di Green a due 
particelle, che descrive appunto la propagazione delle due mentre interagisco- 
no tra loro (oltre che con il resto del sistema), e il prodotto di due funzioni di 
Green a particella singola, che descrivono la propagazione delle due particelle 
senza tener conto della reciproca interazione. 

iL(l,2;3,4) = -G 2 (1,2;3,4) + G(1,3)G(2,4) . (5.14) 

Questa grandezza descrive evidentemente il moto accoppiato delle due par- 
ticelle, possiamo immaginarla come la correzione alia semplice propagazione 
indipendente delle due. I passaggi che dobbiamo or a svolgere per la costruzio- 
ne di un'equazione di tipo Dyson per questa grandezza sono del tutto simili 
a quelli che abbiamo utilizzato per costruire l'equazione di Hedin per la fun- 
zione risposta (o polarizzazione) e quindi l'equazione scritta nel paragrafo 
precedente con la quale non siamo riusciti a raggiungere il nostro scopo sem- 
plicemente perche stavamo lavorando con grandezze a 2 punti. Partiamo con 
l'osservare che la funzione L(l, 2; 3, 4) coincide con la defmizione del termine 
di destra dell'equazione (I3.52p se invece di considerare un potenziale esterno 
a particella singola ne consideriamo uno a due particelle: 

Procediamo quindi come in precedenza utilizzando la relazione (I3.54p 



(5.16) 



5.4 La Time-Dependent DFT 



61 



e definendo, portando avanti l'analogia con quanto fatto in precedenza, una 
funzione di vertice a 4 punti 

Possiamo ora scrivere l'equivalente dell'equazione (I3.65P nel caso che stiamo 
considerando, cioe con un vertice a quattro punti completo 

r 4 (l, 2; 3, 4) = 5(1,3)6(2,4) + 2 ^ . (5.18) 

Utilizziamo infine le regole di derivazione a catena e definiamo (ora si e 
possibile) l'ordine zero della nostra funzione (L (l,2;3,4) = G(l, 4)G(2, 3)) 
insieme con il kernel 

a . (1 , 2;M) ^MM™) , ,, 19) 

ottenendo la seguente equazione: 
L(l,2;3,4) = G(1,4)G(2,3) + 

G(1 ' 6)G(2 ' 5) SG(8j) M4^ (5.20) 

= L (1,2;3,4) + 

L (1,2;5,6)^(5,6;7,8)L(7,8;3,4) . 

Quella appena scritta e l'equazione di Bethe-Salpeter ed e l'equazione utiliz- 
zata nell'approccio many-body per studiare le eccitazioni del sistema. Con- 
frontando infatti le equazioni (15. lip e (I5.20p si vede come la contrazione della 
funzione L(l, 2; 3, 4) dia come risultato la polarizzazione del sistema 

L(l, 2; 1,2) = 11(1,2) (5.21) 

e dunque lo spettro di eccitazione del sistema. 



5.4 La Time-Dependent DFT 

Per poter studiare il problema delle eccitazioni nell'ambito della teoria del 
funzionale densita il primo passo da fare e quello di estendere la validita 
del teorema di Hohenberg-Kohn al caso in cui FHamiltoniana sia dipendente 
dal tempo. Se consideriamo infatti un'Hamiltoniana in cui e presente un 
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potenziale die dipende dal tempo il nostro sistema non presenta piu uno 
stato fondamentale e dunque non e piu possibile utilizzare il teorema di HK. 
Tale estensione e stata elaborata da Runge e Gross nel 1984 e definisce una 
corrispondenza biunivoca (a meno di una costante additiva dipendente dal 
tempo per il potenziale e di una fase dipendente dal tempo per la funzione 
d'onda) tra la densita, la funzione d'onda e il potenziale esterno. A questo 
scopo viene richiesto che sia gia definite uno stato iniziale al tempo t Q e che il 
potenziale esterno sia sviluppabile in serie di Taylor nell'intorno dell'istante 
to stesso. Non siamo pero ora interessati ai dettagli della dimostrazione, per 
la quale, cosi come per gli altri dettagli della teoria del funzionale densita, 
rimandiamo chi fosse interessato alle citazioni della bibliografia. Cio a cui 
siamo ora interessati e mostrare come possa essere utilizzata l'estensione al 
caso dipendente dal tempo della teoria per ottenere lo spettro di eccitazione 
del nostro sistema. In particolare in questa sezione lavoreremo nell'ambito 
della teoria della risposta lineare, per la quale, oltre che nel risultato di Runge 
e Gross, la teoria trova giustificazione formale in un altro lavoro di Gross e 
Kohn. 

5.4.1 Utilizziamo la teoria della risposta lineare nella 
TDDFT 

L'idea per ottenere lo spettro di eccitazione di un sistema a partire dalla TD- 
DFT e semplicemente quella si studiare la relazione tra una perturbazione 
esterna dipendente dal tempo introdotta nel nostro sistema e la variazio- 
ne della densita. Entrambe queste grandezze infatti devono corrispondere a 
quelle reali nell'ambito di questa teoria. 

Consideriamo dunque nuovamente il nostro sistema descritto dall'Hamil- 
toniana 

H = f + W + V ext , (5.22) 

in cui il potenziale esterno non dipende dal tempo e la corrispondente Hamil- 
toniana efficace della DFT 

H DFT = f + V ext + v H [p}+v xc [p] . (5.23) 

Inseriamo quindi una perturbazione 5V ext (t) dipendente dal tempo e ve- 
diamo come varia FHamiltoniana della DFT 

5v,(t) = 5V ext (t) + 5v H + 5v xc , (5.24) 

dove 5vh = vh[p'] — Vjj[p\ e 5v xc = v xc [p'} —v xc [p\. La variazione della densita 
potra essere scritta, all'ordine lineare, come 

6p{l)= Xo S (h2)Sv s (2) , (5.25) 
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dove la funzione di risposta \ KS e stata indicata con l'apice KS in quanto 
descrive la risposta delle particelle del sistema di KS, che ricordiamo sono 
particelle non interagenti. Notiamo come rispetto alia teoria many-body qui 
utilizziamo la funzione di risposta ritardata. 

Poiche siamo interessati alia variazione della densita rispetto al poten- 
ziale esterno trasformiamo l'equazione precedente scrivendo esplicitamente la 
forma del potenziale efficace. In particolare, per consistenza con l'appros- 
simazione in cui stiamo lavorando, linearizziamo i potenziali di scambio e 
correlazione e il potenziale di Hartree rispetto alia variazione della densita. 
In questo modo risulta 

v H [p'} ~ v H [p] + P Sp 
bp 

Vxc [p'}^ Vxc [p} + ^di6p . (5.26) 
dp 

Definiamo ora la grandezza 

f un <?\ - foMQ) Sv xc \p](i) 

U1(1 ' 2) -i#r + i^r ' (5 - 27) 

la sostiamo nell'equazione (I5.25P e, utilizzando la ( 15. 6p . otteniamo 

5p(l)= (x^ s (l, 2)5(2, 3)5(3, 4) 

+ X ^(l,2)/^ c (2,3)x R (3,4))^ t (4) . (5.28) 

Abbiamo cosi scritto al relazione tra il potenziale esterno e la densita uti- 
lizzando gli elementi della formulazione many-body della TDDFT. Poiche in 
particolare sappiamo che i poli della funzione Xo non sono altro che le diffe- 
renze tra gli autovalori di KS possiamo gia vedere come all'ordine zero della 
teoria (scegliendo cioe fnxc = 0) sia in effetti giustificato considerare gli au- 
tovalori di KS per valutare le eccitazioni del sistema. Tale approssimazione 
e d'altra parte spesso insufficiente a dare una descrizione soddisfacente delle 
proprieta fisiche di un sistema qualsiasi. Tutto questo sara poi meglio chiarito 
nella formulazione di Casida della TDDFT (Capitolo 6). 

5.4.2 Un'equazione di Dyson per la TDDFT 

II risultato appena ottenuto pub essere facilmente riscritto sotto forma di 
equazione di Dyson per la funzione di risposta. Un semplice confronto dell'e- 
quazione ( 15. 28p con l'equazione ( 15 .6p ci porta 

X R (l,2) = X KS (l,2)+ X ^(l,3)W3,4) X fi (4,2) . (5.29) 
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Sottolineiamo ancora una volta la differenza tra la teoria many-body e la TD- 
DFT. In quest'ultima infatti siamo riusciti a scrivere un'equazione di Dyson 
per la funzione risposta utilizzando grandezze a due punti poiche a differenza 
della teoria many-body stiamo utilizzando due grandezze diagonali nello spa- 
zio r. II fatto che la funzione risposta sia diagonale in spazio diretto non deve 
pero far pensare che stiamo utilizzando grandezze differenti. Se scriviamo 
infatti la funzione risposta in seconda quantizzazione vediamo che essa e in 
effetti il valor medio di quattro operatori, sia in TDDFT, sia nelle equazioni 
di Hedin, sia nell'equazione di Bethe-Salpeter. La sola differenza e che nello 
sviluppare la teoria nell'ambito del formalismo many-body dobbiamo consi- 
derare anche i termini "non diagonali" mentre per la teoria del funzionale 
densita questo non e necessario. Ancora una volta vedremo come nella for- 
mulazione di Casida, quando proietteremo le grandezze nello spazio degli stati 
anche l'equazione per la TDDFT diventera una grandezza a quattro punti. 

5.5 Approssimazioni 

5.5.1 L'approssimazione adiabatica LDA 

Come per il computo dello stato fondamentale del sistema, anche per il com- 
puto delle eccitazioni la teoria esatta non puo essere risolta e dobbiamo uti- 
lizzare delle approssimazioni. Poiche nell'ambito di questa tesi calcoleremo 
spettri di eccitazioni con la TDDFT indicheremo in questa sezione quella che 
e l'approssimazione utilizzata nell'ambito di tale teoria: l'approssimazione 
adiabatica. 

Tale approssimazione e in effetti l'equivalente dell'approssimazione locale, 
presentata per la DFT, per la variabile temporale. In termini matematici 
consiste nell'utilizzare la seguente relazione 



dove nell'ultimo termine abbiamo indicato la variabile temporale a pedice 
della densita per indicare che non stiamo facendo la derivata funzionale ri- 
spetto a tale variabile ma solo rispetto alia variabile spaziale. L'energia 
di scambio e correlazione E xc dipende cioe solo dalla densita al tempo t e 
nell'approssimazione vengono dunque persi tutti gli "effetti di memoria". 

Questo e particolarmente chiaro nella scrittura del kernel della TDDFT 
nell'ambito di tale approssimazione che risulta 




SA xc [p] 5E xc [p t ] 



(5.30) 




/ xc [p](x 1 t 1 ,x 2 t 2 ) 



6v xc [p](x.i,t 1 ) 
5p(x 2 ,t 2 ) 




d[ptJ( x i) 



(5.31) 



<fal( x 2) 
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Nello scrivere il kernel abbiamo qui trascurato il termine di Hartree, poiche 
infatti il campo di Hartree e statico l'approssimazione adiabatica non avra al- 
cun effetto su tale campo e quindi sul suo kernel che e d'altra parte conosciuto 
in modo esatto 

/ // [p](x 1 ,x 2 ) = 1 f . (5.32) 

|Xl — X2| 

Se mettiamo insieme l'approssimazione adiabatica con l'approssimazione 
LDA utilizzata per l'energia di scambio e correlazione otteniamo la seguente 
espressione per il kernel 

f£ DA (l, 2) = 5(1, 2) (2^^ + P^|^) • (5.33) 

Possiamo infine inserire una parametrizzazione per l'energia di stato omo- 
geneo in funzione della densita scrivendo cosi in modo esplicito il kernel in 
funzione della densita. A questo scopo scriviamo quella che e l'approssi- 
mazione locale dell'energia di scambio, trascurando l'energia di correlazione. 
Riprendendo lo sviluppo fatto nel capitolo 4 e in particolare la (14 .5p otteniamo 

e h x om (pM)) = -rfoP* 1/3 (x) 

f x ALDA (l,2) = -6(1,2) (^f /3 (xi)) , (5.34) 

da cui possiamo osservare l'andamento divergente del kernel nelle regioni in 
cui la densita tende a zero. Inserendo anche il termine di correlazione (14.71) 
e possibile vedere che anche questo diverge in zone in cui la densita tende a 
zero. 



5.5.2 L'approssimazione adiabatica LSDA 

Per il nostro lavoro dobbiamo cosiderare un kernel che dipenda dallo spin 
e dunque anziche utilizzare, nell'ambito dell'approssimazione adiabatica, il 
funzionale LDA utilizzeremo il funzionale LSDA. Avremo dunque un kernel 
definito come 

w^) = w • < 5 - 35 > 

dove le variabili di spin sono come sempre incluse nella notazione 1 = (x l5 1±, a\). 
In approssimazione A-LSDA otteniamo 

/ ] LSDAf f\ 

f* LDA (h 2) = *(ti - h)5(^ - ^ I d€ * c tP ' U 



dp ai 

de L J DA ( P X) , dh L J DA (p,() 



P 7 , , (5-36) 



dp a2 dp ai dp a2 
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dove la derivazione rispetto alia componente a della densita e intesa come 
(utilizzando la regola di derivazione per le funzioni composte) 

di« " dp d( ' { } 

± dipende dalla componente della densita che stiamo utilizzando: + per py, 
- Per Pi- 



Parte II 
Studio di sistemi isolati 



Capitolo 6 



La formulazione di Casida della 
TDDFT 



6.1 Introduzione 

L'idea di Casida e essenzialmente quella di proiettare la formulazione della 
TDDFT nello spazio degli stati o spazio delle configurazioni. In altre parole 
invece di rappresentare le grandezze della teoria in spazio diretto (x, a) o 
in spazio reciproco (k, a) sceglie di rappresentarle rispetto ad una base di 
funzioni d'onda di particella singola. In particolare la base scelta e quella degli 
autostati dell'Hamiltoniana di KS della DFT, in modo che, come vedremo, 
la funzione di risposta libera sia diagonale in tale spazio. La scelta fatta e 
soprattutto mirata alio studio di sistemi isolati, per i quali lo spazio delle 
configurazioni e sicuramente la descrizione migliore poiche la mancanza di 
periodicita del sistema elimina i vantaggi dell'utilizzo dello spazio reciproco. 

In questa nuova rappresentazione l'idea di Casida e in particolare quella 
di studiare come grandezza di riferimento non la funzione di risposta ma la 
polarizzazione del sistema, definita come: 



dove /i e il momento di dipolo elettrico del sistema, mentre E e B sono il 
campo elettrico e il campo magnetico che vengono applicati al sistema. 

6.2 La polarizzabilita del sistema 

Vediamo innanzi tutto come la polarizzabilita e legata alia funzione di risposta 
e quindi alle grandezze fisiche del sistema. A questo scopo consideriamo ora 




J e=b=o 



(6.1) 
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il caso in cui il potenziale applicato sia del tipo 

8V ext (t) = zE z (t) . (6.2) 

Consideriamo quindi il momento di dipolo in direzione x, esso sara definito 
come 

5 fi x = -qSx , (6.3) 

Qui q e la carica elettrica che d'ora in avanti supporremo pari ad uno, men- 
tre 8x sara la variazione del valor medio dell'operatore x = j <i 3 x p(x)x (il 
termine dopo l'uguale e la sua rappresentazione in seconda quant izzazione). 
Ovvero 

fc= / A ^ ' (6 ' 4) 

da cui la componente xz polarizzabilita nello spazio delle frequenze risulta: 

a xz {uj) = - / d x . (6.5) 

j -^ 2 (w) 

La relazione con lo spettro di eccitazione del sistema puo infine essere 
scritta utilizzando la rappresentazione di Lehmann di quest 'ultima insieme 
con la relazione 

(tto|£|*/)<ttj|5|tto) = (*o|2|tti)<*j|£|*o) , (6.6) 

ottenendo: 

(Sj - £ ) 2 - lo 



2(E I -E ){y \x\v I ){y I \z\y ) 



Nelle sezioni seguenti ci occuperemo dunque di riscrivere la TDDFT nella 
nuova base e quindi di calcolare l'espressione della polarizzabilita con tale 
formalismo. Confrontando infine il nostro risultato con l'espressione appena 
scritta potremo ricavare lo spettro di eccitazione e le forze di oscillatore. 



6.3 La TDDFT nello spazio delle configura- 
zioni 

6.3.1 Cambio di base 

Cio che dobbiamo fare e semplicemente scrivere le relazioni tra matrici e 
vettori dello spazio di Hilbert esplicitandone gli elementi in funzione della 
base scelta in modo da avere alia fine delle relazioni tra elementi di matrice, 
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ovvero dei numeri complessi. I dettagli dei passaggi eseguiti sono piuttosto 
noiosi, riportiamo qui semplicemente la regola necessaria per ottenere questo 
risultato. Cio che in pratica si fa e passare dalle formule scritte in spazio (x, a) 
alle formule scritte nello spazio delle configurazioni utilizzando le formule per 
il cambio di base in seconda quantizzazione. La regola da utilizzare e dunque 
la seguente: 



j 

4 = I>}010 , (6.« 



ai 

j 

dove abbiamo scritto un operatore di campo di base % in una nuova base 
j; nel caso in cui la base di arrivo sia continua le sommatorie andranno 
sostituite con degli integrali, ma la regola del cambio di base risultera la 
stessa. Se preferite, poiche stimiamo lavorando con i valori medi e non con gli 
operatori direttamente e possibile utilizzare anche la regola per proiettare una 
grandezza sugli stati, (facendo pero attenzione al fatto che stiamo utilizzando 
grandezze che sono diagonali in spazio x): 

fh,...,jn = J rf3 xi . . . d 3 x n ^(xO^fra) • • • 

/(xx, . . . , x n ) . . . ^* n _ i (x n _ 1 )^ jn (x n ) . (6.9) 
Effettuando il cambio di base per l'operatore densita otteniamo: 

n CT (x,t) = ^l(^ ja (^ a (t)a ja (t) . (6.10) 

Definiamo ora la grandezza 

Pii*(t) = (V KS \al(t)a ja (t)\* KS ) , (6.11) 

che corrisponde ad una densita nello spazio delle configurazioni e le equazioni 
(I5.25P e (I5.28P nel nuovo spazio diventano 

5P ija (h) = J di! xW,hkr{ti - t2)(v s ) hkT (t 2 ) (6.12) 

SPij^h) = J dt' X ija,hkr(tl ~ t 2 )V^(t 2 ) , (6.13) 

dove stiamo sommando sugli indici ripetuti. Scriviamo ora anche il Kernel K 
nello spazio delle configurazioni 

K ija,hkT^ ~ l V - TB — 7TT + Tn — 77T ^ b > 
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e abbiamo tutti gli elementi per riscrivere l'equazione ( I5.28P nella nuova ba- 
se. Prima di scrivere la nuova equazione vogliamo pero anche cambiare la 
base della variabile temporale e passare nello spazio delle frequenze poiche 
in questo modo le relazioni di convoluzione che compaiono diventeranno dei 
semplici prodotti nel nuovo spazio. Applichiamo dunque la trasformata di 
Fourier a tutte le grandezze che compaiono nelle formule e siamo finalmente 
pronti a scrivere il nostro risultato 

SP hkT (u) = 5V^(u) , (6.15) 

dove come sempre sommiamo sugli indici ripetuti. Abbiamo qui direttamen- 
te invertito la funzione risposta in modo da scrivere l'equazione in modo piu 
compatto, in particolare nel farlo abbiamo invertito la funzione x che po- 
trebbe pero essere non invertibile. In realta sappiamo che tale matrice (lo 
vedremo nel prossimo paragrafo) e diagonale nella rappresentazione scelta e 
dunque certamente invertibile. L'unico pericolo risulterebbe costituito dai 
termini diagonali nulli; questi pero possono essere esclusi dalla teoria poiche 
quando x e pari a zero il sistema presenta una variazione di densita nulla 
(cfr. eq (15. 25ft ) e tutte le equazioni scritte non portano alcuna informazione. 



{xfja,hkr( U )) - K ijaMr (u) 



6.3.2 La funzione risposta libera e il kernel 

Prima di procedere studiamo piu in dettaglio la funzione di risposta di un si- 
stema di particelle non interagenti. Questa e infatti conosciuta in modo esatto 
e puo essere facilmente ricavata nello spazio degli stati a partire dalla defini- 
zione di funzione di risposta. Partiamo dalla rappresentazione di Lehmann 
della funzione di risposta (15. 9p scritta pero nella nuova base scelta: 

Xija,hkAU) Z_. U~(Ej- Eq) + IT} 

UJ + (Ej — Eq) + IT) 

II fatto che stiamo considerando la funzione di risposta libera (ovvero quella 
di un sistema non interagente) significa ora che le funzioni d'onda a molti 
corpi sono semplicemente il prodotto di funzioni d'onda a particella singola e 
che quindi gli stati eccitati non sono altro che eccitazioni di singola particel- 
la. Questo si puo vedere chiaramente nella formula poiche, essendo gli stati 
di particella singola ortogonali tra loro, le funzioni d'onda non possono 
differire da \l/ che per una particella passata dallo stato % alio stato j. 
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Come conseguenza di tutto cio la funzione di risposta libera risulta diago- 
nale in questa rappresentazione. Inoltre poiche gli stati eccitati differiscono 
dallo stato fondamentale per un solo autovalore di KS (lo stato i al posto 
di quello j) risultera Ei — E = ef s — ef s . Infine valutando i numeratori 
dell'espressione, se ammettiamo che gli stati possano avere numeri di occu- 
pazione frazionaria, otterremo daH'applicazione dell'operatore di distruzione 
a numeratore esattamente i numeri di occupazione. Possiamo quindi scrivere 
la funzione di risposta libera (la somma sugli stati eccitati scompare perche 
sopravvive solo quello stato eccitato con eccitazione % — > j e per lo stesso 
motivo abbiamo solo uno dei due termini della somma): 

xfja,hkr{ U ) = $<T,T5j,k5i,h- J^R^~KS\ ' ( 6 ' 17 ) 

u - \ e i - e j ) 

Per esprimere il kernel in funzione di grandezze note non possiamo invece 
far altro che scriverlo in funzione del kernel in spazio (x, a) poiche solo in 
questa base siamo in grado di calcolarlo. L'espressione e facilmente scrivi- 
bile utilizzando le regole di cambio base della prima quant izzazione (16. 9 p e 
ricordando che il nostro kernel e diagonale nello spazio e nello spin 

K ijaAkT (u) = J rf 3 xid 3 x 2 ^* (T (xi)^ jCT (xi) J f(' CTiT [p](xi,X2;a;)^ T (x2)^ T (x2) . 

(6.18) 



6.3.3 L'equazione di Casida 

Abbiamo fmalmente tutti gli elementi per scrivere l'equazione che e possiamo 
dire alia base del lavoro di Casida e che non e altro che la formulazione 
dell'equazione della TDDFT nello spazio delle configurazioni: 

(6.19) 

Poiche siamo passati ad una base in cui ne la densita ne il potenzia- 
le esterno sono diagonali avremo a che fare sia con la parte reale che quella 
immaginaria di tali grandezze, poiche come sappiamo le grandezze fisiche ven- 
gono rappresentate da matrici hermitiane i cui termini non diagonali possono 
dunque essere complessi. 

Per mettere in luce la parte immaginaria e quella reale l'idea e quella di 
riordinare la base del nostro spazio nel seguente modo: i < j — > fi a > fj a . Poi 
di dividere l'equazione in due, la prima per la quale consideriamo i termini in 



E 

?„/_ 



fkr f * 



— K, 



hr 



ija,hk 
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cui f ir7 > fj,j e la seconda per cui f ia < fj a . In ciascuna delle due equazioni 
inoltre separiamo in due parti la sommatoria, la prima con f^ T > f kT e la 
seconda con f hr < f kT . Ora invertiamo a coppie il nome di tutti quegli indici 
per cui quello con numero di occupazione minore viene prima in modo da 
ottenere ovunque f ia > f ja e f hr > f kr : 



fhr—fkr>0 



E 

klr 



XXX u ~ ^ hr ~ 6kT > IS (, <\ 

0(T,T<Ji,hOj,k 7 7 J^ijaMrK 1 ^) 

Jkr Jhr 



SPhkrioU 



fhr-fkr>0 

K ija , khT {u)8P khT {u) = 5V^{u) (6.20) 

klr 



fhT-fkr>0 



E 

klr 



0<T, T Oi,hUj,k 



fhr fki 



Kji atk h T {w) 



5P khT (u) 



K Jla . hkT {u)5P hkT {u) = 5V0(co) , (6.21) 

klr 

Possiamo riscrivere le due equazioni in forma matriciale sfruttando il fatto 
che il kernel KHxcfeiti, X2^) e una matrice hermetiana. Otteniamo cosi: 



A(u) B(u) 
B*(u) A*(co) 



C 

-c 



5P(gu) 
8P*(uj) 



SV ext (u) 



dove abbiamo definito le grandezze: 



(6.22) 



C/lT ~~ 6 kT 



Jhr Jki 
Bij a ,hkr{^) = —Kij a ^ k h T (u) 
5cr,T5i,hSj,k 



Kija t hkr(^) 



Cijcr hkr 



fhr fkr 



(6.23) 



Se ora facciamo Fulteriore assunzione che gli orbitali di KS siano reali e 
che il kernel in spazio (x, a) sia una quantita rea otteniamo le relazioni 
A = A* e B = B* e possiamo quindi diagonalizzare la prima matrice da 
sinistra dell'espressione ricavando un'equazione in cui compaiono la parte 



1 I1 che coincide, in accordo a quanto afferma Casida [12) . ad interpretare come infinite 
il tempo di vita delle eccitazioni. 
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reale e quella immaginaria sia del potenziale che della densita: 



A{u) + B(u) 






A(u) - B{uj) 



-c 

-C 



-i%[6P](u) 



-iQ[SV ext ](uj) 



Se ora esplicitiamo l'equazione per la parte reale della densita da quella per 
la parte immaginaria otteniamo 

[(A(u) + B(u)) - u 2 C{A{u) - B{u))-*C] R[6P](u) = 

^[5V ext ] (u) - iuC{A{u) - B(u))- 1 Q[SV ext \ (u) (6.24) 



[(A(u) - B(cu)) - lu 2 C(A(uj) + Bitu))-^] $t[6P](u) = 

Q[SV ext }(tu) - iuC{A{u) + B^Y^V^u) . (6.25) 

Infine nel caso di una perturbazione reale se consideriamo solo la variazione 
della parte reale della densita possiamo riscrivere nel seguente modo la prima 
delle due equazioni: 

$t[5P](oj) = S~ 1/2 {uj 2 - Q{uj)Y l S~ l/2 5V ext {uj) , (6.26) 

dove abbiamo definite le due grandezze 

S(u>) = C{A - B)- l C 

n(co) = S~ 1/2 (A + B)S~ 1/2 . (6.27) 

A partire dall'espressione ( I6.26P la polarizzabilita del sistema puo essere fa- 
cilmente scritta come (scriviamo qui la componente a xz ) utilizzando la (16. 5p 
proiettata nello spazio degli stati: 

a xz (oo) = 2x^S- 1/2 (n(uj) - oo^S-^z , (6.28) 

dove tutte le grandezze vanno ovviamente pensate nello spazio degli stati. II 
simbolo di vettore sopra le grandezze x e z indica che queste sono vettori in 
tale rappresentazione. 

Confrontando questa espressione con la (16. 7p possiamo immediatamente 
osservare che lo spettro di eccitazione del sistema puo essere ottenuto diagona- 
lizzando la matrice Q. L'equazione agli autovalori di cui parlavamo all'inizio 
di questa sezione sara dunque: 

n(w/)F 7 = w|Fj . (6.29) 
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Tale equazione costituisce il risultato principale del lavoro di Casida ed e di 
fatto l'equazione che viene utilizzata per calcolare lo spettro di eccitazioni 
di sistemi finiti nell'ambito della TDDFT. Segnaliamo che il risultato sin 
qui ottenuto e esatto e del tutto generale, in particolare lo stesso formalismo 
potrebbe essere ad esempio applicato partendo, invece che dal sistema di Kohn 
e Sham daH'Hamiltoniana di Hartree-Fock o da una qualsiasi Hamiltoniana 
a particella singola. 

Nell'ambito della TDDFT tale formulazione offre pero l'indubbio vantag- 
gio di essere molto piu semplice. Se infatti consideriamo particelle a funzione 
d'onda reale possiamo facilmente osservare che il kernel risulta simmetrico 
rispetto all'mversione delle coppie di indici o (h, k) e la matrice VL puo 
essere scritta come 

+ 2\j (f~ - f j(T ){e ja - e ia )K ijaihkT (u)^/{f hT - / fcr )(e fcr - e hT ) . (6.30) 

Anche quest'ultimo risultato e esatto. Nessuna ulteriore approssimazione e 
sin qui introdotta. 

Cio che in pratica viene fatto e diagonalizzare questa equazione scegliendo 
una certa approssimazione per l'energia di scambio e correlazione (che rappre- 
senta come sempre il termine approssimato della TDDFT) e di conseguenza 
per il kernel, che ricordiamo e la derivata seconda dell'energia. L'approssima- 
zione che solitamente viene fatta e che utilizzeremo anche noi nelle molecole da 
studiare nell'ambito della tesi, e l'approssimazione adiabatica, gia introdotta 
nel capitolo precedente. Tale approssimazione, nel formalismo introdotto da 
Casida, corrisponde a considerare un kernel statico e quindi una matrice Q 
statica, partendo dalla trasformata di Fourier della (I5.33p . 

6.4 Studio di sistemi a shell chiusa 

I primi studi di Casida sulla TDDFT sono stati applicati in particolare per 
l'analisi di sistemi a shell chiusa. Tali sistemi hanno la particolarita di avere 
la funzione d'onda a spin up identica alia funzione d'onda a spin down (spin 
lungo l'asse z) il che permette di trasformare la matrice Q in una matrice 
a blocchi sfruttando anche le simmetrie risultanti dal considerare l'approssi- 
mazione adiabatica. Per illustrare tali simmetrie scegliamo di utilizzare un 
sistema modello (fig. 16. ip : un sistema con tre stati, che chiameremo i, j e k 
in cui lo stato i e occupato mentre gli stati j e k sono liberi. 
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1 J 



Ik 



Figura 6.1: Descrizione degli stati eccitati per un modello a 3 livelli, 2 
elettroni. 



Scriviamo quindi la matrice f2 per tale sistema: 

/ A n (u) B n (u) a n (uj) \ 

B#(-u) A u {u) P it (u) a u {u) 

a n (-u) M-u) C n (-u) D n (-u) 

\ «u( _w ) Ai.t(- W ) c u( u ) / 

dove gli indici che variano sia sulla righe che sulle colonne sono 1 = ij t, 
2 = ij 4, 3 = ik t, 4 = ik I. Queste sono infatti le uniche eccitazioni possibili 
del sistema. A, B, C, D, a e (3 indicano qui gli elementi della matrice Q e 
possono essere costruiti a partire dalla f )6.30p (e non devono essere confuse con 
la notazione utilizzata in precedenza per le grandezze A, B e C); la scelta di 
nominarle in questo modo e stata fatta per mettere in evidenza le grandezze 
che risultano poi identiche nell'ambito delle approssimazioni eseguite. Ora 
l'approssimazione adiabatica per il kernel della TDDFT fa si che gli elementi 
con argomento u siano identici a a quelli con argomento — u (i.e.: A(u) = 
A(—lu) = A(0)), mentre il fatto che stiamo considerando un sistema a shell 
chiusa fa si che risulti A^ = Ay,, = C44, = e infine = f3^u 



2 Lc identita D^i = c = B^i sono scmpre vere, anche nclla formulazione esatta. 
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In questo modo la matrice risulta presentare gli autovettori nella seguente 
forma: 

/ a \ / c\ I e \ / g \ 

a c —e —g 

b d f h 

\b J \d J \-f J \-h J 

La matrice puo essere di conseguenza ridotta in due blocchi e di fatto possiamo 
pensare di risolvere due sistemi agli autovalori con matrici di dimensione 
dimezzata: 

( A + B a + (3 \ ( A — B a - (3 \ 
\ a + (3 C + D J \ a-/3 C-D ) ' 

Cercare le eccitazioni di un sistema a shell chiusa utilizzando il formali- 
smo di Casida consiste nel cercare di diagonalizzare le due matrici appena 
scritte; troveremo dunque due gruppi di eccitazioni. Nella prossima sezione 
cercheremo dunque di interpretare il significato fisico delle due matrici; a que- 
sto scopo sara necessario passare alio studio delle funzioni d'onda del sistema 
sebbene queste dal punto di vista formale non abbiano alcun significato fisico 
(ricordiamo che sono le funzioni d'onda di Kohn e Sham). 



6.4.1 Interpretazione delle funzioni d'onda degli stati 
eccitati e regole di selezione 

Questa parte come gia detto non trova una giustificazione formale altrettanto 
rigorosa che le sezioni precedenti ma, come afferma lo stesso Casida, risul- 
ta indispensabile per interpretare la teoria stessa. Per farlo partiamo da un 
modello ancora piu semplice di quello considerato in precedenza che e costi- 
tuito da un sistema con due livelli, una libera ed una occupata. In analogia 
a quanto fa Casida possiamo immaginare che questo modello rappresenti la 
molecola dell'Idrogeno. 

Consideriamo dunque la molecola di idrogeno tenendo in considerazione 
solo i primi due stati elettronici % e j. Lo stato % corrisponde all'orbitale sigma 
legante e lo stato j all'orbitale sigma anti-legante. Lo stato fondamentale sara 
il singolo determinante di Slater: 
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gli stati eccitati saranno invece 

$1 = l^t^'tl 

= I^M^I • 

Questi stati possono essere costruiti cercando funzioni che siano autostati 
degli operatori hxs, S 2 e S z . 

Se ora guardiamo al potenziale esterno applicato al nostro sistema vediamo 
che esso e scritto nella forma 

V ext (t)=J2 [ &(x)&(x)^(x,t) 

Possiamo pero riscriverlo come 

'1 



(6.31) 



V ex \t) = / cf x 



2 



[^(x, t) + Vj"*(x, t)] [p t (x) + p,(x)] + 



^[^(x,t) - ^(x,t)][p t (x) - &(x)] 



(6.32) 



e quindi otteniamo: 

y-'(t) = y d 3 x [F(x, t)p(x, i) + B,(x, t)m z (x, t)] . (6.33) 

Da questo possiamo ricavare le regole di selezione del sistema che corri- 
spondono a 

5S = 0, 1 5S Z = (6.34) 

e vediamo quindi che le uniche eccitazioni possibili per il sistema considerato 
sono quelle verso gli stati 

*o *o . (6-35) 
dunque un'eccitazione di singoletto ed un'eccitazione di tripletto. 
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6.4.2 Interpretazione delle matrici di Casida 

Se analizziamo gli stati verso i quali le eccitazioni sono rese possibili dalle 
regole di selezione appena trovate notiamo come essi sono formati dalla somma 
e dalla differenza delle singole eccitazioni tra stati di KS. In altre parole la 
trasformazione che ci farebbe passare dalla base degli stati di KS alia base 
degli stati eccitati e caratterizzata dai vettori 



Se confrontiamo questi con gli autovettori che dividono in due blocchi la 
matrice di Casida capiamo come i due gruppi di eccitazioni che otteniamo sono 
costituiti dal gruppo delle eccitazioni di tripletto e dal gruppo delle eccitazioni 
di singoletto. A questo punto puo essere utile scrivere in modo esplicito come 
sono fatti gli elementi delle due matrici rispetto agli indici di spin. Utilizzando 
nuovamente le simmetrie dei sistemi a shell chiusa otteniamo 



da cui abbiamo evidenziato come lo spin venga trattato nella divisione in due 
blocchi della matrice. In altre parole la matrice viene divisa in due blocchi: 



un blocco per le eccitazioni di singoletto ed un blocco per le eccitazioni di 
tripletto. 

6.5 Studio di sistemi a shell aperta 

La teoria sin qui scritta e di fatto valida ed esatta, (se considerata nella sua 
forma esatta) anche per sistemi a shell-aperta. La matrice Q da diagonaliz- 
zare sara sempre esprimibile nella forma f)6.30p dato che di fatto fino a quel 
punto nessuna ipotesi e stata fatta da questo punto di vista. Poiche d'altra 
parte abbiamo visto come lo studio delle simmetrie rispetto alio spin possano 
semplificare lo studio della matrice per lo studio di sistemi con stato fonda- 
mentale non polarizzato cerchiamo di capire anche nel caso di sistemi con 
stato fondamentale spin polarizzato (polarizzazione collineare) se la matrice 
di Casida puo essere separata a blocchi. Come in precedenza lavoriamo in 
approssimazione adiabatica per il kernel della TDDFT. 

Partiamo nuovamente da un sistema modello con tre bande (fig. 1(3. 2j) 
i, j e k di cui il primo occupato da due elettroni, il secondo solo da un 





Q s — ^tt + 
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Figura 6.2: Descrizione degli stati eccitati per un modello a 3 livelli, 3 
elettroni. 



elettrone ed il terzo vuoto. Nella matrice f2 compariranno gli elementi die 
descrivono le transizioni 1 = jk t, 2 = ij ^, 3 = ik t> 4 = ik \,, abbiamo 
dunque una matrice di dimensione 4x4. In questo caso pero notiamo subito 
che la matrice non presenta piu alcuna simmetria e dunque non puo piu 
essere separata in due blocchi. La simmetria della matrice rispetto alio spin 
sembra dunque scomparsa. D'altra parte sappiamo che anche per sistemi a 
shell aperta valgono le stesse regole di selezione individuate in precedenza, e 
dunque partendo da un sistema come il nostro in cui S = 1/2 e S z = 1/2 
dovremmo avere due possibili blocchi di eccitazioni, quelle che conservano 
lo spin, eccitazioni di doppietto, e quelle per cui AS = 1, eccitazioni di 
quadrupletto. Cerchiamo di capire perche non ritroviamo tali simmetrie nella 
matrice costruita. 

6.5.1 II problema delle eccitazioni multiple 

Partiamo sempre dal sistema modello con tre elettroni. Imponendo come in 
precedenza che lo stato fondamentale e gli stati eccitati siano contempora- 
neamente autostati di h KS S 2 e S z otteniamo le seguenti funzioni d'onda: lo 
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stato fondamentale con un doppietto 

$1/2 = Wit'&U&ftl = 

e i suoi stati eccitati con quattro doppietti 



(f> 2 

^1/2 — 


\iik) 


^1/2 — 




^1/2 — 




(b 2 

^1/2 — 





e uno stato di quandrupletto 

$1/2 = "/jfdtf fc > + \ijk) + \ijk)) . 

Abbiamo qui gia tenuto in considerazione le regole di selezione per lo spin, 
la barra sopra uno stato indica che questo e occupato da un elettrone a spin 
down. Notiamo immediatamente che Faver imposto che le funzioni d'onda 
siano autostati di S 2 fa si che gli stati eccitati di KS debbano mischiarsi con 
uno stato che dal punto di vista delle funzioni d'onda di Kohn e Sham e uno 
stato a eccitazione doppia (\ijk)). Questo poiche in sistemi a shell aperta 
compaiono stati con tre bande semi occupate e come conseguenza l'operatore 
S 2 di fatto mischia tali stati con altri in cui lo spin di due particelle del sistema 
e stato invertito (questo accadeva anche per sistemi a shell chiusa ma in questo 
caso venivano tra loro mischiati solo stati a eccitazione singola, poiche erano 
presenti solo stati eccitati con due bande semi-occupate). Mancando di fatto 
tale stato nella formulazione da noi fatta la matrice non puo presentare la 
simmetria richiesta e dunque essere separata a blocchi. 

Se cercassimo invece di costruire autostati di h KS , S z e S 2 con le sole 
eccitazioni a particella singola otterremmo tre doppietti 



<b 2 

^1/2 — 


\iik) 






*i/a = 





\ijk)) 



e uno stato di tripletto 



$1/2 = -^(fo'fc) + I W) 
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Vediamo dunque come manchi uno dei doppietti, mentre abbiamo un tripletto 
invece di un quandrupletto. Tra le energie di eccitazione che otterremo dia- 
gonalizzando la matrice ci aspettiamo dunque che potranno essere in qualche 
modo considerate solo quelle eccitazioni che conservano S 2 . 

Poiche d'altra parte la teoria sopra formulata e in principio esatta dob- 
biamo immaginare che tale deficienza vada individuata nelFapprossimazione 
eseguita e dunque nelFapprossimazione adiabatica. In effetti mostreremo nel- 
l'ultimo capitolo di questa tesi il modello che Casida propone per correggere 
tale problema e come nell'analizzare il problema lo stesso Casida individui 
proprio nelFapprossimazione adiabatica per il kernel della TDDFT il motivo 
per cui in tale formalismo non siano presenti tutti gli stati ad eccitazione 
doppia. 

La correzione di questo problema, che come vedremo si tradurra nella 
necessita di scartare alcune delle eccitazioni calcolabili diagonalizzando la 
matrice di Casida, esula pero dagli scopi di questa tesi. Per i calcoli che effet- 
tueremo ci accontenteremo dunque di utilizzare la teoria sin qui costruita e di 
diagonalizzare la matrice completa in approssimazione adiabatica. Abbiamo 
in altre parole terminato di illustrare il formalismo che abbiamo implemen- 
tato e che utilizzeremo nei capitoli successivi per lo studio degli stati eccitati 
di molecole con stato fondamentale spin polarizzato. 

L 'ultima sezione di questo capitolo costituisce invece un tentativo di ge- 
neralizzare il formalismo sin qui presentato. Nei calcoli che faremo le genera- 
lizzazioni proposte non verranno considerate. 

6.6 Proposte per estendere la teoria 

In questa sezione infine cercheremo di mostrare come sia possibile estende- 
re la teoria di Casida sin qui descritta a due situazioni che non sono state 
considerate. 

La prima situazione e quella in cui venga considerata una perturbazione 
esterna di questo tipo 

V ext (t) = f d 3 * O x > ^ a (x)Mx) , (6.36) 

in cui cioe compaiano anche termini non diagonali nello spin. 
La seconda situazione e quella in cui 

[H,S z ]^0 , (6.37) 

in cui quindi FHamiltoniana e Foperatore S z non hanno autostati in comune 
e per cui e dunque necessario considerare degli spinori. 
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Per entrambe le estensioni non dovremo far altro die scegliere in modo 
adeguato lo spazio verso cui proiettare la TDDFT in modo da ottenere la 
corretta formulazione. Altra necessita che comparira in entrambe le teorie 
sara quella di lavorare con un kernel non diagonale rispetto alio spin e dunque 
che dovra partire dalla forma completa dell'energia di scambio e correlazione 
come funzionale della matrice densita 

iWy = ^^ • (6-38) 

6.6.1 Includere la possibilita di spin flip 

Considerare una perturbazione in cui siano presenti anche le componenti non 
diagonali rispetto alio spin significa di fatto considerare la possibilita di una 
perturbazione in cui siano presenti oltre ai termini gia individuati per il caso 
diagonale rispetto alio spin anche le componenti lungo le direzioni x e y del 
campo magnetico. La perturbazione puo in altre parole essere scritta come 

V**{t) = J rf 3 x[^(x,t)p(x,t) + B g (x,t)m z (x,t) 

+B y (x,t)m y (x,t) + B x (x,t)m x (x,t)] . (6.39) 

La principale conseguenza di questa scelta e, come gia detto, che ora an- 
che le componenti non diagonali della densita rispetto alio spin avranno una 
variazione differente da zero, questo perche anche nella perturbazione com- 
paiono tali componenti della densita. Scriviamo dunque l'espressione della 
variazione di densita considerando tali componenti: 

5n af) (xt) = J dV J dt'xa/3 )7 <5(xt,xY)V7f (xV) . (6.40) 

Questa in realta e una generalizzazione della (15.291) nel caso in cui si consi- 
derino componenti della densita non diagonali rispetto alio spin. La gene- 
ralizzazione appena fatta ci appare qui del tutto intuitiva e dunque non ne 
riportiamo una dimostrazione. Scriviamo quindi la funzione di risposta ed il 
kernel nello spazio delle configurazioni. 



Xiaj/3,h-ykS = J J rf 3 xrf 3 xV iQ (x)^* /3 (x)Xa/3, 7 5(xt,X / t')^ 7 (x / )^fc5(x / ) 
Kiaj/3,h~/kS = J J d 3 xd 3 xVia(x)^(x0i^ a /3, 7 5(xt,x't / )^ 7 (x')^fe«5(x') 



(6.41) 
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Vediamo dunque die cosi come per la densita anche per la funzione di 
risposta e per il kernel compaiono dei termini non diagonali nello spin. La 
prima differenza evidente che compare in questa situazione e dunque la pre- 
senza di un ulteriore indice di spin in tutte le formule. Prima di procedere 
con il tentativo di scrivere l'equazione di Casida dobbiamo pero fare un'ul- 
tima considerazione sul potenziale scelto. Come sappiamo infatti l'idea di 
Casida parte dal fatto che, riordinando la base scelta possiamo scrivere una 
relazione tra la variazione della parte reale della densita rispetto alia parte 
reale del potenziale esterno in modo da ottenere una matrice che e di dimen- 
sione dimezzata e di tipo hermetiano. La condizione di realta per il potenziale 
ha pero senso in spazio reale e non nello spazio dello spin. Per questo per 
verificare la condizione di realta del potenziale in quest 'ultimo dobbiamo scri- 
vere le trasformazioni dal primo. Per la situazione studiata da Casida queste 
coincidono con la ( 16 .32 j) in cui compaiono solo coefficienti reali e dunque la 
condizione di realta risulta equivalente nelle due formulazioni. Se conside- 
riamo invece un potenziale non diagonale rispetto alio spin le trasformazioni 
risultano essere 

V(x) = l/2(V n + VU) 
B,{x) = l/2{V n -V u ) 
B x (x) = l/2(V n + Vx) 

B v {x) = -i/2(V n - Vlt) • (6-42) 

Compaiono dunque dei coefficienti immaginari e la relazione tra la realta del 
potenziale esterno in spazio reale e in spazio di spin non e piu quella banale. 

In questo caso, in altre parole, cio che possiamo fare e partire dall'equa- 
zione (16.221) e cercare di estrarre le eccitazioni da questa. Una possibility e 
senz'altro costituita dalla scelta di lavorare in approssimazione Tamm-Dancoff 
(TDA), il che consiste nell'ignorare la seconda riga di tale equazione ottenendo 

(A - uC)5P(u) = 5V ext {u) . (6.43) 

Non sara invece piu possibile costruire la matrice Q(u) di Casida. 

Concludiamo il paragrafo segnalando come di fatto l'introduzione di que- 
sto tipo di perturbazione modifichi evidentemente le regole di selezione per 
lo spin. In questo caso avremo: 

AS = 0,1 ;AS Z = 0,±1 . 
6.6.2 Sistemi non collineari 

La teoria di Casida fin qui esposta prevede che l'operatore S z commuti con 
l'Hamiltoniana totale del sistema e che quindi gli autostati dell'Hamiltoniana 
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siano funzioni d'onda con S z fissato. Cerchiamo ora di studiare come potrebbe 
essere riformulata la teoria nel caso in cui questo non accada. La prima 
conseguenza sarebbe quella che le auto-funzioni a particella singola della DFT 
non sarebbero piu degli spinori con una singola componenete differente da 
zero, il che ci permetteva di lavorare con semplici funzioni d'onda, ma spinori 
con entrambe le componenti differenti da zero. Quello che dobbiamo quindi 
fare b molto semplicemente utilizzare le regole del cambio di base gia viste 
per passare dallo spazio r, a alio spazio degli spinori n. Vediamo quindi come 
risulterebbero le equazioni in questo nuovo spazio. Iniziamo con lo scrivere 
la matrice densita 



n a /3 



x,t) = ^0;(xa)0 m (x/3)^(t)L(t) , (6.44) 



dove con la notazione <j) n intendiamo uno spinore a due componenti ip n ^ e 
ipnl- Sottolineiamo come in generale le nuove funzioni ip nc non siano pero un 
s.o.n.c, ma semplicemente la componente dello spinore a spin up e down. Lc 
relazioni di ortonormalita saranno invece date da 

d 3 r</>*(r)0 m (r) = 5 n , m 
I ^ 3 rC(r)^, ff (r) = 6 n , m . (6.45) 

a J 

Definita quindi la nuova grandezza 

Qnmit) = (^Ksfblp^Ks) , (6.46) 

otterremo le due equazioni 

5Qnm{t) = X™ M {t-t')V^{t') (6.47) 

8Qnm{t) = Xnm,lt(t ~ O^f*^) ■ ( 6 - 48 ) 

Proseguendo in modo identico a quanto fatto in precedenza possiamo scrivere 
il kernel nella nuova rappresentazione, esprimere poi in modo esplicito la 
funzione di risposta libera e infine scrivere la nostra equazione per la funzione 
di risposta totale. Non ripetiamo quindi le considerazioni gia fatte ma ci 
limitiamo a scrivere i risultati per il kernel, 

e per la funzione di risposta libera 
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Infine, come in precedenza esprimiamo il kernel nella nuova base rispetto 
alia base rer: 

K nm>lt (t-t') = j rf 3 rrfV0;(r,«)0 m (r,/3)ir QW7 [p](rt,r't')0*(r / , 5)^(r', 7 ) 

a/3, 57 

(6.51) 

In questo modo abbiamo riscritto il formalismo di Casida per un sistema in 
cui l'Hamiltoniana non commuta con lo l'operatore S z 

Rispetto al caso studiato da Casida dovremo ora tenere in considerazione 
un'importanti differenza dovuta al fatto che avremo a che fare ma con degli 
spinori. Poiche infatti gli spinori sono grandezze intrinsecamente complesse 
(per vedere questo e sufficiente provare a scrivere gli autostati delle matrici 
di Pauli) non sara in alcun modo valida la relazione 

C(r, 7 )0 m (r, 5) = n (r, 7 )C(r, 5) (6.52) 

e come conseguenza la matrice del kernel nello spazio degli stati non avra 
piu quelle simmetrie che abbiamo utilizzato per giungere alia formulazione 
della matrice Q ( I6.27p . Come nel caso precedente dunque lo sviluppo della 
teoria dovra partire dall'equazione H6 .22 j) e dunque una possibile soluzione e 
nuovamente quella di lavorare in TDA (Tamm-Dancoff Approximation). 

Un'ultima breve considerazione sulle regole di selezione per il caso non col- 
lineare. Poiche a determinare le regole di selezione e la perturbazione esterna, 
le regole di selezione saranno identiche a quelle sin qui trovate per i sistemi 
collineari. La non collinearita non porta da questo punto di vista niente di 
nuovo se non forse che perde in effetti di senso considerare la conservazione 
della componente z dello spin. 
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Capitolo 7 
Risultati 



7.1 La molecola di BeH 

7.1.1 Calcoli eseguiti 

Come gia detto nell'introduzione, il lavoro di questa tesi e stato focalizzato su 
una modifica dell'implementazione del formalismo di Casida per la TDDFT 
aH'interno del software ABINIT [H] che permettesse di studiare anche sistemi 
a shell aperta e sbilanciati in spin. 

Prima di utilizzare il programma per lo studio di un sistema con stato 
fondamentale a spin polarizzato abbiamo eseguito alcuni test su di un sistema 
che non presenta effetti di polarizzazione rispetto alio spin. A questo scopo 
abbiamo calcolato con il nostro programma [2] lo spettro di eccitazione della 
molecola N2 utilizzando prima l'opzione nsppol = 1, in modo che il software 
utilizzasse l'implementazione gia presente ed inseguito con l'opzione nsppol = 
2 in modo che il software utilizzasse la nuova implement azione includendo in 
modo separato gli stati a spin up e quelli a spin down. 

Tutti i test hanno dato gli stessi risultati per lo spettro di eccitazione del 
sistema. Abbiamo dunque deciso di proseguire il nostro lavoro cercando di 
studiare un sistema che presenti stato fondamentale spin-polarizzato. 

7.1.2 Presentazione della molecola 

La scelta di utilizzare questa molecola per testare il nostro programma e 
motivata dal fatto che questa e la molecola piu semplice che abbiamo trovato 
a presentare un ground-state a spin polarizzato e dal fatto che di questa 
molecola abbiamo a disposizione dati sperimentali [20] e dati teorici ottenuti 
nell'ambito dello stesso formalismo con cui confrontarci [T6l [T7] . 
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Risultati 



spin up spin down 

vacuum 



-0.07 eV 



-2.34 eV 



-4.63 eV 



-8.35 eV 



4o 



3o 



2a 



4a 



In 



3a 



2a 



-0.02 eV 



-1.59 eV 



-2.83 eV 



-7.99 eV 



Figura 7.1: Schema degli autovalori di KS per il BeH calcolati con il 
programma ABINIT (DFT, approssimazione LSDA). Risultati ottenuti per 
acell=50 x 50 x 50 Bohr, ecut=10 Ha, nband=50. Per la classificazione degli 
stati si veda l'appendice A. 



Per lo studio dello spettro di eccitazione siamo partiti dal dato sperimen- 
tale per la distanza d'equilibrio tra i due ioni della referenza [20] e abbiamo 
lanciato il programma ABINIT in due step successivi: un primo step in cui 
vengono calcolati in modo auto-consistente le funzioni d'onda degli stati oc- 
cupati della molecola e un secondo step in cui vengono calcolate anche le 
funzioni d'onda degli stati non occupati, viene costruita la matrice Q e viene 
diagonalizzata la matrice stessa. 

Lo schema per lo stato fondamentale della molecola che abbiamo ottenuto 
e 2(j 2 3cr 1 ; non compare qui lo stato la poiche il nostro programma lavora 
solo con gli elettroni di valenza, inglobando gli elettroni di core in un pseu- 
dopotenziale [H] . II diagramma completo degli autovalori di KS ottenuto dal 
secondo step e rappresentato in figura I7TT1 
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0.0008 



0.001 



Figura 7.2: Effetti di taglia (L=lato dela supercella) sugli autovalori di KS. 
Energia di cut-off, ecut=10 Hartree 



Abbiamo due serie differenti di autovalori per gli elettroni a spin up e quel- 
li a spin down poiche il sistema e spin polarizzato. Dallo schema riportato in 
figura possiamo iniziare ad analizzare innanzi tutto l'energia di ionizzazione 
del sistema, — e HOMO = 4.63 eV, interpretando come tale l'energia dell'ulti- 
mo stato occupato riferita al livello di vuoto. Vediamo che questo valore e di 
molto inferiore al valore sperimentale, —^hqmo = ^-21 eV [20], ma in perfetto 
accordo con il valore riportato da Casida, — ^ HOMO = 4.60 eV [16]. Vedre- 
mo come questo errore, dovuto all'inadeguatezza del potenziale di scambio 
e correlazione LDA, si tradurra nell'impossibilita di individuare alcune delle 
eccitazioni della molecola. 
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7.1.3 Convergenza per gli orbitali di KS 

Come prima cosa studiamo il comportamento degli autovalori di KS la variare 
di due parametri: le dimensioni della supercella e l'energia di cut-off che 
definisce la base di onde piane. Presentiamo innanzi tutto il grafico con la 
convergenza dei risultati rispetto alle dimensioni della supercella (fig. 17.2]) . 
Abbiamo messo in grafico il variare dell'autoenergia di KS rispetto a 1/L 2 
poiche in questo modo siamo in grado di individuare gli stati di vuoto. Uno 
stato di vuoto e infatti caratterizzato da un'energia del tipo 

h 2 k 2 , 

€j = L , 7.1 

dove m e la massa dell'elettrone, j l'indice del livello eccitato e k e il vettore 
d'onda che, all'interno di una cella di dimensioni finite, vale kj = = -p, 
per j vettore a componenti intere, e quindi l'energia di uno stato di vuoto 
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Figura 7.4: Grafico di convergenza pr le auto-energie di KS rispetto al valore 
del cut-off sull 'energia 



risulta 



4vr 2 ^ 2 ?' 2 



da cui vediamo che gli stati di vuoto possono essere individuati come quegli 
stati che hanno una dipendenza di tipo lineare dal quadrato dell'inverso del 
lato della supercella. Nel grafico (fig. I7.2p possiamo cosi individuare come 
stati di vuoto tutti quegli stati che abbiamo rappresentato in blu. Questo 
diagramma giustifica la rappresentazione fatta dello schema dei livelli (fig. 
17. ip . e in particolare permette di individuare il livello di vuoto come il valore 
a cui tendono tali stati (fig. I7.3P e che verra preso come di riferimento della 
scala di energie. 

Individuati gli stati di vuoto dal grafico (fig. 17. 2p possiamo invice osservare 
come gli altri stati, che chiameremo stati di valenza per distinguerli dagli stati 
di vuoto, convergono rapidamente rispetto alle dimensioni della supercella. 
Alio stesso modo possiamo osservare nel grafico (fig. I7.4p come la convergenza 
sia stata raggiunta anche per l'energia di cut-off del sistema. 
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Tipo 


eccitaz. 




TDLDA 


TDLDA 


Presente 


Calcoli 


Dati 


eccit. 


principale 




Gaussian03 


DeMon2k 


Lavoro 


Quant. 


sper. 




Em 




nu 




TDDFT 


Chem. 

m\ 


m 


n 


3(7- > lvr 


t 


2.3657 


2.2479 


2.40 


2.56 


2.48 


s 


3(7- > 4(7 


t 


4.6633 


4.5103 


4.52 


5.51 




s 


3(7— > 5(7 


t 


4.7576 


4.6300 


4.69 


5.61 




n 


3(7- > 2vr 


t 


4.8451 


4.7047 


4.80 


6.31 


6.31 


s 


2(7- > 3(7 


1 


5.1318 


4.8049 


**** 


6.12 




n 


2(7- > lvr 


n 


5.6338 


5.1685 


5.65 


6.74 


-6.7 


s 


3(7— > 6(7 


t 


5.6414 


5.4803 


**** 


6.71 


-6.7 


A 












6.77 


-6.7 


S 














6.71 


s 














7.02 


n 












7.27 


7.28 



Tabella 7.1: Confronto dei dati riportati con alcuni valori teorici di riferimento 
e dati sperimentali 



7.2 Analisi dei dati 

7.2.1 Confronto con i risultati ottenuti da Casida 

Svolti i test di convergenza per lo studio delle autoenergie di KS possiamo 
passare alio studio delle energie di eccitazione ottenute diagonalizzando la 
matrice di Casida (I6.30p e confrontare i risultati ottenuti con quelli a nostra 
disposizione. Negli articoli [TBI ELY] sono riportate almeno le prime sette ener- 
gie di eccitazione e in particolare nel lavoro [T7] possiamo trovare la transizione 
principale di KS associata ad ogni eccitazione (tab. 17. ip 

Facendo riferimento a questi dati possiamo scegliere, tra le eccitazioni 
che troviamo nel nostro file di output le stesse riportate da Casida, ovvero 
quelle a cui e associata la stessa transizione di KS riportata nella referenza, 
ottenendo cosi i risultati riportati nella tabella ( I7.ip . Vediamo che l'accordo 
con i dati riportati da Casida e ottimo, unica eccezione e costituita dalle 
eccitazioni che coinvolgono la transizione 2a— > 3cr J, che nei nostri calcoli, 
per i parametri indicati in tabella (tab. I7.ip risulta di difficile individuazione 
e che per il momento non riportiamo. L'accordo non e invece altrettanto 
buono se confrontato con i valori teorici di riferimento, tratti da un lavoro di 
LD. Petsalakis [19] (valori teorici presi come riferimento da Casida nei lavori 
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0.0008 



0.001 



Figura 7.5: Convergenza delle energie di eccitazione TDDFT rispetto alle 
dimensioni della supercella 



[H1[T7]) e con alcuni dati sperimentali [2"U] . 

Analizzando le eccitazioni ottenute confrontando i dati con quelli riportati 
da Casida abbiamo pero trascurato di considerare due fatti: 

1. L'approssimazione LDA e responsabile della sottostima dell'energia di 
legame dei livelli molecolari occupati rispetto al livello di vuoto. 

2. L'approssimazione adiabatica, come discusso nel capitolo in cui abbiamo 
presentato la teoria di Casida, dovrebbe impedirci di ottenere, in linea di 
principio, alcune delle eccitazioni del sistema e tra le eccitazioni trovate 
saremo costretti a scartare quelle per le quali lo spin totale del sistema 
non viene conservato. 

La prima di queste due considerazioni ci porta a studiare nuovamente il 
problema degli stati di vuoto gia presentato per gli autovalori di KS e vedere 
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Figura 7.6: Grafici di convergenza delle energie di eccitazione verso stati di 
vuoto rispetto alle dimensioni della supercella. 



come questo si riscontri nello studio delle energie di eccitazione in TDDFT. 
Guardando la figura 17751 possiamo osservare che le eccitazioni che coinvolgono 
stati con energia superiore alio stato 4<r presentano un andamento fortemente 
sensibile alle dimensioni della supercella e possono dunque essere interpretate 
come eccitazioni verso stati di vuoto. 

Come in precedenza possiamo studiare l'andamento delle energie di ec- 
citazione rispetto alia grandezza 1/L 2 per mettere in evidenza gli stati di 
vuoto. Nello spettro compaiono delle serie differenti di energie di eccitazione 
a seconda di quale sia lo stato di partenza da cui l'elettrone viene eccitato 
verso stati di valenza e di vuoto (fig. I7.6|7.7p . 

Decidiamo quindi di scartare tutte le eccitazioni fortemente sensibili alle 
condizioni al contorno e, per lo stesso motivo, di non considerare le eccitazio- 
ni che coinvolgono lo stato 4cr, poiche anche quest'ultimo, pur presentando 
un'autoenergia negativa, ha un comportamento molto piu simile a uno stato 
di vuoto (17. 3p che non ad uno stato di valenza0. La scelta di scartare le ec- 



1 In proposito a questa considcrazionc possiamo anche osservare la serie di immagini, 
che rapprcscntano gli orbitali di KS in spazio reale, ramgurate nella parte finale di questo 
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Figura 7.7: Grafici di convergenza delle energie di eccitazione verso stati di 
vuoto rispetto alle dimensioni della supercella. 



citazioni verso stati di vuoto viene in particolare fatta in accordo con i dati 
presentati nel lavoro [18] per lo studio della molecola A^. In tale lavoro infat- 
ti viene mostrato come l'approssimazione LDA sia responsabile, del 
non corretto andamento asintotico del potenziale di scambio e correlazione, 
di uno spostamento verso l'alto e dunque verso il vuoto di tutti i livelli di va- 
lenza. In questo modo si ha la comparsa di una serie di eccitazioni verso stati 
di vuoto che scompaiono studiando il sistema con un potenziale che presenti 
un andamento asintotico corretto e che quindi l'autrice individua come non 
fisiche. 

La seconda delle considerazioni ci porta a guardare in modo piu dettaglia- 
to le eccitazioni non sensibili alle condizioni al contorno e a cercare un metodo 
per individuare quali di queste siano da scartare. A questo scopo scriviamo 
la tabella con le nostre eccitazioni (j7.2lffl Come gia detto scartiamo tutte le 
eccitazioni verso lo stato 4a. Guardando alle altre vediamo che quelle che in- 

capitolo da cui e possibile osservare nuovamente come lo stato 4a abbia un comportamento 
diffcrcntc dagli altri orbitali di valenza e molto piu simile ad uno stato di vuoto. 

2 Abbiamo scelto di lavorare con una supercella di dimensioni maggiori poiche in questo 
modo siamo anche in gradi di ottenere l'eccitazione 2cr 1 3cr 2 che con questi parametri e piu 
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Risultati 



Eccitazione 


energia 


eccitazione 


Configurazione 


Contributo 




TDDFT 


principale 


elettronica 


dell'eccitazione 








spaziale 


principale 


n 


2.40 


3a ->■ Itt t 


2a 2 l7r 1 


99% 


s 


5.54 


2a ->■ 3a | 


2a x 3a 2 


96% 


n 


5.66 


2a -»■ Itt t4- 


2a 1 3a 1 lvr 1 


99% 


n 


7.14 


2a -»■ Itt |I 


2a 1 3a 1 l7r 1 


97% 



Tabella 7.2: Schema completo delle eccitazioni a convergenza, le eccitazioni 
sono state scelte osservando i grafici (fig. 17.51 17.6[ 17.7ft I parametri dei risultati 
riportati sono ecut 10 Hartree, nband 50, acell 70 x 70 x 70 Bohr 



contrerebbero il problema descritto nel capitolo precedente sono le eccitazioni 
in cui abbiamo 3 stati parzialmente occupati. Osserviamo come per questo 
tipo di eccitazione ne compaiono sempre due che presentano lo stesso tipo di 
occupazione degli orbitali spaziali. Di queste due dunque, immaginiamo che 
una sara l'eccitazione verso uno stato con spin totale maggiore e l'altra l'ec- 
citazione verso uno stato con spin totale S = 1/2 come lo stato fondamentale 
del sistema utilizzando il fatto che AS = 0, 1 e rifacendoci a quanto detto nel 
capitolo precedente. 

Ci affidiamo ora alia regola di Hund ed individuiamo l'eccitazione a spin 
maggiore, cioe quella che vogliamo scartare, come quella ad energia minore e 
dunque tra le due eccitazioni II = 5.657 e II = 7.145 scegliamo di assegnare 
alia seconda S = 1/2. 

Osserviamo che, riguardo ai due problemi segnalati, Casida sembra igno- 
rare il primo, probabilmente poiche lavorando con una base localizzata non 
e in gradi di distinguere stati di vuoto da stati di valenza^. Per il secondo 
invece elabora un modo molto piu raffinato di quello da me proposto in questo 
lavoro (facendo riferimento alia regola di Hund) e chiamato "contaminazione 
di spin" [IT] . Guardando proprio alio studio di tale parametro in effetti, tra 
le eccitazioni da lui scelte, individua come maggiormante "contaminata" la 
stessa che noi abbiamo scartato supponendo fosse a spin maggiore. 



facilmente individuabile. 

3 Nella referenza [TB] Casida individua in realta il problema e cerca di correggerlo stu- 
diando le energie di eccitazione ottenute grazie ad un potenziale con andamento asintotico 
corretto. Tale studio non viene poi piu menzionato nella referenza sucessiva [17] . 
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iiuinero di baude 

(c) TDDFT: convergenza risetto al numero 
di stati 



(d) Particolare del grafico (c) 




Figura 7.8: Test di convergenza per le energie di eccitazione rispetto all'ener- 
gia di cut-off per la base di onde piane e rispetto al numero di stati considerati 
per la costruzione della matrice Q 



7.2.2 Altri test di convergenza 

Abbiamo riportato i test di convergenza fatti per il computo delle autoenergie 
di KS e quindi lo studio delle energie di eccitazione rispetto alle dimensioni 
della supercella. Anche per l'analisi delle eccitazioni del sistema sosno stati 
fatti dei test rispetto all'energia di cut-off e rispetto ad un nuovo parametro, 
il numero di stati considerati per la costruzione della matrice Q, che nel 
computo dello stato fondamentale del sistema non compariva. Dai grafici (fig. 
I7.8p possiamo osservare come i risultati siano a convergenza per entrambi i 
parametri. La dimensione della supercella resta, dal nostro punto di vista, 
il paramentro piu importante per gli studi di convergenza delle energie di 
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Risultati 



Eccitazione 


Dati 


Dati 


Presente 


Calcoli 




sperim. [T6] 


sperim. [20] 


lavoro 


quantum 
chem. [19| 


n 


2.484 


2.48 


2.40 


2.56 








5.54 


5.51 


£ 








5.61 


£ 








6.12 


n 


6.317 


6.31 




6.31 


s 




-6.7 




6.71 


n 




-6.7 




6.74 


A 




-6.7 




6.77 






6.71 










7.02 






n 


7.460 


7.28 


7.14 


7.27 



Tabella 7.3: Raffronto tra le energie di eccitazione (in eV) ottenute nel 
presente lavoro con dati sperimentali e teorici presenti in letteratura. 



eccitazione nell'ambito della TDDFT, soprattutto per l'individuazione degli 
stati di vuoto. 

7.3 Considerazioni finali 

Riportiamo ora una tabella in cui confrontiamo le eccitazioni considerate 
valide con i dati di riferimento nella tabella (17. 3p 

Segnaliamo che lo scegliere per un valore ottenuto una determinata ec- 
citazione sia stato fatto in modo arbitrario, non avendo questa volta alcuna 
indicazione sul tipo di eccitazione, salvo la distinzione £, II. Sembra comun- 
que promettente l'ottimo accordo che in questo modo riusciamo ad ottenere 
con i valori di riferimento con una differenza media, per le 3 eccitazioni calco- 
late, di circa 0.1-0.2 eV. Uno schema riassuntivo dei risultati ottenuti e infine 
riprodotto nella figura (fig 17. 9ft 

Concludiamo l'analisi dei dati sottolineando come nel lavoro svolto le due 
approssimazioni eseguite, l'approssimazione LDA (e quindi l'errato andamen- 
to del potenziale di scambio e correlazione nella regione asintotica) per il com- 
puto dello stato fondamentale e l'approssimazione adiabatica per il computo 
del kernel hanno portato come conseguenza il fatto che non siamo in grado 
di individuare la maggior parte delle eccitazioni del sistema. Sui tentativi di 
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LSDA TDLSDA Quantum chem. Val. Sperimentali 



_ 4 ,| 7.28 eV rr 

7.1. 1 c\* II .- " II 

. 7.02 eV £ 

6 77eV 1A 6.71 eV y 

6 74 eV - ; , z - 

= 3n I z, n, a 

6.40 eV 6.71 eV 5£ -6.7 eV 

... , 6.31 eV 211 II 

6.01 eV 4£ 6.31 eV 



6.12 eV 



5.65 eV II 



5.16 eV 



5.54 eV £ 5.51 eV 2E 



2.29 eV 



2.40 eV II 



2.56 eV 



in 



2.48 eV 



n 



Figura 7.9: Raffronto delle energie di eccitazione calcolate come differenze tra 
autoenergie di KS (ordine zero della nostra teoria), come energie di eccitazione 
in TDDFT (approssimazione LSDA/TDLSDA) e valori teorici e sperimentali 
presenti in letteratura [TBI 20J. L'eccitazione che rappresentiamo in modo 
tratteggiato e quella che abbiamo individuato come da scartare. 



migliorare Fapprossimazione LDA abbiamo visto come una soluzione viene 
presentata in particolare nella referenza [TB] o come un approccio possibile 
sia quello di utilizzare l'equazione di Sham-Schluter presentata all'mterno di 
questa tesi (I4.24p . Tale problema riteniamo sia stato gia studiato in molte al- 
tre occasioni nel computo dello stato fondamentale di differenti sistemi e ben 
delineato. Dal nostro punto di vista crediamo possa essere un obiettivo di non 
sostanziale difficolta implementare una delle possibili soluzioni all'interno del 
software utilizzato [T4"j . 

II problema del miglioramento dell'approssimazione adiabatica invece e ad 
oggi uno degli argomenti di cui si occupa attivamente la ricerca. Abbiamo 
dunque scelto di mostrare, nel prossimo capitolo, un possibile approccio che 
permetterebbe di giungere alia soluzione di tale difficolta e che viene illustrato 
nelle referenze [H |£j. Dal nostro punto di vista Fimplementazione di questa 
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Risultati 



soluzione, o di altre possibili forme dipendenti dalla frequenza del kernel della 
TDDFT, resta una sfida di maggiore difficolta. 
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(a) spin:f, isosuperficie a 10 Bohr 3 ' 2 (b) spin:|, isosuperficie 10 Bohr 3 / 2 




(c) spin:t, isosuperficie 1 Bohr 3 / 2 (d) spin:.],, isosuperficie 1 Bohr 3 / 2 




(e) spin:f, isosuperficie 0.1 Bohr 3 / 2 (f) spin:J,, isosuperficie 0.1 Bohr 3 / 2 



Figura 7.10: Stato 2a; Questa serie di figure, cosi come quelle delle pagie- 
ne successive, mostrano different! isosuperfici della parte reale delle funzioni 
d'onda (la parte immaginaria e nulla). E possibile notare la differenza tra 
gli stati che abbiamo interpretato come di valenza (2a, 3a e 4ir) rispetto agli 
stati di vuoto. In particolare le funzioni d'onda di valenza sono localizzate 
nei pressi degli atomi, mentre gli stati di vuoto risentono delle condizioni al 
contorno. Tutte le immagini sono realizzate con una supercella cubica di 30 
Bohr di lato 
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Risultati 




(a) spin:!, isosuperficie —5 Bohr 3 / 2 (b) spin:|, isosuperficie —5 Bohr 3//2 




(c) spin:f, isosuperficie —0.2 Bohr 3 / 2 (d) spin:|, isosuperficie —0.2 Bohr 3 / 2 




(e) spin:f, isosuperficie 0.2 Bohr 3 / 2 (f) spin:4,, isosuperficie 0.2 Bohr 3 / 2 




(g) spin:f, isosuperficie 5 Bohr 3 / 2 (h) spin:|, isosuperficie 5 Bohr 3 / 2 

Figura 7.11: Stato 3a 
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(a) spin:f, isosuperficie —5 Bohr 3 / 2 (b) spin:^, isosuperficie —5 Bohr 3//2 




(c) spin:']", isosuperficie —0.5 Bohr 3 ' 2 (d) spinij,, isosuperficie —0.5 Bohr 3 ' 2 




(e) spin:f, isosuperficie 0.5 Bohr 3 / 2 (f) spin:|, isosuperficie 0.5 Bohr 3 / 2 




(g) spin:f, isosuperficie 5 Bohr 3 / 2 (h) spin:|, isosuperficie 5 Bohr 3//2 

Figura 7.12: Stato In 
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(g) isosuperficie 1.6 Bohr 3 / 2 (h) isosuperficie 1.8 Bohr 3 / 2 

Figura 7.13: Stato 4a "f, questo stato pur presentando autoenergia negativa 
appare sensibile alle condizioni al contorno con un comportamento simile ad 



uno stato di vuoto. 
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(a) isosuperficie —1 Bohr 3 ' 2 (b) isosuperficie 0.7 Bohr 3 / 2 




(g) isosuperficie 1.3 Bohr 3 / 2 (h) isosuperficie 1.5 Bohr 3 / 2 

Figura 7.14: Stato Aa I, vedi il commento per la serie di immagini dello stato 

4a f. 
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Risultati 




(a) isosuperficie —3 Bohr 3 ' 2 (b) isosuperficie —1 Bohr 3 ' 2 




(c) isosuperficie —0.6 Bohr 3 / 2 (d) isosuperficie 0.6 Bohr 3 / 2 




(c) isosuperficie 1.2 Bohr 3 / 2 (f) isosuperficie 1.4 Bohr 3 / 2 




(g) isosuperficie 2 Bohr 3 / 2 
Figura 7.15: Stato 5cr t> esempio del comportamento di uno sato di vuoto. 
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7.4 II kernel della TDDFT 

7.4.1 Introduzione 

Nella presente sezione presenteremo quello die e un possibile approccio per 
superare le difficolta da noi incontrate nell'ambito dell'approssimazione adia- 
batica per il kernel della TDDFT. Tale approccio e esposto nelle referenze 
P || e nasce dal tentativo di scrivere un kernel a due punti per la teoria 
MBPT. Esporremo qui solo la formulazione esatta poiche siamo semplice- 
mente interessati a mostrare come sia possibile superare le difficolta da noi 
incontrate e non a svolgere calcoli espliciti. 

7.4.2 Un kernel a due punti per la teoria molti corpi 

Vogliamo cercare di formulare un'equazione che leghi il kernel della TDDFT 
con il kernel della MBPT in modo analogo a quanto fatto per il potenziale 
di scambio e correlazione e la self energia, eq. (I4.24p . per la teoria dello 
stato fondamentale. A differenza di quanto fatto in quest'ultimo caso, in 
cui si considera la self-energia esatta del sistema per ottenere il potenziale, 
la grandezza many-body che verra presa in considerazione non sara il kernel 
dell'equazione di Bethe-Salpeter ma una sua versione a due punti. 

Per ottenere tale kernel ripartiamo da quanto mostrato nel capitolo 5 
nella costruzione dell'equazione di Bethe-Salpeter. In particolare ripartia- 
mo dell'equazione (15. lip , questa volta pero lavoriamo con la polarizzazione 
irriducibile 

II*(1,2)=G(1,2)G(2,1)+ 

G(l, 3)G(1, 4) G(5, 7)G(6, 8)f (7, 8; 2) . (7.3) 

Come abbiamo gia visto tentare di costruire un'equazione di Dyson a 
partire da questa si incontrano difficolta legate con il fatto che il kernel e 
una grandezza a quattro punti. D'altra parte nel capitolo sulle eccitazio- 
ni abbiamo rinunciato al tentativo di continuare a lavorare con grandezze 
a due punti poiche abbiamo individuato l'importanza fisica dell'interazione 
particella-buca in tale ambito. 

Cercando invece di lavorare comunque con grandezze a due punti possiamo 
definire l'ordine zero della polarizzabilita come Ilo(l,2) = G(1,2)G(2, 1) e 
riscrivere l'equazione moltiplicando l'ultimo termine per n (l, l')!!^ 1 (3, 1') = 
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5(1, 3) ottenendo: 
n*(l, 2) = n (l, 2) + n (l, l')U^(3, l')G(3, 4)G(3', 5) 

^|^G(6,8)G(7,9)f(8,9;2) . (7.4) 

Come in precedenza consideriamo l'equazione di Hedin per il vertice (I3.65P 
che pero riscriviamo utilizzando le regole di derivazione a catena 

f(l;2,3)= 5(1, 2)5(1, 3) + ^77TT 

5V{) (7 5) 

= 5(1,2)5(1,3)+ w . 

Utilizzando infine la definizione di funzione risposta possiamo scrivere un'e- 
quazione di tipo Dyson per la polarizzazione ridotta a due punti 

n* (i, 2) = n (i, 2) + n (i, 3)fJJ(3, 4)n*(4, 2) , (7.6) 

dove abbiamo definito il kernel dell'equazione usando la notazione della refe- 
renza [9] 

/^(l,2) = ^Hl,3)G(3,4)G(5,3)^M . (7.7) 
7.4.3 Un kernel per la TDDFT 

II kernel f^J viene costruito nella MBPT per cercare di ottenere un'appros- 
simazione successiva a quella GW che possa essere inclusa come grandezza a 
due punti [H |9] . Un'altra importante applicazione risulta pero dal fatto che il 
kernel ottenuto puo essere utilizzato per ottenere il kernel f xc della TDDFT 
in modo esatto, almeno in linea di principio [8]. 

Per vedere questo partiamo nuovamente dall'equazione di Dyson per la 
funzione di Green f)3.4ip scrivendo la self energia completa come H* H = V# +£* 
e quindi il potenziale di Hartree come Vh = vks — v x 



J xc- 



G(l, 2) = G (l, 2) + G (l, 3) (5(3, A)(v KS (A) - v xc (A)) + S*(3, 4)) G(4, 2) , 

(7.8) 

scriviamo quindi la corrispondente equazione per la funzione G 1 

G~\l, 2) = G \l, 2) - 5(1, 2)(v KS (2) - v xc (2)) - S*(l, 2) , (7.9) 

e deriviamo rispetto alia, usando il fatto che la variazione della funzione Gq 1 
e nulla: 

SG-Hh2) „, n ,, Sv KS (2) fa„(2) <S'(1.2) 

^3r=-* (1 - 2)( ^(3r-w ) -^(5r ■ (7 - 10) 
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Moltiplichiamo quindi l'espressione risultante per z'G(4, 1)G(2, 4 + ) e usiamo 
la relazione (I3.54p e le definizioni di Il , xks e fxU ottenendo: 

^(3 ! 4)=n (4,2)(^ 1 5 (2,3)-/, c (2,3))-^(4,l)G(2,4 + ) ^ J . (7.11) 

Infine raccogliamo la funzione n e usiamo la definizione introdotta di f^J 
ottenendo cosi una relazione esatta tra il kernel della TDDFT e f%£- 

/«(1,2)= (x^ 1 s (l,2)-no 1 (l,2))+/^(l,2) . (7.12) 

Poiche siamo interessati alia correzione all'approssimazione adiabatica 
dobbiamo ora applicare la trasformata di Fourier all'equazione scritta. Omet- 
tendo per semplicita di notazione gli indice delle variabili spaziali otteniamo 

UM = (xksM - n ~ V)) + f e JJ{") ■ (7.13) 

Quella appena scritta e un'equazione esatta per il kernel della TDDFT che 
puo essere utilizzata per superare l'approssimazione adiabatica utilizzata nel- 
l'ambito di questa tesi. 
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II computo dello spettro di eccitazione e degli stati eccitati di un sistema 
fisico tramite l'utilizzo di software numerico richiede una quantita di risorse 
molto maggiori di quelle necessarie per calcoli di stato fondamentale [7] . Nella 
presente tesi e stato utilizzato un approccio che da questo punto di vista si e 
rivelato piuttosto agile, perlomeno per i test eseguiti, che sono stati lanciati 
su di un singolo calcolatore e che nel peggiore dei casi hanno richiesto un 
tempo di calcolo di qualche ora. I risultati crediamo siano incoraggianti, ma 
la teoria necessita ancora ulteriori sviluppi. 

L'analisi del formalismo di Casida mette in luce come il computo dell'e- 
nergia di stato fondamentale e degli autovalori di Kohn e Sham possa essere 
legato alle energie di eccitazione del sistema preso in consider azione. L'uti- 
lizzo di tale formalismo evidenzia come lo spettro di eccitazione del sistema 
sia ottenuto a partire dalla differenza tra le autoenergie di Kohn e Sham tra- 
mite due correzioni: la prima dovuta al kernel della TDDFT e quindi agli 
effetti di campo medio e di scambio e correlazione, la seconda dovuta al mix 
tra le transizioni di Kohn e Sham necessario per rispettare le simmetrie di 
spin del sistema. Questo secondo punto, in particolare, delinea la differen- 
za principale tra sistemi a spin equilibrato e sistemi con stato fondamentale 
spin-polarizzato. 

Abbiamo messo in evidenza quali siano alio stato attuale dell'arte i limiti 
di questo tipo di approccio legati alle approssimazioni eseguite. Innanzi tutto 
l'approssimazione LDA per il computo dello stato fondamentale del sistema 
porta ad un errore nella stima dell'energia degli stati occupato rispetto al li- 
vello di vuoto e come conseguenza all'impossibilita di individuare alcune delle 
eccitazioni del sistema. II problema era gia stato riscontrato in letteratura 
[THJ EE] e in particolare e dovuto al non corretto andamento del potenziale 
di scambio e correlazione nella regione asintotica. Una possibile soluzione e 
quella di utilizzare un potenziale di scambio e correlazione ricavato a partire 
dalla MBPT, soluzione presentata nel Capitolo 4, equazione ( 14. 24ft ; sempre in 
letteratura possiamo trovare calcoli TDDFT eseguiti con un potenziale cor- 
retto [18]. Rispetto ai calcoli trovati in letteratura, in cui il computo degli 
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stati viene eseguito con software che sviluppano le funzioni d'onda su una 
base di orbitali localizzate, il lavoro della presente tesi e pero stato realizzato 
lavorando con un codice a onde piane, il che ha permesso di individuare in 
modo piu chiaro, dal punto di vista dell'autore, la presenza di stati di vuoto. 

II secondo problema e invece legato all'approssimazione adiabatica per il 
kernel di scambio e correlazione per il compute degli stati eccitati. Questo 
tipo di approssimazione infatti distrugge la simmetria di spin nel caso di un 
sistema con stato fondamentale spin-polarizzato e porta come conseguenza 
alia mancanza di alcune eccitazioni e alia comparsa di eccitazioni che non 
presentano la corretta simmetria di spin. Nell'ambito di questa tesi abbiamo 
proposto, per lo spettro calcolato, un criterio che ci permetta di individuare 
quali siano le energie di eccitazione soggette a tale problema e nella parte fina- 
le, equazione ( 17.1 5J , abbiamo illustrate una possibile soluzione per superare 
il problema incontrato partendo dal kernel dell'equazione di Bethe-Salpeter. 
Rispetto all'approssimazione LDA, la correzione dell' approssimazione adia- 
batica e uno degli argomenti di cui si occupa attivamente la ricerca e, a 
nostra conoscenza, non sono ancora stati eseguiti calcoli teorici su sistemi 
isolati in cui i problemi connessi all'applicazione di tale approssimazione ai 
sistemi spin-polarizzati siano stati superati con successo. In letteratura viene 
comunque proposto un metodo per individuare le eccitazioni che non pre- 
sentano la corretta simmetria di spin [17]. Sempre in letteratura e possibile 
trovare punti di partenza alternativi a quello da noi proposto per superare 
tale approssimazione [IB] . 

All'interno di questa tesi e stato infine mostrato come il formalismo di 
Casida dovrebbe essere modificato per lo studio di sistemi con stato fonda- 
mentale spin polarizzato, polarizzazione non collineare, e per lo studio di una 
perturbazione di carattere piu generale di quella considerata da Casida (parte 
conclusiva del Capitolo 6). 
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Descrizione di molecole 
biatomiche 

A.l Classificazione degli orbitali molecolari 

Gli autostati elettronici di una qualsiasi molecola non possono essere classi- 
ficati sulla base del valore del momento angolare totale, in modo analogo a 
quanto viene fatto per gli atomi isolati, poiche il campo a cui sono soggetti 
gli elettoni non gode di simmetria centrale e dunque il momento angolare non 
e una grandezza conservata del sistema. Per le molecole lineari in generale e 
biatomiche in particolare e pero conservata la proiezione del momento ango- 
lare lungo l'asse che unisce i due atomi; gli stati elettronici possono dunque 
essere classificati in base a questa proiezione: m z . Nelle molecole come negli 
atomi vengono utilizzati come riferimento gli orbitali a particella singola, an- 
che quando la molecola presenta un numero superiore ad 1 di elettroni. Tale 
situazione riflette anche l'impiego di autofunzioni a particella singola nello 
schema teorico della DFT, Kohn e Sham (sistema ausiliario non interagente), 
e nello schema delle teorie di campo medio (Hartree e Hartee-Fock). 

Vengono dunque classificati come orbitali di tipo a quegli orbitali che 
hanno proiezione del momento angolare nulla, orbitali di tipo tt gli orbitali 
con proiezione m z = ±1, orbitali S quegli orbitali con proiezione m z = ±2; 
raramente vengono considerati valori maggiori della proiezione del momento 
angolare [I]. Osserviamo che tutti gli stati con momento angolare avranno 
degenerazione pari ad 1, mentre tutti gli altri stati sono doppiamente degeneri 
in assenza di campi magnetici. 
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A. 2 Classificazione delle eccitazioni molecola- 
ri 

Per quanto concerne la classificazione delle eccitazioni a molti elettroni per 
molecole biatomiche esse verrano di conseguenza indicate in base al valore 
della proiezione del momento angolare totale dello stato finale. Poiche quasi 
per la maggior parte delle molecole lo stato fondamentale presenta m z = Oil 
valore m z dello stato finale coincidera con la variazione dello stesso nell'ec- 
citazione. Un'eccitazione clie coinvolge due stati di tipo a potra essere solo 
un'eccitazione con variazione del momento angolare pari a e sara dunque 
indicata come eccitazione E, dove viene utilizzata la lettera greca maiuscola 
corrispondente alio stato con lo stesso momento angolare. Un'eccitazione che 
coinvolge uno stato a ed uno stato tt puo avere una variazione del momen- 
to angolare pari al solo valore 1 e dunque avremo un'eccitazione di tipo II. 
Possiamo cosi costruire la seguente tabella per alcune eccitazione molecolari. 



Stato 


Stato 


6m z 


tipo di 


degenerazione 


iniziale 


finale 


possibile 


eccitazione 


dell'eccitazione 


a 


a 





S 


1 


a 


TT 


1 


n 


2 


o 


5 


2 


A 


2 


TT 


a 


1 


IT 


2 


TT 


TT 


0,2 


2 E,A 


1,1,2 



Appendice B 



Differenti rappresentazioni di 
un operatore 



In meccanica quantistica una qualsiasi grandezza fisica puo essere rappresen- 
tata da un operatore auto-aggiunto che agisce sullo spazio di Hilbert delle 
funzioni d'onda. Tale operatore puo essere rappresentato in modi differenti 
a seconda che si scelga di utilizzare il formalismo della prima o della seconda 
quantizzazione e che si scelga di rappresentare tale operatore "nel punto" o 
come nella sua forma "completa" . Scopo di quest 'appendice e chiarire la nota- 
zione che abbiamo utilizzato nel presente lavoro; a questo scopo consideriamo 
l'operatore potenziale di Hartree per la TDDFT. 

In prima quantizzazione l'operatore viene rappresentato come 



mentre in seconda quantizzazione la sua rappresentazione nel puntuale e 



il che mette in modo naturale in evidenza il fatto che stiamo lavorando con 
un operatore locale, ovvero diagonale in spazio x. Le formule appena scritte 
chiariscono inoltre quale sia la relazione tra "0h( x , t) e la grandezza %(x, t) 
senza cappuccio. 

L'apice all'operatore scritto in prima quantizzazione "I" serve a distin- 
guerlo dall'operatore scritto in seconda quantizzazione. In effetti 0#(x) va 
inteso in modo differete dall'operatore £>#(x), poiche quest 'ultimo agisce sol- 
tanto nel punto x mentre vjj(x) agisce su tutto lo spazio. Altra differenza 





y,t)= (x|^[p](z,t)|y)^(x)^(y) 




(B.2) 



^[p](x,t)5 3 (x-y)^(x)^(x) 
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fondamentale e data dal fatto che l'operatore in prima quant izzazione e un 
operatore a particella singola, mentre l'operatore in seconda quantizzazione 
e un operatore a molti corpi. Ha dunque senso parlare di operatore puntuale 
e operatore completo solo in seconda quantizzazione, mentre in prima quan- 
tizzazione esiste una forma unica dello stesso, che pero puo essere utilizzata 
solo per lo studio di un singolo elettrone^. 

La forma completa dell' operatore in seconda quantizzazione h invece rap- 
presentata dalla somma sulla base degli operatori di campo della forma pun- 



Mp}(t)= / / rf 3 xd 3 y^[p](x,t)5 3 (x-y)^( x )^( x ) . (B.3) 



Lo scrivere l'operatore completo anziche semplicemente la sua versione pun- 
tuale nasce in modo del tutto naturale in seconda quantizzazione se si cerca 
di cambiare base. Per scrivere infatti un operatore su di una base qualsiasi 
non dobbiamo far altro che proiettare l'operatore su tale base. II risultato e 
il seguente: 



dove evidentemente siamo dovuti partire dall'operatore completo e non dalla 
sua versione puntuale per ottenere le componenti in spazio se fossimo 

partiti dalla sua forma puntuale avremmo ottenuto un potenziale u^[p](x, t), 
rappresentazione che non e ne in spazio (x) ne in spazio delle configurazioni. 



1 L'estensione ad un operatore a un corpo per n elettroni sarebbe dato dalla somma di 
n operatori differenti 

2 Notiamo come la variabile tcmporalc venga trattata in modo diffcrentc da quelle spa- 
ziali. La diffcrenza nasce dal fatto che il formalismo scritto abbia come punto di partcnza 
il formalismo hamiltoniano che tratta in modo differente lo spazio ed il tempo. 





MpK*) = I>M*)[p]IM^ 
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Note finali 



II lavoro della presente tesi e stato realizzato grazie al supporto che mi e 
stato dato da molte persone che fanno parte del Nanoquanta Network of 
excellence [22]. L'esperienza di quest'anno di lavoro mi ha permesso infatti 
di conoscere, almeno in parte, la realta della ricerca all'interno di differenti 
universita europee e grazie al mio relatore, il prof. Giovanni Onida, sono stato 
messo in contatto con differenti progetti di livello europeo e anche mondiale 
che sono alia base del mondo stesso della ricerca. 

Innanzi tutto il progetto Nanoquanta appunto, progetto fondato nell'am- 
bito del "Sixth Framework Programme" della Comunita Europea e che coin- 
volge dieci centri di ricerca dislocati in sette stati del'Unione Europea. Sono 
venuto in contatto con questo progetto soprattutto partecipando, in qualita 
di spettatore, al "Nanoquanta Workshop" tenutosi a Houffalize (Belgium) 
nel settembre 2006, dove ho potuto assistere alia presentazione di numerosi 
articoli e lavori di ricerca tra i quali quello dello stesso prof. Mark Casida, 
autore del formalismo che e stato studiato all'mterno di questa teso Quindi 
il progetto "ETSF" (European Theoretical Spectroscopy Facility [23]), che 
potremmo definire il futuro del progetto Nanoquanta e che si pone come 
obiettivo quello di sviluppare una serie di strumenti teorici e sperimentali, al- 
l'interno dei centri di ricerca che ne fanno parte, che permettano di descrivere 
le proprieta spettroscopiche e non solo di un qualsiasi materiale. 

Grande importanza hanno rivestito inoltre per questa tesi, cosi come la 
rivestono per la comunita scientifica, i software open sources, che vengono 
sviluppati e condivisi gratuitamenta9, tra i quali Linux [21] , che costituisce il 
sistema operativo sul quale ho costantemente lavorato e sul quale funziona in 
particolare Abinit [H] , Quest'ultimo e un software che permette di realizzare 
calcoli di struttura della materia di differente tipo e costituisce appunto il 

3 Grazie ai contatti con lo stesso Casida e stato possibile ottenere un preprint del lavoro 
che egli sta attualmente svolgendo e che abbiamo quindi utilizzato per un confronto con i 
risultati della presente tesi. 

E pensare che Bill Gates e diventato uno degli uomini piii ricchi del modo vendcndo 
un prodotto di cui esiste una valida alternativa distribuita in modo gratuito! 
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codice a cui ho in parte lavorato e che mi ha permesso di ottenere i risultati 
che sono riportati all'interno di questa tesi. In particolare ho partecipato, 
come sviluppatore, al Workshop Abinit tenutosi a Lige (Belgium) nel Gennaio 
2007 mostrando durante una breve presentazione orale le modifiche da me 
implementate nel programma e soprattutto venendo a contatto con una parte 
della comunita mondiale degli sviluppatori e degli utenti del programma. 

Un punto di riferimento importante e stato infine costituito dai professori 
Xavier Gonze, leader del progetto Abinit, e Gian-Marco Rignanese con i 
quali ho lavorato per circa cinque mesi a Louvain la Neuve (Belgium) grazie 
al sostegno del progetto Erasmus che ha permesso il mio soggiorno all'estero 
mettendo in contatto l'Universita degli Studi di Milano con L'UCL (Universit 
Catholique de Louvain), al contributo economico del progetto CNISM, che 
ha sostenuto il presente progetto di ricerca e infine, nuovamente, al prof. 
Giovanni Onida che mi ha messo in contatto con l'universita del Belgio. 
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